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PREFACE 


This is the final report on the IBM SAR-Landsat Image Registration Study performed 
under NASA contract NAS6-2827. It summarizes the effort and results of the 
entire study, and it is submitted in accordance with Item 2 of Article IX of the 
study contract. 
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Section 1 


INTRODUCTION AND SUMMARY 


1 . 1 PURPOSE 

This is the final report on IBM's SAR-Landsat Image Registration Study. It 
describes the work performed under contract NAS6-2827 and satisfies Item 2 of 
Article IX .of that contract. Some results of additional related studies outside 
the scope of this contract performed by IBM under its Independent Research and 
Development (IRAD) program are also described m this report. 


1.2 SCOPE OF THE STUDY 

The primary purpose of this study was to develop techniques and capabilities to 
assemble SAR/Landsat data sets. To accomplish this, two SAR data sets were 
extensively analyzed. Suitable algorithms and techniques were developed to 
register SAR data with corrected Landsat-MSS data. Then an example of registered 
SAR and Landsat-MSS data was produced. The result's of this processing are presented 
later in this report. 

A second purpose of this study was to contribute to the writing of the SAR/Landsat 
System Plan. This document, published by Wallops Flight Center, was edited by 
Purdue/LARS and jointly authored by Wallops Flight Center, Purdue/LARS, IBM 
Corporation, and Goodyear Aerospace Corporation. 

A third purpose of this study was to design a software system that will be able 
to produce registered SAR and Landsat data sets. If it is approved and funded by 
NASA, the system will be implemented at Purdue/LARS in FY 1979. The Landsat 
processing portion of the system was designed in detail. Documents (Appendixes A 
and B) describing this design were produced. 


1.3 SUMMARY OF EFFORT 

During the course of this study which began in April 1977, the following results 
were accomplished: 

a. One Landsat-MSS frame was geometrically corrected (using geodetic 
control points) to a 50.8-meter pixel spacing. Two subimages of this 
frame were corrected to a 25.4-meter pixel spacing. 

b. SAR data sets were registered to each of the 25.4-meter corrected 
Landsat subimages. 

c. An informal report that represents IBM's contribution to the SAR/Landsat 
System Plan was written and submitted to Purdue/LARS. 
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d. The Landsat-MSS portion of the SAR/Landsat Data Merging System was 
designed in detail. A document describing the system was written. 

e. Support of NASA planning activities was provided as required by NASA. 

An IRAD program of related investigations outside the scope of this contract was 
defined and completed. Some of the results ^are included in this report. 


1 . 4 CONCLUSIONS 

The primary conclusion to be drawn from the results of this study is that digital 
methods provide a viable technology for accurate registration of SAR and Landsat-MSS 
data. Existing software, developed previously by IBM, was used to produce an 
example of registered SAR and Landsat-MSS data. Residual RMS errors at the 
registration control points of this data set were on the order of 55 meters. 


1 . 5 RECOMMENDATIONS 

The results of this study lead to the following recommendations for future NASA 

actions: 

a. NASA should implement the SAR/Landsat Data Merging System at Purdue/LARS. 
This will provide the means by which NASA and other users can create a 
data base of registered SAR and Landsat-MSS data. 

b . Experiments should be conducted to determine how the addition of a 
registered SAR channel to the Landsat-MSS channels contributes to the 
information extraction process. 

c. NASA should fund a study to determine the differences (if any) involved 
in registering digital images from Seasat-1 SAR data with corrected 
Landsat-MSS images . Techniques for modeling geometric distortions 
should be investigated thoroughly in this study. 

d. The subject of registration of image data from different types of 
instruments (that has been shown feasible by this study) should be 
pursued with regard to other instruments. Sensors from both aircraft 
and spacecraft devices should be studied. 

1.6 NEW TECHNOLOGY 

No new technology was discovered during the course of this contract. 
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Section 2 


LANDSAT PROCESSING 


2.1 SUMMARY OF LANDSAT PROCESSING 

The first step in. the registering of aircraft SAR image data and Landsat-MSS data 
was the geometric correction of the Landsat data to a Universal Transverse Mercator 
(UTM) map projection. This included the creation of both digital and enhanced 
photographic images that were registered to a UTM map. The processing involved 
was organized into the following parts : 

a. Preprocessing 

b. Control-Point Location 

c. Scene Correction , 

d. Digital Filtering 

e. Film Plotting and Photographic Processing 

Each of these is described below or In Section 4, following a brief summary of 
Landsat-MSS digital data and its sources. 


2.2 SOURCES OF LANDSAT IMAGERY 

A multispectral scanner (MSS) is part of the payload of each of the first three 
Landsat spacecraft. In addition, inclusion of an MSS device in the Landsat-D 
payload is currently planned by NASA. The launch dates of the Landsat satellites 


follow: 
a . 

Landsat-1 , 

July 1972 

b. 

Landsat-2 , 

January 1975 

c . 

Landsat-3 , 

March 1978 

d. 

Landsat-D, 

Scheduled m 1981 


The MSS instruments on these spacecraft image the surface of the earth m four 
(Landsat-1 and Landsat-2) or five (Landsat-3 and Landsat-D) spectral bands. This 
is done simultaneously through the same optical system. The first four bands 
operate in the solar-reflected spectral region, and the fifth band operates in 
the thermal (emissive) spectral region, as follows: 


a. 

Band 

1, 

0.5 

to 

0.6 

micrometers 

b. 

Band 

2, 

0.6 

to 

0.7 

micrometers 
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c. Band 3, 0.7 to 0.8 micrometers 

d. Band 4, 0.8 to 1.1 micrometers 

e. Band 5, 10.4 to 12.6 micrometers (Iandsat-3 and Landsat-D only) 

MSS computer compatible tapes (CCTs) may be purchased from EDC. If a user does 
not know which scenes he requires, he may request a computer geographic search to 
obtain a listing of available scenes in his area of interest. To place an order, 
to inquire about the availability of data, or to establish a standing order, a 
user may contact: 

User Services Unit 

EROS Data Center 

Sioux Falls, South Dakota 57198 

Phone 605/594-6511, extension 151 

Users within the NASA community may be able to obtain MSS data directly from 
NASA-GSFC. 

A new processing system at EDC will become operational m 1978. Digital MSS data 
sets obtained prior to this are in a two-pixel-interleaved (X) format. The new 
system will produce data sets in either band-sequential (BSQ) or band-interleaved 
(BIL) format. 

The new EDC system will provide two levels of processing for MSS data: partial 

processing and full processing. Partially processed data has been radiometrically 
corrected. That is, decompression, gain, and offset adjustments have been applied 
to the image data. Fully processed data has, m addition, been geometrically 
corrected. That is, the image data has been resampled by either cubic convolution 
or nearest neighbor techniques to present the data m one of several possible map 
projections. Corrections due to spacecraft altitude, ground truth, etc., have 
been applied. 

A complete description of the Landsat program, its available products, and how to 
acquire them may be found in the Landsat Data Users Handbook, Revised . A new 
version of this document is being prepared by NASA Goddard Space Flight Center 
(GSFC) and USGS EROS Data Center (EDC) and is expected to be available in the 
latter part of 1978. 


2.3 PREPROCESSING 

Preprocessing consisted of all work of a preparatory nature that was performed 
before the major parts of the geometric correction process were started. It 
involved reorganization, evaluation, and radiometric correction of the image 
data . 

The first step was the reformatting of the digital image data from the X format 
to the BSQ format required by the existing software at IBM. There were some I/O 
errors on the input CCTs , but the reformatting was eventually completed for two 
MSS scenes: 2579-14535 (the State of Delaware and the Eastern Shore area of 

Maryland) and 2579-14541 (the southern part of the Delaware-Maryland-Virginia 
peninsula) . 
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The reformatted data was histogrammed to obtain knowledge h about the distribu- 
tion of the pixel intensities in each spectral band. The histogram program was 
run once for each spectral band of each of the two Landsat-MSS scenes. The 
output from this program consisted of absolute frequency distribution, cumulative 
frequency distribution, arithmetic mean, standard deviation, and root of the mean 
of the squares (RMS). Sample output is shown in Figure 2-1. 

Using the mean and standard deviation of the pixel intensity distribution of the 
image data and knowledge of the radiometric characteristics of the IBM Drum 
Scanner/Plotter , ^ gain and bias coefficients were computed for each spectral band 
of each scene. The data from spectral band 3 (0.7 to 0.8 micron wavelength) of 
each scene was radioraetrically adjusted using the gam and bias coefficients in a 
first-degree polynomial transformation. This radiometric adjustment was merely a 
contrast stretch that enabled the plotting device to produce a high-quality 
image. The two adjusted images were plotted on film (see Section 4.2), and the 
resulting images were visually evaluated. The pictures appeared to be of good 
quality. 

At this point, at the request of MSA, the Landsat processing was limited to the 
Delaware scene (2579-14535). All further landsat processing was done on that 
scene. 

It is well known that Landsat- 1 and Landsat-2 MSS data often exhibit an undesirable 
striping pattern that occurs on a six-scan-lme basis. This striping is especially 
evident in spectral band 1 (0.5 to 0.6 micron wavelength). A striping reduction 
program that uses a sweep mean and standard deviation equalization algorithm^ was 
used on all four spectral bands of the Deleware scene (2579-14535). A gam and 
bias correction function was used to perform this radiometric correction. A new 
gain and bias were calculated for each image-scan-lme in such a way that the 
mean and standard deviation of each corrected image scan line equaled those of 
the uncorrected six lines of that scan-mirror sweep. The resulting image data 
still displayed some striping, especially in the coastal regions. However, it 
was an improvement over the uncorrected data. The radiometrically corrected data 
was used as input to all of the remaining Landsat processing. 


2.4 CONTROL-POINT LOCATION 

In order to perform precise geometric correction of Landsat-MSS data, it was 
necessary to use scene data to define some of the error corrections. 2,3, 4, 5, 6, 7, 8, 9 
In particular, UTM map coordincates and input-image pixel coordinates of several 
control points were needed. Twenty-five control points were found. 
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Figure 2-1. Landsat-MSS Histograms (sheet 2 of 4) 
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Figure 2-1. Landsat-MSS Histograms (sheet 3 of 4) 
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Figure 2-1. Landsat-MSS Histograms (sheet 4 of 4) 
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Several years ago, IBM produced a mosaic image of the State of New Jersey. The 
bottom scene (1079-15133) in that mosaic is of the same ground area as the Delaware 
scene. An attempt was made to correlate the control points from the older scene 
to the current scene. Six control points were located in this manner, using a 
digital correlation program. Several points in the older scene either were not 
found or had a correlation peak that was too small. The remaining control points 
were found manually. 

The first step in manual control-point location was to select some potential 
control points. This was done by viewing the film image and locating features 
that were visible both on the image and on a map. Next, computer generated 
shadeprints were created for each potential control point. A shadeprint is a 
display of a subimage on a computer listing as a two-dimensional array of printed 
characters, on which gray levels are simulated by overstrike printing (see 
Figure 2-2). The shadeprints and maps were then manually compared to determine 
corresponding points in both pixel and map coordinates. Finally, the set of 
control-point locations was checked for errors by an iterative process of evalu- 
ating the geometric transformation (between corrected output space and uncorrected 
input space) derived from a given set of control points and mapping those control 
points through that transformation. Points at which the geometric errors were 
largest were rechecked for location errors, and the process was repeated for the 
refined set of points. 


2.5 SCENE CORRECTION 

Scene correction of Landsat-MSS data is the accurate geometric correction of the 
data to a UTM map projection. It involved two steps: determining the geometric 

transformation between the corrected and uncorrected images, and resampling 
(interpolating) the uncorrected data to the data samples in the corrected data 
set. 

The geometric transformation used for the Delaware scene used fifth-degree least- 
squares polynomials according to algorithms and techniques 10 developed by IBM 
under contract NAS5-21716. Twenty-five well-distributed control points were used 
to generate the image error models. The error models were evaluated at a nine-by-nine 
array of points, and the polynomial models were obtained by fitting to these 
eighty-one points. 

2 10 

Cubic convolution ’ resampling incorporating a high-frequency correction tech- 
nique was used to remove some horizontal geometric errors. Resampling was done 
at two different pixel lattices. First, the full frame was resampled to a hori- 
zontal and vertical pixel spacing of 50.8 meters. Second, two subimages (one of 
the Salisbury, MD area and one of the Cambridge, MD area) were resampled to a 
spacing of 25.4 meters. These particular spacings were chosen to give images of 
1:500,000 and 1:250,000 scales w’hen the data was plotted on the IBM Drum Scanner/ 
Plotter with a nominal pixel spacing of .004 inches in both directions. The 
Salisbury SAR data was spaced at 25 meters, and this led to the selection of 
25.4 meter spacing for the Landsat subimages. 

Copies of all digital data sets that represent processed Landsat images were sent 
to NASA/WFC or to Purdue/LARS as requested by NASA. 
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Section 3 


SAR PROCESSING 


3-1 SUMMARY OF SAR PROCESSING 

Once the landsat-MSS data was geometrically corrected to register with a UTM 
map, the SAR data was geometrically corrected to register with the corrected 
Landsat data. This included the creation of Noth digital and enhanced photographic 
images having the same geometry as the corrected Landsat images. The processing 
involved was organized into the following parts : 

a. Preprocessing 

b. Temporal Registration 

c. Digital Filtering 

d. Film Plotting and Photographic Processing. 

Each of these is described below or in Section 4. 


3.2 PREPROCESSING 

Ten aircraft SAR data sets were provided at the beginning of this contract. 

There were two different pixel spacings of each of five ground areas . All 
preprocessing done on this data was for the purpose of generating film images. 
Histograms were created for each data set. That information was used to create 
radiometrically adjusted (i.e., contrast-stretched) data suitable for film 
plotting. That data was plotted on film, and various photographic products were 
created from the resulting negatives. These ten images were of very marginal 
quality. The data was extremely noisy, and landmarks were very difficult to 
distinguish. A SAR data set containing Salisbury, MD was selected to be the 
subject of further study. 

A second SAR data set was obtained by scanning a black-and-white print of a SAR 
image of Cambridge, MD. The scanning was performed on the IBM Drum Scanner/Plotter . 
A nominally square aperture of size .004 inches was used. The pixel spacing was 
.004 inches in both the horizontal and the vertical directions. Although there 
was still some noise evident in the data, the quality of the image was considerably 
better than that of the Salisbury image. There were several lines in the image 
that appeared to be noise. They are believedto be a result of scratches in the 
negative used to produce the print from which the digital data was scanned. 

The digital data produced by the scanning process was subjected to several 
preprocessing steps to prepare it for further study. Histogrammmg, contrast 
stretching, and film plotting were done as described above to verify the quality 
of the digital data. 
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3 . 3 TEMPORAL REGISTRATION 


Temporal registration experiments were conducted with both the Salisbury and the 
Cambridge images. One version of temporally registered data for each of these 
data sets was created under IRAD funding. The Cambridge image was then subjected 
to further study and a second version of temporally registered data was created. 
The techniques used to register the SAR data were identical in each case and are 
described below. 


3.3.1 Processing Steps for Temporal Registration of SAR Image Data 

The first step was to select potential registration control points (RCPs) m 
both the corrected Landsat data and the uncorrected SAR data. This was done by 
viewing pictures made from each data set and locating features that were visible 
in both pictures. For the Salisbury image, this was very difficult to do. The 
high noise content of the SAR data resulted m few good features being visible. 
Selecting potential control points was much easier for the Cambridge image, 
since the noise content was much lower. However, there were still some regions 
of the image in which control points could not be found. 

The second step was to generate shadeprints (i.e., display of the subimage with 
a computer line printer as a two dimensional array of printed characters, on 
which gray levels are simulated by overs trike printing) for both the uncorrected 
SAR and the corrected Landsat-MSS data. One shadeprmt was generated for each 
potential control point in each of three data sets: SAR, MSS band 2, and MSS 

band 4. Some features visible in one MSS spectral band were not visible in 
another band. Bands 2 and 4 were used because they are frequently the most 
useful bands for control-point location. 

The third step was to manually determine control-point locations (i.e., corre- 
sponding pixel coordinates in each image space) by viewing the shadeprints. This 
was very difficult for the Salisbury image. The noise content rendered most of 
the SAR shadeprints useless. For the Cambridge image, this problem was not as 
severe. Consequently, many more control points were located in the Cambridge 
image. 

The fourth step can be called refinement of control-point locations. It involved 
checking the control-point locations for errors . This was done by obtaining the 
geometric transformation resulting from a set of control points, evaluating the 
geometric transformation at the control-point locations (m corrected Landsat-MSS 
space), and computing the differences between the observed and computed values 
of the control points (in uncorrected SAR space). These differences, or residual 
errors, were then studied to see if any errors in control-point location could 
be discovered. (Frequently the residual error for a mislocated point is larger 
than that of the correctly located points.) The process was iterated until a 
reasonable set of control points was obtained. 

For Salisbury, a total of 14 control points were used. For Cambridge, a total 
of 73 control points were used. In this particular example, many of these 
control points were clustered in sets of two or three similar locations. The 
important factor seems to be the distribution of the control points . Wherever 
there was a relatively large region of the SAR image in which there were no 
control points, registration accuracy observed m the photographic images seemed 
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to be lower. This was evident in the Salisbury image and in the first two 
processed versions of the Cambridge image. 

The fifth step was the calculation of the coefficients of the geometric transfor- 
mation that were used to correct the SAR data, least-squares polynomials in two 
variables of degree five were used. Each point (y,x) in corrected Landsat-MSS 
space was mapped to the point (v,u) in uncorrected SAR space. Here v = H(y,x) 
and u = V(y,x), where H and V were fifth-degree polynomials. Each polynomial 
contained 21 terms and the sum of the exponents in any term was less than or 
equal to five. (See Section 3.3.2 for a description of the least-squares algorithm.) 

The sixth, and final step was the resampling of the SAR image data. Cubic 
convolution^’ 1® resampling was used. 


3.3.2 Description of Least-Squares Fitting Algorithm 

A direct polynomial fit to the registration control points was used to define 
the transformation between corrected Landsat space and uncorrected SAR space. 

This process involved the following steps: 

a. Locating a set of RCPs in both the corrected Landsat data and the 
uncorrected SAR data, as discussed above. 

b. Selecting an appropriate polynomial form to use as a geometric transfor- 
mation. 

c. Determining the coefficients of the polynomials by a least-squares 
fit. 

The resulting polynomials were used by the resampling program to create a corrected 
SAR image. 

It can be very difficult to locate RCPs in the SAR data. The high noise content 
in the Salisbury, MD image (see EVALUATION OF TECHNICAL RESULTS section below) 
made precise location- of RCPs impossible for that image. The Cambridge, MD 
image had much lower noise content, and RCP location was much easier. The 
nature of SAR data itself contributes to the RCP location problem. The radiometry 
of a feature in an SAR image is usually totally different to that of the same 
feature in a Landsat image. In addition, shadowing or apparent widening of 
edges is present in SAR data, but not in Landsat data. 

A manual technique for RCP location was used by IBM m its experiments with the 
Salisbury and Cambridge SAR data sets. Potential RCPs were displayed using a 
computer printout in which overstrikes are used to simulate gray levels. 

(These are called shadeprrnts . ) RCPs were located by visually comparing Landsat 
and SAR shadeprints . 
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To discuss the least-squares technique used to obtain the polynomials that 
mapped corrected Landsat space into uncorrected SAR space, let 


< ‘ y i ,X i' > = llorizontal and vertical pixel coordinates of 1 th RCP in 
corrected Landsat space. 

( v i> u i) = horizontal and vertical pixel coordinates of i th RCP in 
uncorrected SAR space. 
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The columns of R will represent the coefficients of two polynomials. That 
is. 


o n / r 

H(y,x) = H. + H x + ELx + H x + H c x + W.x' 3 + 

1 4 o y 3 6 

H ? y + H g yx + Hgyx 2 + H^yx 3 + H n yx 4 + 

H 12 y2 + H 13 y2x + H 14 y2x2 + H 15 y2x3 + 

3 3 3 2 

H 16 y + H 17 y x + H 18 y x + 

H 19 y4 + H 20 y4jt + 

H 21 y5 

21 

“Z H i f i (y - x) 

i~l 

21 

V(y,x) = z v i f i (y ’ x) 
i = 1 
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The polynomials H and V approximate the mapping 


(y i > x ^) ' *- ( v i ,u i> for i=l,2,...,n 

m the least square sense. That is, 


H(y i ,x ± ) « v± 


> for i-1 ,2, . . . ,n 


V (y . , x . ) a u . 
11 i 


The coefficients of H and V are found by solving the matrix equation T = W R 
for R. The solution is 

R = ( W T W ) _1 W T T (1) 

The equations given above are for full, fifth-degree, bivariate polynomials. This 
is the form selected to correct the Cambridge image. For the Salisbury image, 
only 12 RCPs were found. Full, third-degree, bivariate polynomials were used 
for that image. 

In each case, the form of the polynomials was determined empirically. Several 
polynomials were found by performing the least-squares fit, and the resulting 
residual errors were computed. The ideal situation in doing a direct fit to 
RCPs is to have a highly overdetermmed system (equation l) and to have very low 
residual errors. This would indicate that the polynomials were good models of 
the geometric distortions. This did not happen for the two images in question. 

For the Salisbury image, the small number of RCPs led to a system that was only 
slightly overdetermined. Polynomials of degree less than that used led to unaccept- 
ably high residual errors . The Cambridge image had enough RCPs located to give a 
very overdetermined system, and the corresponding residual errors were acceptably 
small . 

Full bivariate polynomials of degrees three, four, five,- six, and seven were 
considered for the processing of the Cambridge SAR image. Those of degree five 
appeared to fit the control points best, especially near the edges of the SAR 
image data. The results of the least-squares fitting, using the final set of 
73 control points, are shown in Section 5. 
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Section 4 


IMAGE PROCESSING 


Some of the processing performed on image data under this contract was identical 
for both SAR and Landsat data. These standard image-processing techniques 
include : 

a . Enhancement 

b. Film Plotting and Photographic Processing 

\ 

c. Resampling Methods 
which are discussed below. 


4.1 ENHANCEMENT 

It is possible to modify Landsat digital data in such a way as to enhance certain 
aspects of the image. Two such aspects are contrast and edge definition. The 
purpose of both types of enhancement is to make certain features more visible in 
a photographic rendition of an image. An enhanced image ususally represents a 
less accurate reconstruction of the true data than the original Landsat digital 
radiometric data. However, enhanced images generally appear to be more pleasing 
to view and to contain more information. 


4.1.1 Contrast Enhancement 

The usual purpose of contrast enhancement is to enable an image display device 
(such as a film plotter, a CRT display, or a computer line printer) to utilize 
its full dynamic range capability. Consider a 256-level plotter being used with 
64-level Landsat-MSS data. Unless some contrast stretching is done to the data, 
the resulting image will be too dark (since only the lowest fourth of the dynamic 
range of the plotter is being used) . Another use for contrast enhancement is to 
bring out the information in a particular region in an image. 

A standard method of stretching the contrast of a digital image is to apply a 
first-degree polynomial function to the image data. That is, let 

C(i, j) = g I(i,j) + b 
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where 


th th 

I(i,j) = intensity value of i sample in j line in the input data 
g = gain coefficient, 

b = bias coefficient. 

th. t h. 

C(i,j) = intensity value of 1 sample in j line in the 
contrast-stretched data. 

It is assumed that the values C(i,j) have been clipped to fit the dynamic range 
of the display device. 

The gain and bias coefficients can be determined through a very simple statistical 
technique : 



b = M - g m 

where 

m = mean of the subimage of input data that is of interest. 

s = standard deviation of the subimage of input data. 

M = desired mean of the subimage of output data. 

S = desired standard deviation of the subimage of output data. 

The values of the statistical parameters of the subimage of output data are 
determined empirically. However, reasonable values may be estimated by letting 
M be the midpoint of the effective dynamic range of the image display device and 
by letting S be one sixth the length of the effective dynamic range of the 
device . 

For example, consider a plotter that accepts pixel values m the range 0,1,..., 255 
and that plots the same gray level for all values in the range 0,1,..., 55. In 
the range 56,57 ,... ,255 assume that the response (gray level versus pixel value) 
is approximately linear. The effective dynamic range of this hyopthetical 
plotter is 56 ,57 , . . . ,255 . Then M=155 and S=33 should be reasonable values to 
obtain a good plot with this plotter. Experience would enable these values to 
be refined. 
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4.1.2 Digital Filtering for Edge Enhancement 

Edge enhancement is an artificial sharpening of the features of a digital image. 

It is often used to accentuate boundaries of features or to bring out linear 
features of an image. Pictures of edge enhanced images are generally more 
appealing than those of unenhanced images (although this is very subjective). 

Digital edge enhancement can be performed by a digital filter (discrete convolution) 
function of the form 


3 3 

E(i,j) = E E I (i+m-2 , j+n— 2) F(m,n) 

m=l n=l 

where 

th til 

I(i,j) = intensity value of l sample in j line in the input data 

F(*,*) =3x3 matrix of filter weights 

*th. t"ii 

E(i,j) - intensity value of i sample m j line in the edge-enhanced 
data . 

The particular algorithm depends only on the matrix of filter weights 

To perform edge enhancement, the filter matrix F must define some kind of high-pass 
filter. Four examples of Laplacian-type , high-pass filters are: 


F 1 - 

0 -1 o' 

-1 5 -1 

F 2 = 

"-1 -2 -1 
-2 13 -2 


1 

o 

1 

o 


-l -2 - 1 . 

F 3 = 

- - 

-1 -1 -1 

-1 9 -1 

F 4 

i -2 r 

-2 5 -2 


-1 -1 -1 _ 


, 1-2 1 _ 


Each of these examples has the property that the sum of the weights is one. 

This results in an image that is the sum or the original image plus a Laplacian- 
filtered image. Therefore, the resulting radiometry of an edge-enhanced image 
will be somewhat comparable (statistically) to that of the original image. 
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4.1.3 Enhancement Processing 

The Salisbury subimage, corrected to 25.4-meter pixel spacing, was digitally 
filtered for edge enhancement in two ways. First, a Laplacian filter was used on 
the data and the result added to the image data. Numerically, this amounted to 
multiplying a given pixel value by five and subtracting each of its four (horizontal 
and vertical) neighboring pixel values (filter F^ above). The resulting image 
contained enhanced edges and also contained the intensity distribution of the 
original data. 

The second method involved using a Laplacian filter on the data and adding a 
constant to the result. This amounted to multiplying a given pixel value by four 
and subtracting each of its four neighboring pixel values. The resulting image 
contained only edge information. 

In addition, contrast enhancement (as described in Section 4.1.1) was performed 
on all Landsat and SAR images that were created under this contract. The purpose 
of this processing was to obtain digital data that was compatible with and well-suited 
to the film plotting device that was used to make photographic negatives. 


4.2 FILM PLOTTING AND PHOTOGRAPHIC PROCESSING 

Digital images were converted to film negatives on the IBM Drum Scanner/Plotter 1 . 
From these negatives, two basic photographic products were made: black-and-white 

contact prints, and false-color contact prints. The false-color prints used MSS 
spectral bands 1, 2, and 4 to give yellow, magenta, and cyan colors on the prints 

Copies of all photographic products have been delivered to Dr. Harold Maurer 
(NASA/WFC). Some pictures were also sent to Paul Anuta (Purdue/LARS) for his use 
in a related effort. 

In addition, 35 mm slides and other transparencies were prepared and sent to 
Dr. Harold Maurer as requested. 


4.3 RESAMPLING METHODS 

One of the requirements normally imposed upon image processing when the image is 
m digital form is that the output pixel lattice be regular, or equispaced, in 
the output space. This is needed for two reasons: compact data storage, and the 

limitations of many film recorders. If the output pixel lattice is not regular, 
then additional information must be carried to specify pixel location. The last 
requirement precludes the simple repositioning of the original image pixels. The 
only way to change geometry and specify the target pixel location, too, is to 
"resample" the orignal image. Resampling here consists of calculating the location 
of a particular target pixel in the original image and interpolating over the 
surrounding original pixels to determine the output intensity. During processing, 
the target pixels are considered output and the original pixels are considered 
input. Figure 4-1 shows the relationship between input and output pixels. 

In the following discussion, the calculation of the output pixel location in the 
input image is presumed. 
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INPUT 


OUTPUT 



Figure 4-1. Locating Output Points in Input Space 


Cubic convolution is a family of resamplers that approximate the SINC, or SINE (x)/x 
function. The SINC function is the theoretical perfect resampler, but it requires 
an infinite number of terms. Cubic convolution substitutes a truncated piecewise 
cubic approximation to the SINC function, so that the amount of processing is feasi- 
ble for large image data sets. Cubic convolution is a one-dimensional process that 
must be repeated to provide for two-dimensional resampling. This is illustrated in 
Figure 4-2, where four horizontal resamplings provide the values for a final vertical 
resampling. The final vertical resampling then provides the output intensity. This 
is for a so-called four-point resampler. 



Figure 4-2. Two-Way Cubic Convolution Resampling 


The equation of the classic four-point cubic- convolution resampling algorithm 
used for this contract is 

i out = + ^-r^vw] + th-i.p} + L (2) 

An experimental program that evaluates this equation five times for each output 
pixel and that uses 32-bit floating point arithemetic was used to resample both 
Landsat and SAR data. 
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For SAR data, each of the four one -dimensional horizontal resamplings used the 
same value for D in Equation 1. That is, D = as m Figure 4-2. For Landsat-MSS 
data, each one-dimensional resampling used a distinct value for D. This is due 
to the use of the resampling program to remove high-frequency geometric errors 
that occur m the horizontal direction m MSS data. 
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Section 5 


TECHNICAL EVALUATIONS 


In this section, the evaluations of the geometry of the corrected Landsat-MSS 
image and the corrected Salisbury SAR image, and the final version of the corrected 
Cambridge SAR image are described. 


5.1 LANDSAT PROCESSING EVALUATION 

Under contract NAS6-2827, IBM processed Landsat-2 MSS scene number 2579-14535. 

The processing performed was scene correction. That is, the image was geometric- 
ally corrected to a Universal Transverse Mercator (UTM) map projection using 
geodetic control points to model scene dependent errors (attitude and altitude) . 
Digital data at two scales was produced. First, the entire scene was processed, 
resulting in an image with pixels at a spacing of 50.8 meters in both the horizontal 
and the vertical directions. Second, two subimages of the scene were processed 
with pixels spaced at 25.4 meters in both directions. 

The geometric accuracy of the corrected Landsat images was evaluated by checking 
the errors at the locations of the geodetic control points (GCPs) that were used 
to correct the image. For each GCP, a nominal input-space location was found by 
the following transformations: 

a. Mapping the UTM coordinates into output-space pixel coordinates using 
the elementary rotation, translation, and scale-change transformation 
which exists between the two coordinate systems 

b. Mapping the output-space coordinates into input-space pixel coordinates 
using the mapping polynomials that define the geometric transformation 
for the scene. 

These nominal input-space locations were compared with the corresponding observed 
input-space locations that were found manually from the input data. The resulting 
errors m meters were found. The RMS of these errors, as shown in Table 5-1, was 
77.3 meters. 

The corrected digital data was recorded on film. The resulting photographic 
image products were of good quality. 


5.2 SAR IMAGE PROCESSING EVALUATION 


5.2.1 Salisbury SAR Image 

A subset of the SAR image containing Salisbury, MD was chosen to be temporally 
registered to the corrected Landsat data from MSS scene 2579-14535. This SAR 
data had very high noise content which hindered the finding of registration 
control-point (RCP) locations. Due to the sparsity of suitable visible features 
in the SAR data, only 14 RCPs were located. 
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Table 5-1* Results of Accuracy Check for Landsat Data 


Ln 

i 

N> 




Measured 

UTM Map Coordinates 


Observed 

Input-Space 

Coordinates 

Nominal 

Input-Space 

Coordinates 

Input-Space 
Differences 
In Pixels 

Radial 
Error In 
Meters 

GGP- 

NORTHING- 

•EASTING 

ZONE 

YOB& 

- XG6S 

YNOM 

- XNOM- 

DELY- 

DELX- 

R5SM- 

1 

4368380 

599200 

18 

361.00 

729.00 

359.69 

730.32 

1 .3081 

- 1 .3200 

1 28 .5 1 7 

2 

4319160 

370820 

18 

91.00 

1 4 02 . 00 

90.57 

1402.00 

0.4348 

-0.0004 

24.843 

3 

4264533 

450858 

ie- 

1618.00 

1 623.00 

1618.81 

1622.84 

-0.8133 

0.1637 

48.244 

4 

4281320 

465530 

18 

13b3.0O 

1 6 25. OO 

1883. lO 

1 625.06 

-0 .0990- 

—0.0551 

7.1-4S- 

5 

4265817 

477486 

18 

2155.00 

l 764. 00 

2156.48 

1764.25 

-1 .4 757 

-0.2456 

86 .5 37 

6 

4269440 

407370 

18 

945.00 

1 9 15.00 

944.69 

1917.11 

0.3063 

-2. 1149 

1 63 .4 1 3 

7- 

4249186 

379330 

-18 

554 .00 

2236.00 

554 .46 

2236.31 

-0 .4646 

-0.3072 

36.014 

8 

4235348 

41 9500 

18 

1303.00 

2304 .00 

1303.69 

2303.83 

-0.6935 

0.1690 

41 .826 

9 

4238490 

420600 

18 

1308 .00 

2 2 62.00 

1308.49 

2262.53 

-0 .4887 

-0.5264 

50 .1 8 1 

10 

424321 0 

370290 

18 

425.00 

2333.00 

-425 .45 

2332.68 

— O *4 459 

9.31 72- 

-35.783- 

It 

424 1 985 

375 72-0 

1-6 

524.00 

2 3 33.00 

524.59 

2333.89 

-0.5905 

—0. 8859 

77.8 53 

12 

4390075 

399 120 

18 

261.00 

46S. 00 

260.61 

465.84 

0.3899 

-0.8435 

70 .425 

13 

43«2940 

408410 

18 

454.00 

529.00 

453.25 

529. SI 

0.7472 

-0.51 34 

58.959 

14 

4377045 

413810 

18 

573.00 

-be6.oo 

572 . 66 

587.82- 

-0.-337-3- 

— 1.-82 18 

-1 45 .569 

-IS 

-4350880 

38 7920 

16 

244*00 

971 . 00 

243.51 

972. 12 

0 .4946 

-1.1209 

93 . 1 67 

16 

4378540 

449 91 0 

18 

1185.00 

4 78. 00 

1183.88 

478.70 

1.1 172 

-0.6984 

84 .4 7 0 

17 

4381505 

451450 

18 

1196.00 

4 36.00 

1 197.01 

438.67 

-1 .0054 

-0.6745 

78.4 52 

18 

4372495 

453 54 0 

18 

1273.00- 

543.00 

-12 72 .57 

543. 26 

0 .4295 

-0.2598 

-32 .022- 

19 

4367605 

453665 

16 

1298.00 

599.00 

1298.92 

599.61 

-0.9172 

-0.61 15 

71 .356 

20 

4379950 

425480 

18 

760.00 

523.00 

760 .06 

522.99 

-0.0573 

0.0066 

3.314 

21 

4373880 

422460 

18 

736.0 0 

6 04 . 00 

735.36 

604.63 

0.6420 

-0.6303 

61 .952 

22 

4357310 

493810 

18 

2026.00 

626.00 

2025.99 

626.-98 

0.0 14-fr 

—0.9799 

77.609 

23 

43 19348- 

516 745 

18 

2589.00 

1 0 31 .00 

2588 .68 

1031.95 

0.3244 

-0.9549 

77.866 

24 

43 17520 

61 1 74 0 

18 

2510.00 

1 0 67 . 00 

2510.85 

1 066.76 

-0 .8480 

0.2375 

61 .979 

2b 

4273070 

494460 

18 

2414.00 

1652.00 2414.46 1652.81 

MEAN DF EACH COLUMN t 
RMS OF RSSM : 

-0 .4622 
-0.0726 

-0.8053 

-0.5790 

69 .0 32 

67.261 

77.301 



Full, bivariate, cubic, least-squares polynomials were used to define the trans- 
formation between corrected and uncorrected SAR data. The RSS residual errors 
(in pixels) of the least-squares fit at the RCP locations are shown in Table 5-2. 
These results seem reasonable. 

The resulting image that was produced using the least-squares fit summarized in 
Table 5-2 had reasonable registration at the 14 RCP locations. However, some 
severe misregistration is evident at locations between the RCPs For example, US 
Route 50 east of Salisbury is clearly visible in both the SAR and the Landsat 
data. In the corrected SAR image, that road is displaced to the north of its 
proper location. RCPs above (lake) and below (airport) the misregistered road 
are properly registered. 

The conclusion that can be drawn is that the distortions m the SAR data are too 
severe to be modeled by cubic polynomials. Since the maximum number of terms in 
a least-squares polynomial is directly related to the number of points being 
fitted, it would be difficult to obtain significant improvement in the correction 
accuracy of the Salisbury SAR image. 


5.2.2 Cambridge SAR Image 

Exactly the same procedure as was used on the Salisbury image was used on the 
Cambridge image. However, 73 control points and fifth-degree polynomials were 
used for the least-squares fitting process. The results from the least-squares 
fit are shown in Table 5-3. The mean RSS error of about 5 pixels is quite accept- 
able, since the pixels in the Cambridge image are much smaller (about 10 meters) 
than were those m the Salisbury image. 

A visual check of false-color composite photographic images of this SAR image 
merged with corrected Landsat-MSS data indicated that the SAR data were suitably 
registered to the Landsat data. The visual check is evidently necessary. The 
first two attempts at correcting the Cambridge SAR image used some incorrectly 
located control points, but a similar mean RSS error of about 5 pixels at the 
control-point locations was attained in each case. Each of these first two 
versions was easily found to be misregistered by visual inspection of false-color 
composite pictures. 

It can also be seen from the photographic products that the SAR data added infor- 
mation that was not present in the Landsat data. For example, highly detailed 
information in the Cambridge urban area was present in the SAR data, but not in 
the Landsat data. 
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Table 5-2. Results of Accuracy Check for Salisbury SAR Data 
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Table 5-3. Results of Accuracy Check for Final Cambridge SAR Data 


LSA5T-5QUARE5 FIT EVALUATION 8V SUBROUTINE LS22EVS 


RESULTS OF LEAST SQUARES I 

Observed 

Output-Space Coordinates 
300*0000 246.0000 

339.000 0 261.0000 

235.0000 290.0000 
856 *0000 
8TI *0000 

606.0000 

264 .0000 

660.0000 
493.0000 


750.0000 

359.0000 

496.000 0 

531.0000 

575 .0000 

633.0000 
701 .0000 

293. 0000 

339.0000 

387.0000 


1039.0000 

1015.0000 

1096.0000 

1 1 12.0000 

1 147.0000 
1 X 16.0000 



743.0000 

833.0000 

912.0000 

1 020.0000 

889.0000 

397.0000 

1037.0000 

1 1 78.0000 


Observed 

Input-Space Coordinates 

132.0000 


254.000 0 

369.0000 
1 14.0000 

1696.0000 

1795.0000 

1229.0000 

534.0000 

1254.0000 

1062.0000 

1797.0000 

1779.0000 
408 
717 
830 
957 

1120 
1271 .000 
302 
4 20 
579 

1659.0000 

1731 .0000 

1222.0000 

951 .OOOo 

1 1 64 
621 
102 


1073.0000 

1165.0000 

1283.0000 

1541 .0000 

1513.0000 

1009.0000 


1920.0000 

2016.0000 
2059.0000 


1410.0000 

1^81 .0000 

1601 .0000 

1817.0000 
82.0000 

588 .0000 

683.0000 

232.0000 

1501 .0000 

1493.0000 
1864 .OOOO 

583.0000 
274 .OOOO 

282.0000 
1 33.0000 
451 .OOOO 

38.0000 

1253.0000 


94.0000 

127.0000 

1299.0000 

1670.0000 
1 769.0000 

1723.0000 

1270.0000 

1069.0000 

227.0000 

1337.0000 


2334.0000 

2341.0000 

2439.0000 
2292 .DOOO 

800.0000 


1870 .DOOO 
1851 .0000 
1660.0000 
1298.0000 
76.0000 

308.0000 

629.0000 

031.0000 
I 030.0000 


1329.0000 

1506.0000 

1713.0000 

2038.0000 

1760.0000 

1840.0000 

2454.0000 

2452.0000 
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Input-Space Coordinates 
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6.0857 5.6744 

6.99)3 1.6406 

-2*9369 3.1659 

-4*4091 — 0.1377 

S . 2485 4.3831 

4.5510 4.6211 

-6*2489 -0.8225 

-3.6395 -4.8179 

-1.0979 4.3368 

6.158b -7.9611 

3.3519 4.3629 


ID NUMBER or POINT KITH LAWCEST RESIOlML EHROM 


RETURNING NORMALLY FROM SUHROUTIH*. L*^2EV 
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Section 6 


SAR/LANDSAT PROCESSING CAPABILITY 


tinder this contract, a SAR/Landsat Data Merging System (SLDMS) was designed. The 
documentation produced during this effort is included in the appendixes of this 
report: 

a. Appendix A - Functional Description of the SAR/Iandsat Data Merging 
System 

b. Appendix B - Detailed Design of SAR/Landsat Data Merging System 
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Section. 7 


SUPPORT TO NASA PLANNING ACTIVITIES 


7 . 1 DOCUMENTATION 

IBM generated documentation, as required by NASA, to support its planning activi- 
ties. In particular, inputs were provided to support the 1978 RTOP on SAR/Landsat 
Data Merging. 

i 

7.2 PHOTOGRAPHIC PRODUCTS FOR PRESENTATIONS 

At various times during the course of this contract, IBM generated photographic 
products of the SAR and Landsat-MSS data being investigated. Some of these 
products were produced in color transparency, 35 mm slide, or lantern slide 
format for use as visual aids for oral presentations by NASA investigators. 
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Section 8 


RESULTS AND CONCLUSIONS 


8.1 FEASIBILITY OF SAR/LANDSAT REGISTRATION 

The results from the experiment on the Cambridge SAR image indicate that temporal 
registration of SAR data with Landsat-MSS data is both feasible (from a technical 
viewpoint) and useful (from an information- content viewpoint) . Multi-sensor 
temporal registration is likely to be a very desirable process in the future. It 
will enable users to make use of varied sources of information to solve their 
problems. This experiment shows the viability of using SAR data m this way. 


8.2 CONTROL-POINT LOCATION RESULTS 

This investigation confirmed that the greatest difficulty in registering aircraft 
SAR data to corrected Landsat-MSS data is control-point location. The difficulties 
occur m the following areas: 

a. Feature selection 

b. Distribution of control points 

c. Comparing Landsat and SAR features 

d. Correctness of control-point locations. 

The results and conclusions regarding control-point location that are implied by 
this study are described below. 

The selection of features that will serve as good control points is made difficult 
by the differences in SAR and MSS data . Features that are usually high-quality 
control points in the MSS data (such as road intersections) were often unusable in 
the SAR data. Features of normally lower quality in MSS data (such as agricul- 
tural field boundaries) had to be used. This problem is compounded by the fact 
that features often looked quite different in the data from the two sensors . For 
example, the shadowing in the SAR data did not occur in the MSS data. 

The experiments on the Salisbury and Cambridge SAR images clearly demonstrated 
the importance of having a set of control points that is well-distributed throughout 
the image. In every case tried during this contract (and the related IRAD work), 
there were geometric errors in any relatively large areas that contained no 
control points. 

The experience on this contract indicates that SAR and Landsat-MSS data are 
unsuitable for automatic computer correlation of digital control-point data. It 
is clear that the gray-level data cannot be compared by computer, due to the 
different response characteristics of the MSS and SAR images . The image data 
could be processed by an edge-detection algorithm prior to machine correlation, 
but the SAR shadowing would likely cause this method to be ineffective. However, 
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the high noise content of the SAR data makes edge detection very difficult, if 
not impossible. The rotational difference between the SAR and Landsat-MSS data 
precluded any attempt at machine correlation of control points during this study. 

All of the problems normally associated with checking for errors in control-point 
location are present in the SAR/Landsat registration problem. However, extreme 
care must be taken to avoid obvious errors. All points near the edges should be 

double checked, since errors in these points can easily go undetected. Errors in. 

control-point coordinates caused the first two versions of the corrected Cambridge 
SAR data to have severe geometric errors. Most of these control-point location 
errors were on points near the border of the SAR image. When a set of control 

points that contained only correctly-located points was used, the resulting corrected 

image had good registration with the corrected Landsat-MSS image. 
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A. 1 GENERAL SYSTEM REQUIREMENTS 


A. 1.1 Functional 

The SAR/Landsat Data Merging System (SLDMS) will have the capability of pro- 
ducing SAR and Landsat MSS images that are registered to each other and to a 
Universal Transverse Mercator (UTM) map projection. A SLDMS user will accom- 
plish this by processing the Landsat data to his own specifications and then 
processing the SAR data to register with the corrected Landsat data. He then 
can use the corrected SAR and Landsat data for multispectral classification and 
other information extraction processes. 

The specific Landsat processing functions that the SLDMS will have are radio- 
metric correction and geometric correction. Radiometric correction is defined 
to be a striping reduction process that is performed on MSS data that has been 
radiometrically calibrated but not resampled. Geometric correction is defined 
to be the resampling of MSS data to a UTM map projection. The user will be able 
to specify both the spacing between pixels in the output UTM image and the 
orientation of the image (i.e., the direction of north relative to the image 
scan lines) . 

The specific SAR processing function that the SLDMS will have is registration. 
Registration is defined to be the resampling of SAR data to have the same 
geometry as a corrected Landsat-MSS data set. That is, the SAR data will be 
resampled to the same lattice of pixels as was the geometrically corrected 
Landsat-MSS data. 


A. 1.2 Hardware 


The SLDMS will run on the IBM 370/148 computer at LARS. The only hardware 
requirements for the system are the tape drives needed for image data sets and 
sufficient virtual memory for the resampling program. The system will be 
designed to minimize the use of disk space by using tapes wherever it is reason- 
able to do so. 

The resampling program will require an internal input buffer large enough to 
store the maximum number of input image lines that will be needed to create one 
output image line. This is a function of the rotational difference between the 
input and output images. Therefore, the memory requirement varies from image to 
image . 


A. 1.3 Software 

The SLDMS will be written in Fortran IV and IBM 370 Assembler Language. Gen- 
erally, Fortran will be used wherever it is reasonable to do so. However, 
Assembler Language will be used for those programs that would be significantly 
more inefficient if written in Fortran. 

The programs will use some local LARS programs to perform certain standard 
functions. For example, the program TAPOP will be used to perform image data 
set I/O. 
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A. 2 IMPLEMENTATION AT LARS 


A. 2.1 Landsat-MSS 

Landsat-MSS processing software which will be implemented at LARS will be 
defined in this section. Full use of techniques, algorithms, and software 
developed under the Master Data Processor (NAS5-22999) and the Landsat-C Return 
Beam Videcon (RBV) Software and Interactive Ground Control Point (GCP) System 
(NAS5-23790) contracts will be made in implementing the SLDMS. In particular, 
the Control Point Library Copy Tape, which is a back-up copy of NASA’s control 
point library, will be used by the SLDMS. 

The Landsat-MSS software will consist of five programs: 

a . Reformatting Program 

b. Automatic Control-Point Location Program 

c. Manual Control-Point location Program 

d. Geometric Transformation Program 

e. Resampling Program. 

The inputs, functions, and outputs of each of these programs are described 
below. 

There are three kinds of input MSS data that can be processed by the SLDMS: 

a. Uncorrected Data -- This consists of four bands of MSS data in the X 
format. The image data has been radiometrically corrected and line- 
length adjusted, but no other corrections have been applied. Ancil- 
lary data includes no geometric transformation. 

b. Partially Processed Data — This consists of four or five bands of MSS 
data in the BSQ format. The only correction applied to the image data 
is radiometric correction. Ancillary data includes geometric trans- 
formation information. 

c. Fully Corrected Data -- This consists of four or five bands of MSS 
data in the BSQ format. Both radiometric correction and geometric 
correction have been applied to the data. 

The SLDMS will be able to geometrically correct and perform a striping reduction 
on both uncorrected and partially processed data. It will also be able to 
'resample fully corrected data in order to change the pixel spacing. 

As shown m Figure A-l, there are three paths through the software. The first 
path involves the reformatting, automatic control-point location, geometric 
transformation, and resampling programs. This path would normally be used for 
an uncorrected MSS scene for which a corresponding control-point library exists. 
The result from this path is a fully corrected <MSS scene with user-selected 
orientation and pixel spacing. This would be the most desirable path for pro- 
cessing uncorrected data. It can also be used to correct partially processed 
data. If no control-point library exists for an uncorrected MSS scene that is 
to be processed, there are two alternate paths through the system that can be 
used. 




Figure A-l. Software Paths 
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The second path involves the reformatting, geometric transformation, and resamp- 
ling programs. This path would normally be used for an uncorrected MSS scene 
for which no control-point library exists or for a partially processed MSS 
scene. The result from using this path on an uncorrected MSS scene would be a 
systematically corrected MSS scene with user-selected orientation and pixel 
spacing. From a partially processed scene, a fully corrected MSS scene that has 
been either systematically corrected or scene corrected (depending on the ori- 
ginal ancillary data) can be produced with this path through the software. 

The third path involves the reformatting, manual control-point location, geo- 
metric correction, and resampling programs. This path would normally he used 
for an uncorrected MSS scene for which no corresponding control-point library 
exists or for a fully corrected MSS scene. The result from this path would be a 
fully corrected MSS scene in either case. For a fully corrected input scene, 
only the reformatting and the resampling portion of the manual control-point 
location program are used. In this case, the result is a fully corrected scene 
with different pixel spacing. 


A. 2. 2 Aircraft SAR 

The aircraft SAR processing software which will be implemented at LARS will be 
defined in this section. The techniques and algorithms were developed by IBM 
under contract NAS6-2827. 

The aircraft SAR software will consist of three programs: 

a. Manual Control-Point Location Program 

b. Geometric Transformation Program 

c. Resampling Program. 

These programs have the same functions as the corresponding Landsat programs. 

The Geometric Transformation Program is a least-squares fitting program. It 
will be a modification of an existing LARS program. This modification will be 
performed by LARS, and it is not included m this document. 

The other two programs are minor modifications of the corresponding Landsat 
programs. That is, the MSS high-frequency corrections must be removed in order 
to process SAR data. These are included m this document. 

The input to this portion of the SLDMS consists of digital image data from 
aircraft SAR instruments. The first step is the location of control points m 
both the SAR data and the corrected Landsat data. Then a geometric transfor- 
mation that maps points (X,Y) in the corrected Landsat image to points (U,V) in 
the SAR image is obtained in the form of 21-term fifth-degree polynomials u(X,Y) 
and v(X,Y) such that U = u(X,Y) and V - v(X,Y). Finally, the SAR data is 
resampled using the above geometric transformation. The result is a corrected 
SAR data set that is registered to the corrected Landsat data set. 
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A. 2. 3 Seasat SAR 


The only difference between the Seasat and aircraft SAR processing software is 
the Geometric Transformation Program. The Manual Control-Point Location Pro- 
grams and the Resampling Program used for aircraft SAR data will also be usable 
for Seasat SAR data (assuming that it is 8-bxt data). 

The Geometric Transformation Program will use the algorithms defined in the 
Seasat SAR Investigation described above. 

The input to the Seasat SAR portion of the SLDMS is digital image data from the 
Seasat SAR instrument and the corresponding Seasat ancillary parameters. The 
output is a corrected SAR data set that is registered to the corrected Landsat 
data set. 


A. 3 USER CONSIDERATIONS 

Users of the SAR/Landsat Data Merging System (SLDMS) have several choices to 
make when processing a Landsat data set. These user options are discussed in 
this section. 

In some cases, a user will have a choice of using partially processed or fully 
corrected input MSS data. Since fully corrected data is already resampled to a 
standard map projection, further resampling may not be necessary. A common 
reason for further resampling would be the need to have the corrected Landsat 
data at a different pixel spacing than the 57 meters (horizontal and vertical) 
of the fully corrected data set. Although resampling the data only once is 
clearly desirable, a recent study (Bi-resampled Data Study; Final Report for 
Contract NAS5-23708; R. Benner, W. Young; IBM Corp., March 1977) has indicated 
that a second resampling will degrade the data only slightly andwill not 
essentially change multispectral classification results. 

When processing uncorrected MSS data and some partially processed MSS data, a 
SLDMS user must decide whether to use geodetic control points in determining the 
geometric transformation between the corrected and uncorrected spaces. If a 
systematic correction will provide satisfactory geometric accuracy to meet a 
user's needs and if no control-point library exists for a particular scene, then 
the expense of control-point location can be avoided by using systematic cor- 
rection. 

The SLDMS will have the capability to resample with nearest neighbor or cubic 
convolution. Although cubic convolution is widely regarded as a better algo- 
rithm than nearest neighbor, it is considerably more expensive to perform on a 
general purpose computer. Computer expense for resampling is directly propor- 
tional to the number of pixels that are being created. 

Pixel spacing m the corrected Landsat data is a user option in the SLDMS. MSS 
input data is spaced at about 57 meters horizontally and at about 79 meters 
vertically. Fully corrected MSS data produced by the MDP will be spaced at 57 
meters in both directions. A particular spacing may be chosen because the SAR 
data has a similar spacing. A second factor that a user must consider when 
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choosing the pixel spacing is resampling cost. The number of output pixels that 
are created in the resampling step is inversely proportional to the pixel spac- 
ing. 

Another user option for Landsat MSS processing is the orientation of the image. 
Orientation is the angle between North and an image scan line at some point m 
the image. The standard orientation is such that the spacecraft velocity vector 
is approximately perpendicular to the image scan lines. That is, input-space 
scan lines are nearly parallel to output-space scan lines. Orientation affects 
the resampling step. The memory required for the resampling program is directly 
proportional to the angle between the input-space scan lines and the output- 
space scan lines. 


A. 4 FUNCTIONAL REQUIREMENTS OF SLDMS PROGRAMS 


A. 4.1 Reformatting Program 

The reformatting program converts Landsat MSS computer compatible tapes (CCTs) 
into LARS Multispectral Image Storage Tape (MIST) format. The inputs, func- 
tions, and outputs of this program are stated in Table A-2. It is important to 
note that the reformatting program performs different functions on each of the 
different kinds of input data. The different possibilities are shown in 
Table A-l . 


Table A-l. Reformatting Program Functions 


Uncorrected 

Partially Processed 

Fully Corrected 

Convert to MIST 

Convert to MIST 

Convert to MIST 

Remove LCC pixels 



Detector histograms 

Detector histograms 


Striping red. tables 

Striping red. tables 


Striping reduction 

Striping reduction 


Extract ancillary 

Extract ancillary 

Extract ancillary 


A. 4.2 Automatic Control-Point Location Program 

The automatic control-point location program is used to locate the input-space 
coordinates of a set of points whose geodetic coordinates are known. The output 
from this program is a set of pixel coordinates of control points. These 
coordinates are specified to a fraction of a pixel, and the location accuracy is 
assumed to be within one tenth of a pixel for a good control point. 

The specific functions of this program are shown m Table A-3. Basically, the 
program extracts suitable prototype control points from a tape, creates cor- 
responding search areas from an image tape, correlates the data, and determines 
the search-area coordinates of the control points. The tape containing the 
control points is a Control Point Library Copy Tape, which is a back-up copy of 
NASA's control-point library. It is also possible to obtain prototype control 
points from an image tape. 
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Table A-2. Reformatting Program Functions 


Input 


Process 


1. MSS 

computer tapes • 

1. 

Convert uncorrected MSS data 
from X format to MIST format. 

a . 

Uncorrected 

2. 

Convert partially processed 

b. 

Partially processed 


MSS data from BSQ format to 
MIST format. 

c. 

Fully corrected 

3. 

Convert fully corrected MSS 

2. User's requests. 


data from BSQ format to MIST 




format. 



4. 

Remove line-length correction 
(LLC) pixels from uncorrected 




data . 



5. 

Compute detector histograms . 



6 . 

Compute striping reduction 
tables . 



7. 

Perform striping reduction. 



8. 

Extract ancillary data from 
input tapes . 


Output 

1. MIST format tape. 

2. Ancillary data set. 

3. Printer listing. 
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4. 

5. 

6 . 


Program Functions 


Process 

Output 

Extract prototype control 

1. Enlarged ancillary data set 

points from: 

containing control point 


locations . 

a. Control Point Library 


Copy Tape 

2. Printer listing. 

b. MIST format image tape 



Extract search areas from 
image tape. 

Horizontally resample search 
areas (and control points, if 
from image tape) to correct 
high-frequency errors . 

Obtain correlation surface 
for each control point. 

Fit correlation surface to 
find control point location 
to sub-pixel accuracy. 

Evaluate quality of each 
control point. 


A. 4. 3 Manual Control-Point Location Program 

The basic function of the manual control point location program is to process 
image data in a way that will enable a user to determine control-point coordin- 
ates manually. The specific functions of this program are listed in Table A-4. 

In order to assist a user to determine control-point coordinates manually, this 
program produces two kinds of output. First, it generates a resampled image 
data set for each control point. Second, it displays a control-point data set 
on a computer listing by simulating pixel gray levels with printer overstrikes. 

The resampling that is performed for a control-point location task has two 
purposes. First, horizontal geometric distortions that are high-frequency in 
nature are removed during the resampling process. These distortions are listed 
in Table A-5. Second, a change of scale is performed. This is done to enlarge 
the control-point subimage so that a user may determine the control-point 
coordinates to a sub-pixel accuracy level. 


A. 4. 4 Geometric Transformation Program 

The geometric transformation program has the- following basic functions for 
uncorrected or partially processed scenes: 

a. To define the transformation between the corrected output space and 
the high-frequency-corrected input space. 

b. To assess the accuracy of the geometric transformation. 

The resulting geometric transformation is added to the ancillary data set that 
corresponds to the scene. It is specified as a finite, tabular function. A set 
of output-space grid points is defined, and the corrected input-space coordin- 
ates are determined. This grid-point correspondence is the primary output of 
the program. 

Table A-4. High-Frequency Horizontal Geometric Errors 


Line length variations 
Earth rotation errors 
Sampling delay errors 
Band-to-band offset 
Mirror velocity errors 


The specific functions of this program are shown in Table A-5 . For uncorrected 
data for which no control points exist, models that correct all known systematic 
errors are used. If control point locations are available for an uncorrected or 
a partially processed data set, models that correct all known systematic and 
scene-dependent errors are used. These error models are then used to determine 
grid point correspondence. 




A- 10 


Table A-5. Manual Control-Point Location Program 


Input 

Process 


Output 

1. MIST format tape (from 

1 . Extract small subimages to 

1 . 

MIST format tape. 

Reformatting) . 

computer memory. 

2. 

Shadeprints on listing data 

2. User's requests. 

2. Resample small subimages 


set . 


including: 

3. 

Printer listing. 


a . Scale change 


b. High-frequency error 
correction 


3. Histogram small subimages. 

4. Obtain translate tables 
using histograms. 

5. Print shadeprint on listing 
with overstrikes. 



In the case of scene correction, this program also provides an assessment of the 
geometric errors. This is done using the covariance matrix of the attitude/ 
altitude fitting process, as described m section LANDSAT MSS IMAGE TO MAP 
REGISTRATION. 

This program can also create a grid point correspondence from the ancillary data 
for a partially processed scene. 


A. 4. 5. Resampling Program 

The resampling program creates a geometrically corrected image data set from an 
input image data set. It can resample the data using nearest-neighbor or cubic- 
convolution resampling algorithms. 

The resampling program can process the following kinds of data: 

a. Uncorrected or partially processed Landsat MSS data. 

b. Fully processed landsat MSS data. 

c. SAR image data. 

The program operates in essentially the same way for all three kinds of data. 
However, when uncorrected or partially processed Landsat MSS data is resampled, 
the program removes horizontal, high-frequency, geometric errors as part of the 
resampling process. 

The specific inputs, functions, and outputs are shown m Table A -6. A check- 
point-restart capability is included because cubic-convolution resampling is a 
computationally bound process. 
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Table A-6. Geometric Transformation Program 


Input 


Process 

Output 

1. Ancillary data set. 

2. 

Systematic error modeling 

1. Enlarged ancillary data set 



for uncorrected MSS data. 

containing grind-pomt 

2. User's requests. 

2. 

Scene error modeling for 

correspondence . 



uncorrected or partially 
processed MSS data. 

2. Printer listing. 


3. 

Output space definition. 



4. 

Interpolation grid point 
creation from error models . 



5. 

Interpolation grid point 
creation from ancillary data 
(for partially processed 
scenes) . 



6. 

Geometric error assessment 
from error models. 




Table A-7. Resampling Program 


Input 

1. Ancillary data set from 
geometric transformation 
program. 

2. MIST format tape from 
uncorrected or partially 
processed scene. 

3. User’s requests. 

Checkpoint restart data set. 


Process 


1 . Set up interpolation grid 
point arrays 

2. Horizontal, high-frequency 
error correction 

3. Cubic convolution resampling. 

4. Nearest neighbor resampling. 
Checkpoint restart capability. 


Output 

1. MIST format tape. 

2. Printer listing. 

3. Checkpoint restart data set. 


4. 


5. 
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B. 1 INTRODUCTION 


This document describes in detail the software that comprises the SAR/LANDSAT 
(SLDMS) . Its primary purpose is to specify the software in sufficient detail 
so that a competent, image-processing programmer can write the SLDMS programs 
It is assumed that the programmer will be somewhat familiar with the IBM 
existing image-processing software and that he (or she) will use or modify 
existing code wherever possible. 

A second purpose of this document is to provide NASA with enough knowledge to 
be able to understand the characteristics of the SLDMS. 

The design of each SLDMS program is described by a module hierarchy diagram 
and a set of program specifications (one for each module). In general, the 
program specifications are independent of the particular programming language 
chosen for implementation. 
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B.2 SYSTEM DATA SETS 


The SLDMS software uses various input and output data sets. These data sets are 
described in this section at a system level. A cross-reference of the data sets 
and programs used in the SLDMS is contained in Table B-l. Each of these data 
sets is described below in some detail. While it is expected that the formats 
of some data sets will be changed during the SLDMS implementation effort, the 
descriptions in this document should give a reasonable indication of the amounts 
and kinds of data that will be used by the SLDMS . 


B.2.1 X-Eormat CCT 

This tape, purchased from EROS Data Center (EDC) or obtained from NASA/GSFC, con- 
tains uncorrected MSS image data. There are five files on this tape. The format 
is completely described in tbe following: 

"Generation and Physical Characteristics of the Landsat 1 and 2 MSS 
Computer Compatible Tapes," Goddard Space Flight Center, November 1975, 

NASA Publication X-563-75-223 . 


B. 2.2 BSQ-Format CCT 

Partially processed MSS data and fully corrected MSS data can be produced on CCTs 
by the MDP, The formats of these tapes are defined m the following Interface 
Control Documents between the Image Processing Facility and EDC Digital Image 
Processing System for Landsat: 

"Fully Processed Multispectral Scanner Computer Compatible Tape (CCT-PM)," 
IBM/FSD, Contract NAS5-22999 

"Partially Processed Multispectral Scanner Computer Compatible Tape 
(CCT-AM) IBM/FSD, Contract NAS5-22999 


B.2. 3 BIL-Format CCT 

No documentation on BIL-format CCTs is currently available. 

B.2. 4 MIST-Format CCT 

The MIST format will be used for all image files created by the SLDMS. It is a 
BIL format in which one record contains data from all spectral hands of a given 
line. The format is completely described in the following: 

"LARSYS System Manual," Pages 5-52 through 5-58. 


i. 
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Table B-l. Data Set/Program Cross Reference 


Data Set 




Program 

MRFT 

ACPL 

MCPL 

GTRN 

RSPL 

X-Format CCT 

I 





BSQ-Format CCT 

I 





BIL-Format CCT 

I 





MiST-Format CCT 

0 

I 

10 


10 

CPL Copy Tape 


I 




Checkpoint Tape 





10 

Ancillary Data Set 

0 

10 


10 

I 

Shadeprint Data Set 



0 



User-Request Data Set 

I 

I 

I 

I 

I 

Temporary Data Set 

X 
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B.2.5 CPL Copy Tape 


The Control-Point Library (CPL) Copy Tape provides a copy of either the opera- 
tional sub-library or the delta sub-library of a library residing on a given MDP 
system. Its use in NASA's Control-Point Library Building System (CPLBS) is to 
transfer the operational library to the other system, to restore the operational 
library, to update operational libraries on both systems, or to restore the delta 
library. 

The CPL Copy Tape will provide prototype control points (windows) to users of 
the SLUMS. For those Landsat MSS scenes that have corresponding control points 
on the tape, an SLDMS user will be able to perform precise geometric correction 
(to his particular scale and orientation) with an automatic process. 

The format of the CPL Copy Tape is currently under review at IBM. The format 
described below is not necessarily the actual format that will he used. 

The tape contains in the following order: 

a. Library Control Record 

b. All directory records 

c. All control-point records 

The directories and control points are copies of the corresponding CPLBS disk 
files. In general, the data will require more thaa one tape reel. Each reel 
will contain the Library Control Record followed by the remaining data. The 
Library Control Record will insure that the reel is part of the sequence. 

Detailed descriptions of the above records follow. 

The directory file and the control-point files on the CPL Copy Tape contain 
data blocked at five records per block. Record lengths and block lengths are 
shown in Table B-2. In the case of the control-point file, block size is not a 
multiple of record length. Record formats are shown in Tables B-3, B-4, and B-5. 


Table B-2. Control Point Library Copy Tape Record Formats 


Record Type 

Record Length 

Block Length 

Maximum # Records 

Library Control Record 

400 

400 

1 per reel 

Directory Record 

360 

1800 

6250 

Control-Point Record 

1244 

6232 

125000 
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TABLE B-3. 


PARAMETER 

DMCFFNAM 

DMCFPROG 

DMCFDATE 

DMCFTIME 

DMCFSTAT 

DHCFVALD 

DMCFSYSI 

DMCFEEC# 

DHCFGPNM 

DMCFLBNM 

DMCFSLNM 

BMCFLBBT 

DMCFINSV 

DMCFTNNM 

DMCFDRSV 

DMCFDRNM 

DMCFCPSV 

DMCFCENM 


BMCFDRA# 

EMCFDR#D 

DMCFCPA# 

DMCFCP#D 

DMCFSLUP 

DHCFCTS# 

DMCFRINP 

DMCFRT# 

DMCFRT#! 

DMCFBT#A 

DMCFE T#B 


DMCFPMUR 

DMCFPBSV 

BMCFUINP 

DMCFUT# 

DMCFUT#I 

DMCFUT# A 

DMCFUT#B 

DMCFES# 

DMCFTAPC 


LIERARY CCNTBOL RECORD 


ZfSS:2& 


DESCRIPTION 

ALPHANUMERIC SYSTEM NAME OF FILE 
NAME OF PROGRAM THAT WROTE THIS RECORD 
DATE RECORD WAS LAST WRITTEN 
TIME RECORD WAS LAST WRITTEN 
RECORD STATUS 

SUE-LIBRAPY VALID INDICATOR 

SYSTEM ON WHICH THIS SUB-LIERARY RESIDES 

CONTROL RECORD NUMBER IN DMTL 

GROUP NAME ( PRODUCTION OR DEVELOPMENT ) 

LIBRARY NAME ( MSS OR RBV ) 

SUB-LIBRARY NAME ( OPERATIONAL OR DELTA ) 

DATE LIERARY INITIALIZED 
SPAPE 

INDEX FILE SNT VALUE 

INDEX FILE NAME 

DIRECTORY FILE SNT VALUE 

DIRECTORY FILE NAME 

CONTPOL POINT FILE SNT VALUE 

CCNTECL POINT FILE NAME 

ALIGNMENT 

SPARE 

NEXT AVAILABLE DIRECTORY RECORD NUMBER 
NUMBER OF DELETED DIRECTORY RECORDS 
NEXT AVAILABLE CONTROL POINT RECORD NUMBER 
NUMBER OF DELETED CONTROL POINT RECORDS 
SPARE 

DATA WRITTEN TO SUB-LIERARY SINCE LAST COPY TAPE? 
TAPE NUMBER OF LAST COPY TAPE PRODUCED 
RE STORE-IN -PROG BESS ( YES OR NO ) 

SEQUENTIAL TAPE # OF TAPE USED IN LAST RESTORE 
SYSTEM ID ( A OR B } 

SEQ. TAPE # OF LAST RESTORE TAPE FROM SYSTEM A 

SEQ. TAPE # OF LAST RESTORE TAPE FROM SYSTEM B 

ALIGNMENT 
SPARE 

NEXT PAGE NUMBER FOR PMUR REPORT 
PERFORMANCE BATA FILE SNT 
UPDATE-IN-PROGRESS ( YES OR NO ) 

SEQUENTIAL TAPE # OF TAPE USED IN LAST UPDATE 
SYSTEM ID ( A OR B 1 

SEQ. TAPE # OF LAST UPDATE TAPE FROM SYSTEM A 

SEQ. TAPE # OF LAST UPDATE TAPE FROM SYSTEM B 

SPARE 

REEL NUMBER 

CONTENTS OF TAPE (DELTA, OPERATIONAL, PERFORMANCE) 

ALIGNMENT 

SPARE 


COLUMNS FORMAT 


1 

- 

8 

C8 

9 

- 

12 

C4 

13 

— 

16 

C4 

17 

— 

22 

C6 

23 

— 

24 

X2 

25 

— 

26 

Y2 

27 

- 

28 

12 

29 

— 

30 

12 

31 

— 

32 

C2 

33 

— 

34 

C2 

35 

— 

36 

C2 

37 

— 

44 

C8 

45 

— 

48 

A4 

49 

- 

50 

12 

51 

— 

58 

C8 

59 

— 

60 

12 

61 

— 

68 

C8 

69 

— 

70 

12 

71 

- 

78 

C8 

79 

- 

80 

A2 

81 

— 

84 

A4 

85 

— 

88 

14 

89 

— 

92 

14 

93 

— 

96 

14 

97 

— 

10 0 

14 

10 1 

- 

10 4 

A 4' 

105 

— 

106 

Y2 

107 

— 

10 8 

12 

109 

— 

110 

Y2 

111 

— 

112 

12 

113 

— 

114 

12 

115 

- 

116 

12 

117 

- 

11 8 

12 

119 

— 

120 

A2 

121 

- 

124 

A4 

125 

- 

128 

14 

129 

— 

130 

12 

131 

— 

13 2 

Y2 

133 

- 

134 

12 

135 

- 

136 

12 

137 

— 

138 

12 

139 

— 

140 

12 

141 

— 

144 

A4 

145 

— 

146 

12 

147 

— 

15 8 

Cl 2' 

159 

— 

160 

A2 

161 

- 

40 0 

A 2 40 
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TABLE B-4. CONTROL POINT LIBRARY DIRECTORY RECORD 


PARAMETER 

DM BREN AH 
DMDPPEOG 
DMDR DATE 
DMD5TIHE 
D MDRSTAT 
DM DR CHAN 

DMDRSNSR 

D MDRFATH 
DMDFPOWN 
DMDRQUAD 
DMDRREGT 
DMDRMISN 
DMDROBSV 
DMDRUPDL 
DMDRMODL 
DMDRCPCT 

DMDRCIDS 

DMDRCRBS 

D BDRLRUl 

DMDRLBU2 

DMDRLRD3 

DMDRLRU4 

DM DR SR U1 

DMDRSRD2 

DMDRSEU3 

DMDPSRD4 

DMDRL PHI 

DMDRLRH2 

DMDRLPH3 

DMDRLHH4 

DMDRSRH1 

DM DP SR H2 

DKDPSEH3 

DMDRSEH4 

DMDRLWU 

DMDRSWU 

DMDRLWB 

DMDRSWB 

DMDPNGLA 

DMDRNLON 

DMDEBETA 

DMDCCNCT 

DMDCDDAT 
DM DCUUID 
DMDCUSYS 
DMECUTAP 
DMDCMDAT 
D MDCMUID 
DMDCMSYS 
DMDCMTAP 


DE SCRIPTION 

ALPHANUMERIC SYSTEM NAME OF FILE 
NAME OF PEOGEAM THAT WROTE THIS RECORD 
DATE RECORD WAS LAST WRITTEN ( YDDD ) 

TIME RECORD WAS LAST WRITTEN ( HHBMSS ) 

RECORD STATUS 

CHAIN POINTER - RELATIVE RECORD NUMBER 

SPARE 

SENSOR 

SPARE 

WES PATH NUMBER 
WES BOW NUMBER 
EBV IMAGE QUADRANT 
REGISTRATION TYPE 
LANDSAT MISSION 

IMAGE OBSERVATION TIME: DDDDHHMMS 

UPDATE LEVEL OF CONTROL POINT SET 

MODIFICATION LEVEL OF CONTROL POINT SET 

# OF ENTRIES IN CONTROL POINT ID AND POINTER LISTS 

SPARE 

CONTROL POINT ID NUMBER LIST 
CONTROL POINT POINTER LIST 

INPUT-SPACE LINE (UTM/PS) - CORNER 1 (EBV ONLY) 

INPUT-SPACE LINE (UTM/PS) - CORNER 2 (RBV ONLY' 

INPUT-SPACE LINE (UTM/PS) - CORNER 3 (RBV ONLY' 

INPUT-SPACE LINE jUTM/PS[ - CORNER 4 (RBV ONLY) 

INPUT-SPACE SAMPLE (UTM/PS) - CORNEE 

INPUT-SPACE SAMPLE (UTM/PS) - CORNEE 2 

INPUT-SPACE SAMPLE (UTM/PS) - CORNER 3 

INPUT-SPACE SAMPLE (UTM/PS) - CORNEE 4 

INPUT-SPACE LINE (HOM) - CORNER 1 (RBV ONLY) 

INPUT-SPACE LINE (HOM) - CORNER 2 (RBV ONLY) 

INPUT-SPACE LINE (HOM) - CORNER 3 (EBV ONLY 

INPUT-SPACE LINE jHOM) - CORNER 4 (RBV ONLY) 

INPUT -SPACE SAMPLE (HOM) - CORNER 
INPUT-SPACE SAMPLE (HOM) - CORNER 2 
INPUT-SPACE SAMPLE (HOM) - CORNER 3 
INPUT-SPACE SAMPLE (HOM) - CORNER 4 
INPUT-SPACE LINE POSITION (UTM/PS) - WRS CENTER 
INPUT -SPACE SAMPLE POSITION - WRS CENTER 
INPUT-SPACE LINE POSITION (HOM) - WRS CENTER 
INPUT-SPACE SAMPLE POSITION (HOM) - WRS CENTER 
GEODETIC LATITUDE OF NADIR 
LONGITUDE, OF NADIR 
BETA ANGLE OF REFERENCE IMAGE 
CONTROL POINT ID NUMBER COUNTER 
SPARE 

DATE OF UPDATE FOR LATEST UPDATE LEVEL 

ID OF USER WHO ENTERED UPDATE 

ID OF SYSTEM WHERE UPDATE PERFORMED 

SEQUENTIAL TAPE NUMBER OF UPDATE TAPE 

DATE OF UPDATE FOR LATEST MODIFICATION LEVEL 

ID OF USER WHO ENTERED MODIFICATION 

ID OF SYSTEM WHERE MODIFICATION PERFORMED 

SEQUENTIAL TAPE NUMBER OF UPDATE TAFE 

SPARE 


COLUMNS FORMAT 


1 

— 

8 

C8 

Q 

— 

12 

C4 

13 

— 

16 

C4 

17 

— 

22 

C6 

23 

— 

24 

X2 

25 

— 

28 

14 

29 

— 

32 

A4 

33 

— 

34 

12 

35 

— 

36 

A2 

37 

— 

38 

12 

39 

— 

40 

12 

41 

— 

42 

12 

43 

— 

44 

12 

45 

— 

46 

C2 

47 

— 

56 

CIO 

57 

— 

58 

12 

59 

— 

60 

12 

61 

— 

62 

12 

63 

— 

64 

A2 

65 

— 

124 

3012 

125 

- 

24 4 

3014 

245 

— 

24 6 

12 

247 

— 

248 

12 

249 

— 

250 

12 

251 

- 

252 

12 

253 

— 

254 

12 

255 

— 

256 

12 

257 

- 

258 

12 

259 

• 

260 

12 

261 

— 

262 

12 

263 

— 

264 

12 

265 

- 

266 

12 

267 

- 

268 

12 

269 

• 

270 

12 

271 

— 

27 2 

12 

273 

— 

27 4 

12 

275 

— 

276 

12 

277 

— 

27 8 

12 

279 

— 

280 

12 

281 

- 

282 

12 

283 

— 

28 4 

12 

285 

- 

292 

D8 

293 

— 

300 

D8 

301 

— 

30 8 

D8 

309 

— 

310 

12 

311 

- 

312 

A2 

313 

— 

316 

14 

317 

- 

324 

C8 

325 

— 

32 6 

12 

327 

- 

328 

12 

329 

- 

33 2 

14 

331 

- 

34 0 

C8 

341 

- 

342 

12 

343 

— 

344 

TO 

345 

- 

360 

A16 
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TABLE B-5. 


PARAMETER 

DMCPFNAM 
DMCPPROG 
DMCPDATE 
DMCP TIME 
DMCPSTAT 
DMCPSNSR 

DMCPPATH 
DMCPRCWN 
DMCPQUAD 
DMCPREGT 
DMCPMISN 
DMCPOBSV 
DHCPDPDL 
DMCPMODL 
DMCPCNCT 
DMCPWTYP 
DMCPLOLI 
D MCPLOSM 
DMCPPXMN 
DMCPSHFT 
DMCPPXVR 
DMCPMMAG 
DMCPLSBR 
DMCPLSBC 

DMCPWNDW 

DMCPCOFS 

DMCPSPBD 

DMCPLXEC 

DMCPLYEC 

DMCPLZEC 

DMCPEEOT 

BMC SCOBS 
DMCSGIEC 
DMCSMPSC 
DMCSNCDE 
DMCS7ACR 
DMCSOOTS 
DMCSEDGE 
DMCSHIMG 
DMCSNCOR 
DMCSPEAK 
EMCSCURV 
DMCSSPEK 
DMCSSCRV 
BHCSLOLT 
DMCSLCLN 
DMCSELEV 
DMCSEBAT 
EMCSERCT 
DMCSUSER 
DMCSCDAT 
EMCSMCIB 
BMCSCSYS 
DMCSCTAP 


LIERARY CONTROL POINT RECORD 


okigmal Pap „ 

i 00 - R quality 


DESCRIPTION 

ALPHANUMERIC SYSTEM NAME OF FILE 
NAME OF PFOGRAM THAT WROTE THIS RECORD 
DATE RECORD WAS LAST WRITTEN ( YDDD ) 

TIME RECORD WAS LAST WRITTEN ( HHMMSS ) 

RECORD STAIRS 

SENSOR 

SPARE 

WPS PATH NUMBER 
WRS PCW NUMEER 
RBV IMAGE QUADRANT 
REGISTRATION TYPE 
LANDS AT MISSION 

IMAGE OBSERVATION TIME: DDDDHHMMS 
UPDATE LEVEL OF CONTROL POINT SET 
MODIFICATION LEVEL OF CONTROL POINT SET 
CONTROL POINT ID NUMBER WITHIN SET 
WINDOW TYPE ( H, C, F. OS P ) 

FEATURE LOCATION IN WINDOW IN LINES 

FEATURE LOCATION IN WINDOW IN SAMPLES 

MEAN OF WINDOW PIXEL VALUES {DECIMAL AFTER BIT 8) 

SHIFT FACTOR ( PO ) 

VARIANCE OF WINDOW PIXEL VALUES X WINDOW SIZE 
MEAN MAGNITUDE OF CONDITIONED WINDOW (PI) 

# LEADING SIGN BITS IN LARGEST ROW FFT COMPONENT (P2 

# LEADING SIGN BITS IN LARGEST COL FFT COMPONENT (P3 

SPARE v 

WINDOW ARRAY ( 32-BY-32 ARRAY OF 8-BIT PIXELS ) 

ORDER OF SELECTION BY EDIT: 1 - N 

SPECTRAL BAND NUMBER 

GEOCENTRIC RECTANGULAR X COORD OF CONTROL POINT 
GEOCENTRIC RECTANGULAR Y COORD OF CONTROL POINT 
GEOCENTRIC RECTANGULAR Z COORD OF CONTROL POINT 
VELOCITY OF CONTROL POINT DUE TO EARTH ROTATION 
SPARE 

SOURCE IMAGE ID: A. EDDD.HHMMS 
GENERIC IDENTIFICATION CODE 
MAP SCALE CODE 
NODE; DESCENDING/ASCENDING 

# IMAGES IN WHICH CP VIOLATED 7-ACRCSS RESTRICTION 

# IMAGES IN WHICH EST. CP LCC. OUTSIDE AVAIL. DATA 

# IMAGES IN WHICH LOC. OF PEAK TOO NEAR EDGE OF SA 

# IMAGES FOR WHICH CORRELATION WAS ATTEMPTED 
NUMBER OF . SUCCESSFUL CORRELATIONS 

MEAN OF PEAK CORRELATION VALUES OBSERVED 
MEAN, OBS. VALUES OF MIN. CURVATURE OF COR SURF. 

PEAK CORRELATION VALUE 

CURVATURE OF CORRELATION SURFACE AT PEAK 
FEATURE LOCATION IN GEODETIC LATITUDE 
FEATURE LOCATION IN LONGITUDE 

ELEVATION OF FEA.TUPE {METERS ABOVE REF ELLIPSOID) 
AXIS 90 % ERROR ELLIPSE EST. ALONG TRACK (METERS) 

AXIS 90% ERROR ELLIPSE EST. ACROSS TRACK (METERS) 

ID OF USER WHO CREATED THIS POINT 
CREATION DATE - YDDD {JULIAN DAY) 

NASA CP ID # FROM MAP CP DESCRIPTION DATA TAPE 

ID OF SYSTEM ON WHICH CREATED 

ID OF TAPE WHICH ENTERED INTO OPER'L LIB. 

SPARE 


COLUMNS FORMAT 


1 

— 

8 

C8 

9 

- 

12 

C4 

13 

— 

16 

C4 

17 

- 

22 

C6 

23 

— 

24 

X2 

25 

— 

26 

12 

27 

— 

28 

A2 

29 

— 

30 

12 

31 

— 

32 

12 

33 

— 

34 

12 

35 

— 

36 

12 

37 

— 

38 

C2 

39 

— 

48 

CIO 

49 

— 

50 

12 

51 

— 

52 

12 

53 

— 

54 

12 

55 

- 

56 

C2 

57 


60 

E4 

61 

- 

64 

E4 

65 

— 

66 

12 

67 

— 

68 

12 

69 

— 

72 

E4 

73 

- 

74 

12 

} 75 

— 

76 

12 

) 77 

— 

78 

12 

79 

— 

80 

A2 

81 

— 

1104 


1105 

— 

1106 

12 

1 107 


1108 

12 

1109 

— 

11 16 

D8 

1117 

- 

1124 

D8 

1125 

— 

1132 

D8 

1133 

— 

11 4C 

D8 

1141 

— 

1144 

A 4 

1 145 

— 

1156 

Cl 2 

1157 

— 

1160 

C4 

1161 

— 

1164 

14 

1165 

— 

1166 

12 

1167 

— 

1168 

12 

1169 

— 

1170 

12 

1171 

- 

1172 

12 

1173 

— 

1176 

14 

1177 

— 

1180 

14 

1181 

— 

11 84 

E4 

1185 

— 

1188 

E4 

1189 

— 

1192 

E4 

1193 

— 

11 96 

E4 

1197 

— 

1204 

D8 

1205 

— 

1212 

D8 

1213 

— 

1216 

E4 

1217 

— 

12 20 

E4 

1221 

- 

1224 

E4 

1225 

— 

1232 

C8 

1233 

— 

12 36 

14 

1237 

— 

1238 

C2 

1239 

- 

1240 

12 

1241 

— 

1242 

12 

1243 

- 

1244 

A2 
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B.2.6 Checkpoint Tape 


The Checkpoint Tape is created and used by the Resampling Program (RSPL) . Its 
purpose is to save data whenever a checkpoint is taken. Whenever a computer 
failure occurs, this data set can Be used to restart RSPI at the last checkpoint. 
This reduces the amount of lost computer time in the event of a failure during 
the expensive resampling program. 

The format of the Checkpoint Tape will be defined during the implementation of 
the SLDMS. 


B.2.7 Ancillary Data Set 

Every MSS scene that is processed by the SLDMS will have an associated Ancillary 
Data Set (ADS) . It contains information that completely describes the associated 
scene as well as information produced by one SLDMS program for use by another 
program. 

The ADS will be a sequential data set consisting of 80-byte records. It will 
usually reside on a direct-access device, although it may be a tape data set. 

It is organized into sections, each consisting of a title record followed by one 
or more data records. All data in a single section is functionally related. 
Figure B-l shows the layout of the ADS. 



Figure B-l. Organization of Ancillary Data Set 


Most of the data on the ADS is expected to be in printable character form. This 
will permit SLDMS users to easily determine what data is on an ADS by viewing it 
at a terminal or printing it with a standard utility. The basic formats of the 
two types of ADS records are shown in Tables B-6 and B-7. 

The formats of the data cards m each "section" will be determined during the 
implementation of the S1DMS. Preliminary definitions of all ADS "sections" 
follow in Tables B-8 through B-20. 
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TABLE B-6. ANCILLARY DATA SET; TITLE CABD FORMAT 


PARAMETER 

DESCRIPTION 

COLUMNS 

FORMAT 

TITLE 

A UNIQUE DESCRIPTIVE IDENTIFIER THAT DEFINES A 

"SECTION" OF THE ADS (IN PRINTABLE CHARACTER 
FORM). 

1 - 60 

FORM) 

C60 

NR EC 

NUMBER OF DATA RECORDS IN THIS "SECTION". 

61 - **0 

no 

NCARD 

SEQUENCE NUMBER OF THIS CARD WITHIN THE ADS (NOT 
REQUIRED TO BE USED) . 

73 - 80 

18 


TABLE B-7. ANCILLABY EATA SET: EATA CABD FORMAT 


PARAMETER DESCRIPTION 

ELANK CHARACTER ( X * 40 • ) . 

DATA FIELDS (IN PRINTABLE CHARACTER FORM IF 
POSSIBLE) . 

NCASE SEQUENCE NUMBER OF THIS CARD WITHIN THE ADS (NOT 

REQUIRED TO BE USED) . 





ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE B-8. FORMAT OF ANCILLARY DATASET 


X-FORMAT CCT HEADER DATA; 

SCENE ACQDATE 

NADLAT NADLON FCLAT 


CCT 


SI AT DATA: 
NATEST 
ATTTIM1 
ATTTIM2 
ATTTIB3 
ATTTIM4 
ATTTIM5 
ATTTIM6 
ATTTIM7 
ATTTIM8 
ATTTIM9 


NALEST 

AMSROL1 

AMSEOL2 

AMSROL3 

AMSRCL4 

AMSROL5 

AMSROL6 

AMSR0L7 

AMSROL8 

AMSROL9 


ALFER 
AMSPIT1 
AMSPIT2 
AMSPIT3 
AMSPIT4 
AM SPITS 
AMSPIT6 
AMSPIT7 
AMSPIT8 
AMSPIT9 


ACQSITE DIRECT/RECOFDED 

FCLON 


DELVV 
AMSYAW1 
AM SYAW2 
AMSYAE3 
AMSYAR4 
AMSYAR5 
AMSYAN6 
AMSYAR7 
AMSYAR8 
AMSYAR9 


AMSNLT 
AITTIM1 
ALTTIM 2 
ALTTIB3 
ALTTIH4 
AITTIS5 
ALTTIM 6 
ALTTIM7 
ALTTIM8 
AITTIM9 
ALTTIH10 
ALTTIM11 


AMSNLN 
AMS ALT 1 
AMSALT2 
AMS ALT 3 
AMSALT4 
AMS ALTS 
AMSALT6 
AMSALT7 
AMSALT8 
AMS ALT 9 
AMSALT10 
AMS ALT11 


3 

13 


DETECTOR -BY-DETECTOR HISTOGRAM: 


B1D1V0 

B1D2V0 

B1D3V0 

B ID 1 VI 
• 

B1D2V1 

• 

B1D3V1 

• 

• 

B1D1V2S5 

• 

B1D2V255 

• 

B1D3V255 

B2D1V0 

B2D2V0 

. B2D3V0 

B2D1V1 

• 

B2D2V1 

• 

B2D3V1 

« 

* 

B2D1V255 

B2D2V255 

B2D3V255 

B3D1 VO 

B3D2V0 

B3D3V0 

B3D1V1 

« 

B3D2V1 

• 

B3D3V1 

• 

• 

B3D1V255 

• 

B3D2V2S5 

• 

B3D3V2S5 

B4D1V0 

B4D2V0 

B4D3V0 

B4D 1 V 1 
• 

B4D2V1 

« 

B4D3V1 

• 

• 

B4D1V255 

• 

B4D2V2S5 

• 

B4D3V255 

B5D1V0 

B5D2V0 


BSD1V1 

• 

B5D2V1 

• 


« 

B5D1V255 

• 

BSD2V2S5 


DETECTOR-BY- 

DETECTOR STATISTICAL 

B1D1M 

E1D2M 

B1D3M 

B1D1S 

B1D2S 

B1D3S 

B2D1 M 

B2D2M 

B2D3M 

B2D1S 

B2D2S 

B2D3S 

B3D1M 

B3D2M 

B3D3M 

B3D1 S 

B3D2S 

B3D3S 

B4D1M 

B4D2M 

B4D3M 

B4D1S 

B4D2S 

B4D3S 

B5D1 M 

E5D2M 

B5D3M 

B5D1 S 

B5D2S 

B5D3S 

OUTPUT -SPACE 

DEFINITION: 

HSPACO 

VSPACO 

HEADING 

FDH 

FCV 

FCLAT 


1281 


B1D4V0 

B1D5V0 

B1D6V0 

B1 VO 

B1D4V1 

• 

B1D5V1 

« 

B1D6V1 

• 

B1V1 

• 

• 

B1D4V255 

• 

B1D5V2S5 

♦ 

B1D6V255 

• 

B IV 255 

B2D4V0 

B2D5V0 

B2D6V0 

B2V0 

B2D4V1 

• 

B2D5V1 

• 

B2D6V1 

• 

B2V1 

• 

• 

B2D4V255 

• 

B2D5V2S5 

• 

B2D6V255 

• 

B2V255 

B3D4V0 

B3D5V0 

B3D6V0 

B3V0 

B3D4Y1 

• 

B3D5V1 

« 

B3D6V1 

• 

B3V 1 

■ 

« 

B3D4V255 

• 

B3D5V255 

• 

B3D6V255 

* 

B3V255 

B4D4V0 

B4D5VO 

B4D6V0 

B4V0 

B4D4V1 

• 

B4D5V1 

• 

B4D6V1 

• 

B4V1 

• 

• 

B4D4 V255 

• 

B4D5V255 

• 

B4D6V255 

• 

B 4V 2 5 5 




B5V0 
B5V 1 

• 




• 

B5V255 


PARAMETERS: 


B1D6M 

11 

B1D4M 

B1D5M 

B1 M 

B1D4S 

B1 DBS 

B1D6S 

BIS 

B2D4M 

B2D5M 

B2D6M 

B2M 

B2D4S 

B2D5S 

B2D6S 

B2S 

B3D4M 

B3D5M 

B3D6M 

B3H 

B3D4S 

B3D5S 

B3D6S 

B3S 

B4D4M 

B4D5M 

B4D6M 

B4M 

B4D4S 

B4D5S 

B4D6S 

B4S 

B5D4M 

B5D5M 

B5D6M 

B5M 

B5D4S 

E5D5S 

BSD6S 

B5S 

OSWIDP 

OSLENP 

MAPPROJ 

3 

FCLON 

FCMAPl 

FCMAP2 

FCMAP3 


B-10 



TABLE ‘B-8. FORMAT CF ANCILLARY DATASET (CONTINUED) 


ERROR MODELS 


ROLLO 

ROLL1 

ROLL2 

ROIL3 



PITCHO 

PITCH1 

PITCH2 

PITCH3 



YAWO 

YAW 1 

YAW2 

YAW3 



DELHO 

DELHI 





Cl? ASSESSMENT: 





OSH1 

OSV1 

CTPH1 

CTPV1 

AESTH1 

AESTV1 

0 SH2 

• 

0 SV2 

• 

CTPH2 

• 

CTPV2 

• 

AESTH2 

« 

AESTV2 

• 

• 

0 SHN 

• 

0 SVN 

• 

CTPHN 

• 

CTPVN 

• 

AESTHN 

• 

AESTVN 

HEADER DATA: 





SCENE 

ASC/DSC 

WES PTH 

WRSROW 



SEN SOR 

MISSION 

ORBIT 

ADSTAT 

NUMDET 

ISWID 

WRSLIN 

WRSPIX 

CPY EAR 

CPDAY 

CPHOUR 

CPMIN 

MEANN 

MFANC 

GCFLG 

G CDF LG 

HCFLG 

FCDFLG 

ISLIN 

IFCRMAT 

ILTYPE 

ILCNT 

RSTYPE 

MAPROJ 

WRSOFF 

IDJUST 

UBANDS 

MBNDNO 



FASSIYPE 

CAITYPE 





' ANCILLARY 

RECORD NUMBER 1 DATA 




ISWID 

ISLEN 

ISPACH 

ISPACV 



OSWID 

OSLEN 

OSPACH 

OSPACV 

SCALTM 

ISWIDM 


MIRPO 


MIRR1 


MIRR2 


MIPR3 


MAXMIR 




SCSKTW 


TBTWN 


T ACTIVE 


A 


B 


ECCON 

SD1 1 

SD 1 2 

SD 1 3 

SD1 4 

SD15 

SD16 

SD21 

SD22 

SD23 

SD24 

SD25 

SD26 

SD31 

SD32 

SD33 

SD34 

SD35 

SD3 6 

SD41 

SD42 

SD43 

SD44 

SD45 

SD46 

SD51 

SD52 

BB0FF1 

BBOFF2 

BEOFF3 

BBOFF4 

1 ANCILLARY 

RECORD NUMBER 2 DATA 

- 



WESERM 

WES ORB 

WRSLAT 

WRSLON 

FCXCTIM 



NADLAT 


NADLON 




NADIRX 


NADIRY 


NADIR Z 


BETA 


INADLIN 


INADSPL 


EELVV 


ERVEL 


ERPARM 

ROLLO 

ROLL 1 

ROLL2 

ROL13 



PITCHO 

PITCH 1 

PITCH2 

PITCH3 



YAWO 

YAW1 

YAW2 

YAW3 



DELHO 

DELHI 





NUMGCP 






; POINTS FROM ANCILLARY RECORDS 

3 - 10: 



NHOR 

NVER 

{ K 

= NHOR B 1 V 1 NVER ) 


HI H 2 H3 

H4 H5 

H6 H7 H8 

H 9 Hi 0 

H1 1 HI 2 Hi 3 

HI 4 HI 5 

117 HI 8 HI 9 

H20 H21 

H22 H23 H24 

H2 5 H26 

H27 H28 H29 

H30 H31 

• • • 

• • 

• • • 

• « 

« a a 

♦ a 


HK 


( NOTE: STORED BOW BY ROW ) 


VPS 


POINTS FROM ANCILLARY RECORDS 3 - 10: 

NHOR NVER ( L = NHCR B1V1 NVER ) 

VI V2 V3 V4 V5 V6 V7 V8 V 9 V10 VII V12 V13 VI 4 
Vi 7 VI 8 V 19 V20 V21 V22 V23 V24 V25 V26 V27 V28 V29 V39 


VI 5 
V31 


VI 6 
V32 


VL 


( NOTE: STORED COLUMN BY COLUMN ) 


CCT ANCILLARY RECORD 10 DATA: 

■ WBSPIX WBSOFE ORIENT 


5 

N + 1 
8 

CPSEC 

12 


11 


K + 1 
L+1 
2 


B-ll 



ORIGINAL PAGE IS 
OF POOR QUALITY 

TABLE B-9. ANCILLARY EAT A SET FIELDS { X-FORMAT CCT HEADER DATA ) 


PARAMETER DESCRIPTION COLUMNS FORMAT 


TITLE 

X— FORMAT CCT HEADER DATA: 

1 

— 

60 

C60 

NREC 

NUMBER OF DATA RECORDS IN THIS "SECTION" ( * 3 ) 

61 

— 

70 

no 

— 

BLANK CHARACTER 



1 

Cl 

SCENE 

RASA SCENE IDENTIFIER 

2 

- 

11 

CIO 

ACODATE 

DATE OF ACQUISITION OF SCENE 

12 


21 

cio 

ACQSITE 

SITE OF ACQUISITION OF SCENE 

22 


31 

CIO 

DIRPEC 

DIRECT OR RECORDED DATA INDICATOR 

32 


41 

CIO 

— 

BLANK CHARACTER 



1 

ci 

NADLAT 

LATITUDE OE SPACECRAFT NADIP ( RADIANS ) 

2 

— 

11 

FI 0 

NADLON 

LONGITUDE OF SPACECRAFT NADIR ( RADIANS ) 

12 

— 

21 

F10 

FCLAT 

LATITUDE OT FORMAT CENTER ( RADIANS ) 

22 

— 

31 

FIG 

FCLON 

LONGITUDE OF FORMAT CENTER ( RADIANS ) 

32 

- 

41 

F10 


TABLE B-1C 

. ANCILLARY DATA SET FIELDS ( 

CCT SIAT DATA ) 





PARAMETER 

DESCRIPTION 



COLUMNS 

FORMAT 

TITLE 

CCT SIAT DATA: 



1 

— 

60 

C69 

NREC 

NUMBER OF DATA RECORDS IN THIS 

"SECTION" { = 13 ) 

61 

— * 

70 

110 

— 

BLANK CHARACTER 





1 

Cl 

NATEST 

LATITUDE OF SPACECRAFT NADIR 

( 

RADIANS ) 

2 

— 

11 

F 10 

NALSST 

LONGITUDE OF SPACECRAFT NADIR 

( 

RADIANS ) 

12 

- 

21 

P10 

ALFER 

EARTH ROTATION ANGLE f RADIANS 

) 

22 

— 

31 

n o 

DELVV 

NORMALIZED VELOCITY CHANGE 



32 

— 

41 

no 

A MSN IT 

LATITUDE OF SPACECRAFT NADIR 

( 

RADIANS ) 

42 

— 

51 

F 10 

AMSNLN 

LONGITUDE OF SPACECRAFT NADIR 

( 

RADIANS ) 

52 

— 

61 

P10 

— 

E1ANK CHARACTER 





1 

Cl 

ATTT IM 1 

TIME OF ATTITUDE MEASUREMENT 

1 


2 


11 

F10 

AMSFCll 

AMS RCII MEASUREMENT 1 



12 

— 

21 

F 10 

AM SPIT 1 

AMS PITCH MEASUREMENT 1 



22 

— 

31 

F 1 0 

AMSYAH1 

AMS YAH MEASUREMENT 1 



32 

- 

41 

E 1 0 

A1TTIM1 

TIME OP ALTITUDE MEASUREMENT 

1 


42 

- 

51 

•fI n 

AM SALT 1 

AMS ALTITUDE MEASUREMENT 1 



52 

— 

61 

FI 0 

• 

• 

• 

• 

• 

• 




• 

• 

• 


• 

« 

— 

BIANK CHARACTER 





1 

Cl 

ATTT IM 9 

TIME OF ATTITUDE MEASUREMENT 

9 


2 

- 

11 

FI C 

A MSEC! 9 

AMS ROLL MEASUREMENT 9 



12 

- 

21 

?1 rt 

AMSPIT9 

AMS PITCH MEASUREMENT 9 



22 

- 

31 

F10 

AMSYAH9 

AMS YAH MEASUREMENT 9 



32 

- 

41 

F10 

ALTTIM9 

TIME OE ALTITUDE MEASUREMENT 

9 


42 

— 

51 

FI 0 

AMSALT9 

AMS ALTITUDE MEASUREMENT 9 



52 

— 

61 

FI C 

— 

BLANK CHARACTER 





1 

Cl 

— 

SPARE 



2 

— 

51 

FI 0 

A1TTIM10 

TIME OF ALTITUDE MEASUREMENT 

10 


42 

— 

51 

F10 

AMSALT10 

AMS ALTITUDE MEASUREMENT 10 



52 


61 

F1C 

_ 

BLANK CHARACTER 





1 

Cl 

- 

SPARE 



2 

- 

51 

FI 0 

IAITTIM11 

TIME OF ALTITUDE MEASUREMENT 

11 


42 

- 

51 

FI "! 

AMSALT11 

AMS ALTITUDE MEASUREMENT 11 



52 


61 

F 1 0 


B"12 



TABLE B-11 


ANCILLARY DATA SET FIELDS ( DETECTOR-BY-DETECTOR HISTOGRAMS ) 


PARAMETER DESCRIPTION COLUMNS FORMAT 


TITLE 

NPEC 

DETECTOR-BY -DETECTOR HISTOGRAMS: 

NUMBER OF DATA RECORDS IN THIS "SECTION” 

( = 1281 ) 

1 

61 

- 

60 

70 

C60 

110 

— 

BLANK 

CHARACTER 







1 

Cl 

B1D1V0 

NUMBER 

OF 

BAND-1, DETECTOR-1 

PI XELS 

WITH 

VALUE 

0 

2 

- 

11 

110 

B1D2V0 

NUMBER 

OF 

BAND-1, DETECTOR-2 

PIXELS 

WITH 

VALUE 

0 

12 

— 

21 

no 

B1D3VC 

NUMBER 

OF 

BAND-1, DETECTOE-3 

PIXELS 

WITH 

VALUE 

c 

22 

- 

31 

no 

B1D4VC 

NUMBER 

OF 

BAND-1, DETECTOR-4 

PI XELS 

WITH 

VALUE 

c 

32 

— 

41 

no 

B1D5V0 

NUMBER 

OF 

BAND-1, DETECTOR-5 

PIXELS 

WITH 

VALUE 

0 

42 

• 

51 

no 

B1 D6V n 

NUMBER 

OF 

BAND-1, DETECTOR-6 

PIXELS 

WITH 

VALUE 

o 

52 

- 

61 

no 

B1 VO 

NUMBER 

OF 

BAND-1 PIXELS WITH 

VALUE 1 

0 



62 

- 

71 

no 

- 

BLANK « 

CHARACTER 







-? 

Cl 

B1D1 VI 

NUMBER 

OF 

BAND-1, DETECTOR-1 

PIXELS 

WITH 

VALUE 

1 

2 

- 

1i 

no 

B1D2V1 

NUMBER 

OF 

BAND-1, DETECTOR-2 

PIXELS 

WITH 

VALUE 

1 

12 

— 

21 

no 

B1D3V1 

NUMBER 

OF 

BAND-1, DETECTOR-3 

PIXELS 

WITH 

VALUE 

1 

22 

— 

31 

no 

B1D4V1 

NUMBER 

OF 

BAND-1, DETECTOR-4 

PIXELS 

WITH 

VALUE 

1 

32 

— 

41 

no 

B1D5V1 

NUMBER 

OF 

BAND-1, DETECTOR-5 

PIXELS 

WITH 

VALUE 

1 

42 

— 

51 

no 

B1D6V1 

NUMBER 

OF 

BAND-1, DETECTOR-6 

PIXELS 

WITH 

VALUE 

1 

52 

— 

61 

I m 

B 1 VI 

NUMBER 

OF 

BAND-1 PIXELS WITH 

VALUE 

1 



62 

- 

71 

no 

■ 

« 

• 


• 

« 

• 







• 

• 

• 


• 

• 

• 

- 

BLANK 

CHARACTER 







1 

Cl 

B1D1V255 

NUMBER 

OF 

BAND-1, DETECTOR-1 

PIXELS 

WITH 

VALUE 

255 

2 

— 

11 

no 

B 1D2V255 

NUMBER 

OF 

BAND-1, DETECTOR-2 

PI XELS 

WITH 

VALUE 

255 

12 

— 

21 

no 

B1D3V255 

NUMBER 

OF 

BAND-1, DETECTOR-3 

PIXELS 

WITH 

VALUE 

255 

22 

— 

31 

no 

B1D4V255 

NUMBER 

OF 

BAND-1, DETECTOR-4 

PIXELS 

WITH 

VALUE 

255 

32 

— 

41 

no 

B1D5V255 

NUMBER 

OF 

BAND-1, DETECTOR- 5 

.PIXELS 

WITH 

VALUE 

255 

42 

— 

51 

1 10 

B1D6V255 

NUMBER 

OF 

BAND-1, DETECTOR-6 

PIXELS 

WITH 

VALUE 

255 

52 

— 

61 

no 

B1V255 

NUMBER 

OF 

BAND-1 PIXELS WITH 

VALUE 255 



62 

- 

71 

no 

m 

* 

« 


* 

• 

• 







• 

m 

• 


• 

• 

• 

— 

BLANK CHARACTER 







1 

Cl 

B4D1V255 

NUMBER 

OF 

BAND-4, EETECTOR-1 

PIXELS 

WITH 

VALUE 

255 

2 


11 

T 10 

E4D2V255 

NUMBER 

OF 

BAND-4, DETECTOR-2 

PIXELS 

WITH 

VALUE 

255 

12 

- 

21 

no 

B4D3V2S5 

NUMBER 

OF- 

BAND- 4, DETECTOR-3 

PI XELS 

WITH 

VALUE 

255 

22 

- 

31 

no 

B4D4V255 

NUMBER 

OF 

BAND-4, DETECTOR-4 

PIXELS 

WITH 

VALUE 

255 

32 

— 

41 

i n 

B4D5V255 

NUMBER 

OF 

BAND-4, DETECTOR-5 

PIXELS 

WITH 

VALUE 

255 

42 

— 

51 

no 

B4D6V255 

NUMBER 

OF 

BAND-4, DETECTOR-6 

PI XELS 

WITH 

VALUE 

255 

52 

— 

61 

no 

B4V255 

NUMBER 

OF 

BAND-4 PIXELS WITH 

VALUE : 

255 



6-2 

- 

71 

no 

- 

BLANK < 

CHARACTER 







1 

Cl 

B5D1V0 

NUMBER 

OF 

BAND-5, DETECTOR-1 

PIXELS 

WITH 

VALUE 

0 

2 

— 

11 

no 

B5D2V0 

NUMEER 

OF 

BAND- 5, DETECTOR-2 

PIXELS 

WITH 

VALUE 

0 

12 

— 

21 

no 

“ 

SPARE 







22 

— 

61 

no 

e5vo 

NUMBER 

OF 

BAND- 5 PIXELS WITH 

VALUE 0 



62 

- 

71 

110 

• 

• 

• 


• 

• 







• 

• 


• 

* 

- 

BLANK CHARACTER 







1 

Cl 

B5D1 V255 

NUMBER 

OF 

BAND-5, DETECTOR-1 

PIXELS 

WITH 

VALUE 

255 

2 

— 

11 

no 

B5D2V2 55 

NUMBER 

OF 

BAND-5, DETECTOR-2 

PIXELS 

WITH 

VALUE 

255 

12 

— 

21 

no 


SPARE 







22 

— 

61 

no 

B5V255 

NUMBER 

OF 

BAND-5 PIXELS WITH 

value : 

255 



62 

— 

71 

no 


B-13 



origWal page is 

OF POOR QUALITY 


TABLE B-12. ANCILLARY DATA SET FIELDS { DETECTOR-BY-DETECTOR STATISTICAL DATA ) 


PARAMETER DESCRIPTION COLUMNS FORMAT 


TITLE 

DETECTOE-BY-EETECTOR STATISTICAL PARAMETERS: 

1 

— 

60 

C60 

NR EC 

NUMBER OF DATA RECORDS IN THIS "SECTION" 

( = 11 ) 

61 

— 

70 

no 

— 

BLANK CHARACTER 




1 

cl 

B1D1M 

MEAN OF BAND-1, DETECTOR-1 PIXELS 


2 

— 

11 

F 10 

B1D2M 

MEAN OF BAND-1, DETECTOR-2 PIXELS 


12 

— 

21 

F 10 

B1D3M 

MEAN OF BAND-1, DETECTOR-3 PIXELS 


2 2 

• 

31 

FI 0 

E1D4N 

MEAN OE BAND-1, DETECTOR-4 PIXELS 


32 

— 

41 

F 10 

B1D5M 

MEAN OF BAND-1, DETECTOR-5 PIXELS 


42 

— 

51 

F10 

B1D6M 

MEAN OF BAND-1, DETECTOR-6 PIXELS 


52 

- 

61 

F10 

B1M 

MEAN OF BAND-1 PIXELS 


62 

— 

71 

F 10 

— 

BLANK CHARACTER 




1 

Cl 

B1D1 S 

STANDARD DEVIATION OF BAND-1, DETECTOR-1 

PIXELS 

2 

— 

11 

FID 

B1D2S 

STANDARD DEVIATION OF BAND-1, DETECTOR-2 

PIXELS 

12 

- 

21 

F10 

B1D3S 

STANDARD DEVIATION OF BAND-1, DETECTOR-3 

PIXELS 

22 

— 

31 

F10 

B1D4S 

STANDARD DEVIATION OF BANE-1, DETECTOR-4 

PIXELS 

32 

— 

41 

FI 0 

B1D5S 

STANDARD DEVIATION OF BAND-1, DETECTOR-5 

PIXELS 

42 

— 

51 

no 

B1D6S 

STANDARD DEVIATION OF BAND-1, DETECTOR-6 

PIXELS 

52 

— 

61 

FID 

BIS 

STANDARD DEVIATION OF BAND-1 PIXELS 


62 

— 

71 

no 

• 

* 

• 

• 

• 

« 



• 

• 

• 


• 

« 

• 

— 

BLANK CHARACTER 




1 

Cl 

B5D1M 

BEAN OF BAND-5, DETECTOR-1 PIXELS 


2 

— 

11 

FI 0 

B5D2M 

BEAN OF BAND-5, DETECTOR-2 PIXELS 


12 

— 

21 

F 10 

B5D3M 

MEAN OF BANE-5, DETECTOR-3 PIXELS 


22 

— 

31 

F10 

B5D4M 

MEAN OF BAND-5, DETECTOR-4 PIXELS 


32 

- 

41 

FI 0 

B5D5M 

MEAN OF BANE-5, DETECTOR-5 PIXELS 


42 

— 

51 

F 10 

B5D6M 

MEAN OF BANE-5, DETECTOR-6 PIXELS 


52 

- 

61 

F 1 0 

B5M 

MEAN OF BAND-5 PIXELS 


62 

— 

71 

F10 

— 

BLANK CHARACTER 




1 

Cl 

B5D1 S 

STANDARD DEVIATION OF BAND-5, DETECTOR-1 

PIXELS 

2 

— 

11 

F 1 0 

B5D2S 

STANDARD DEVIATION OF BANE-5, DETECTOR-2 

PIXELS 

12 

- 

21 

FI 0 

B5D3S 

STANDARD DEVIATION OF BAND-5, DETECTOR-3 

PIXELS 

22 

— 

31 

FI 0 

B5D4S 

STANDARD DEVIATION OF BAND-5, DETECTOR-4 

PIXELS 

32 

— 

41 

F 10 

B5D5S 

STANDARD DEVIATION OF BAND-5, DETECTOR-5 

PIXELS 

42 

— 

51 

Fin 

B5D6S 

STANDARD DEVIATION OF BAND-5, DETECTOR-6 

PIXELS 

52 

— 

61 

F10 

B5S 

STANDARD DEVIATION OF BAND-5 PIXELS 


62 


71 

FI 0 
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TAELE B-13. ANCILLARY DATA SET FIELDS ( OUTPUT-SPACE DEFINITION ) 


PARAMETER 


DESCRIPTION 


TITLE OUTPUT-SPCE DEFINITION: 

NR EC NUMBER OF DATA RECORDS IN THIS "SECTION" ( 


BLANK CHARACTER 

HSPACO DISTANCE BETWEEN OUTPUT-SPACE PIXEL CENTERS 

VS PACO DISTANCE BETWEEN OUTPUT-SPACE PIXEL CENTERS 

HEADING SPACECRAFT HEADING ANGLE 

OSWIDP WIDTH OF OUTPUT SPACE (PIXELS) 

OSLENP LENGTH OF OUTPUT SPACE (PIXELS) 

MAPPROJ MAP PROJECTION 


FCH 

FCV 

FCLAT 

FCLON 

FCMAP1 

FCMAP2 

FCMAP3 


BIANK CHARACTER 

FORMAT CENTER COORDINATE (HCRIZONTAI) 
FORMAT CENTER COORDINATE (VERTICAL) 
FORMAT CENTER LATITUDE (RADIANS) 
FORMAT CENTER LONGITUDE (RADIANS) 
FORMAT CENTER MAP COORDINATE 
FORMAT CENTER MAP COORDINATE 
FORMAT CENTER MAP COORDINATE 



COLUMNS 

FORM A1 


1 

— 

60 

C60 

3 ) 

61 

— 

70 

no 




1 

Cl 

METERS) 

2 

— 

11 

F10 

METERS) 

12 

— 

21 

F 17 


22 

- 

31 

F 1 0 


3 2 

— 

41 

FIO 


42 

— 

51 

FIG 


52 

— 

71 

C20 




1 

Cl 


2 

— 

11 

F10 


12 

— 

21 

?10 


22 

— 

31 

FIG 


32 

- 

41 

FIO 


42 

— 

51 

FI" 


52 

- 

61 

FIG 


62 

- 

71 

F10 


TABLE B-14. ANCILLARY DATA SET FIELDS ( ERROR MODELS ) 


PARAMETER 


DESCRIPTION 


COLUMNS FORMAT 


TITLE 

EPROP MODELS: 





1 

— 

60 

C67 

NREC 

NUMBER OF DATA RECORDS IN 

THIS "SECTION" ( = 

5 ) 

61 


70 

im 

— 

BLANK CHARACTER 






1 

Cl 

rollo 

CONST ANT -TERM 

COEFEICIENT 

OF 

ROLL MODEL 


2 

— 

11 

v 1 0 

ROLL 1 

DEGREE- 1- TERM 

COEFFICIENT 

OF 

ROLL MODEL 


1 2 

— 

21 

F 10 

ROLL2 

DEGREE-2 -TERM 

COEFFICIENT 

OF 

ROLL MODEL 


22 

- 

31 

F 10 

ROLL 3 

DEGREE- 3- TERM 

COEFFICIENT 

OF 

ROIL MODEL 


32 

— 

41 

FI 0 


BLANK CHARACTER 






1 

Cl 

PITCHO 

CON ST ANT- TERM 

COEFFICIENT 

OF 

PITCH MODEL 


2 

— 

11 

FlO 

PITCH1 

DEGEEE-1-TERM 

COEFFICIENT 

OE 

PITCH MODEL 


12 

— 

21 

Tl 1 ' 

PITCH2 

DEGREE- 2- TERM 

COEFFICIENT 

OF 

PITCH MODEL 


22 

- 

31 

Fio 

PITCH 3 

DEGREE- 3- TERM 

COEFFICIENT 

OF 

PITCH MODEL 


32 

— 

41 

FlO 

— 

BLANK CHARACTER 






1 

Cl 

YAWO 

CONSTANT-TERM 

COEFFICIENT 

OF 

YAW MODEL 


2 

— 

11 

FIO 

YAW1 

DEGREE- 1 -TERM 

COEFFICIENT 

OF 

YAW MODEL 


12 

— 

21 

"Pi 0 

YAW2 

DEGREE- 2- TERM 

COEFFICIENT 

OF 

YAW MODEL 


22 

— 

31 

FIO 

YAW3 

DEG REE -3 -TERM 

COEFEICIENT 

OF 

YAW MODEL 


32 

— 

41 

FIO 

— 

BLANK CHARACTER 






1 

Cl 

DELHO 

CON STANT-TERM 

COEFFICIENT 

OF 

ALTITUDE-CHANGE 

MODEL 

2 

— 

11 

FIO 

DEIH1 

DEGREE-1-TERM 

COEFFICIENT 

OF 

ALTITUDE -CHANGE 

MODEL 

12 


21 

fio 
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TABLE B-15. ANCILLARY DATA SET HELDS ( ERROR ASSESSMENT ) 


PARAMETER DESCRIPTION COLUMNS FORMAT 


TITLE 

ERROR ASSESSMENT: 




1 

_ 

60 

C60 

NREC 

NUMEER OF DATA RECORDS IN THIS "SECTION 1 ' ( 

= N+1 ) 

61 

— 

70 

no 

— 

BLANK CHARACTER 






1 

Cl 

OSH1 

OUTPUT-SPACE COORDINATE OF POINT 1 




2 

— 

11 

FI 0 

OSV1 

OUTPUT-SPACE COORDINATE OF POINT 1 




12 

— 

21 

FI 0 

CTPH1 

CORRECTED TANGENT-PLANE COORDINATE 

OF 

POINT 

1 

22 

— 

31 

F 10 

CTPV1 

CORRECTED TANGENT -PLANE COORDINATE 

OF 

POINT 

1 

32 

- 

41 

F10 

AESTH1 

ERROR ESTIMATE FOR POINT 1 




42 

— 

51 

F10 

AESTV1 

ERROR ESTIMATE EOR POINT 1 




52 

— 

61 

F 10 

• 

• 

• 

* 

% 

• 





• 

• 

• 


• 

• 

m 

— 

BLANK CHARACTER 






1 

Cl 

CSHN 

OUTPUT-SPACE COORDINATE OE POINT N 




2 

— 

11 

F 1 0 

OS VN 

OUTPUT-SPACE COORDINATE OF POINT N 




12 

— 

21 

F 1 0 

CTPHN 

CORRECTED TANGENT-PLANE COORDINATE 

CF 

POINT 

N 

22 

— 

31 

F 1 0 

CTPYN 

CORRECTED TANGENT-PLANE COORDINATE 

OF 

POINT 

N 

32 

— 

41 

F10 

AESTHN 

ERROR ESTIMATE EOR POINT N 




42 

— 

51 

FI 0 

AEST7N 

ERROR ESTIMATE FOR POINT N 




52 

— 

61 

F 1 0 


ORIGINAL PAGE IS 
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TABLE 'B- 16. ANCILLARY DATA SET FIELDS { CCT HEADEH DATA ) 


PARAMETER DESCRIPTION COLUMNS FORMAT 


TITLE 

CCT HEADER DATA: 


1 


60 

C60 

NR EC 

NUMBER OF DATA RECORDS IN THIS "SECTION" ( = 8 ) 

61 


70 

no 

— 

BLANK CHARACTER 




1 

Cl 

SCENE 

NASA SCENE IDENTIFIER 


2 

— 

11 

CIO 

ASCDSC 

ASCENDING OS DESCENDING NODE 


12 

— 

21 

CIO 

WRSPTH 

WRS PATH NUMBER 


22 

— 

31 

no 

WRSROW 

WRS FOW NUMBER 


32 

- 

41 

no 

- 

BLANK CHARACTER 




1 

Cl 

SENSOR 

SENSOR IDENTIFIER 


2 

— 

11 

cn 

MISSION 

MISSION IDENTIFIER 


12 

- 

21 

CIO 

ORBIT 

SPACECRAFT ORBIT NUMBER 


22 

— 

31 

no 

ADSTAT 

ACTIVE DETECTOR STATUS 


32 

— 

41 

X10 

NUMDET 

NUMBER OF ACTIVE DETECTORS 


42 

— 

51 

no 

isrid 

WIDTH OF INEUT SPACE (PIXELS) 


52 

• 

61 

no 

— 

BLANK CHARACTER 




1 

Cl 

WRSLIN 

LINE NUMBER OW WRS FORMAT CENTER 

2 

- 

11 

no 

WRSPIX 

PIXEL NUMBER OF WRS FORMAT CENTER 

12 

- 

21 

no 

CP YEAR 

CENTER PICTURE EXPOSURE TIME | 

(YEAR) 

22 

— 

31 

no 

CPDAY 

CENTER PICTURE EXPOSURE TIME ! 

JULIAN DAY OF YEAR) 

32 

- 

41 

no 

CPHOUP 

CENTER PICTURE 'EXPOSURE TIME | 

HOUR) 

42 

— 

51 

no 

CPHIN 

CENTER PICTURE EXPOSURE TIME l 

MINUTES) 

52 

“ 

61 

F 10 

— 

BLANK CHARACTER 




1 

cl 

MFANN 

NUMBER OF MAJOR FRAMES OF ANNOTATION DATA 

2 

- 

11 

no 

M FAN C 

NUMBER OF MAJOR FRAMES OF ANCILLARY DATA 

12 

- 

21 

no 

GCFLG 

GEOMETRIC CORRECTION FLAG 


22 

— 

31 

cn 

GCDFLG 

GEOMETRIC CORRECTION DATA FLAG 

32 

— 

41 

CIO 

PCELG 

RADIOMETRIC CORRECTION FLAG 


42 

— 

51 

CIO 

ECDFLG 

RADIOMETRIC CORRECTION DATA FLAG 

52 

— 

61 

CIO 

— 

BLANK CHARACTER 




i 

Cl 

I SUN 

NUMBER CF MAJOR FRAMES (LINES) 

OF IMAGE DATA 

2 

- 

ii 

no 

IFORMAT 

IMAGE DATA FORMAT 


12 

— 

21 

CIO 

ILTYPE 

TYPE OF INTERLEAVING 


22 

— 

31 

CIO 

ILCNT 

LINE INTERLEAVING COUNT 


32 

— 

41 

no 

FSTYPE 

TYPE OF RESAMPLING APPLIED 


42 

— 

51 

CIO 

MAFECJ 

MAP PROJECTION APPLIED 


52 

“ 

61 

CIO 

— 

BLANK CHARACTER 




1 

Cl 

WRSOFF 

WRS OFFSET EROM FULLY PROCESSED IMAGE CENTER 

2 

- 

11 

no 

IDJUST 

IMAGE DATA JUSTIFICATION 


12 

- 

21 

CIO 

UBANDS 

NUMBEWR OF USABLE SPECTRAL EANDS 

22 

— 

31 

no 

MBNDNC 

MSS SPECTRAL BAND NUMBER 


32 

— 

41 

no 

— 

BLANK CHARACTER 




1 

Cl 

PASSTYPE 

TYPE OF PASS ( DAY / NIGHT ) 


2 

— 

11 

CIO 

CA1TYPE 

MODE CF CALIBRATION WEDGE 


12 

— 

21 

CIO 
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TABLE B-17. ANCILLARY LATA SET FIELDS ( CCT ANCILLARY RECORD NUMBER 1 DATA ) 


PARAMETER 

DESCRIPTION 


COLUMNS 

FORMAT 

TITLE 

NBEC 

CCT ANCILLARY RECORD NUMBER 1 DATA: 
NUMBER OF DATA RECORDS IN THIS "SECTION" 

( = 12 ) 

1 

61 

- 60 
- 70 

C60 

no 

ISWID 

ISLEN 

ISPACH 

ISPACV 

BLANK CHARACTER 

WIDTH OF INPUT SPACE (PIXELS) 

LENGTH OE INPUT SPACE (PIXELS) 

HOB. SPACING BETWEEN INPUT-SPACE PIXELS 
VEE. SPACING BETWEEN INPUT-SPACE PIXELS 

(METERS) 

(METERS) 

2 

12 

22 

32 

1 

- 11 
- 21 

- 31 

- 41 

Cl 
FI 0 
FI 0 
Fin 
F10 


BLANK CHARACTER 

OSWID WIDTH OF OUTPUT SPACE (PIXELS) 

OSLEN LENGTH OF OUTPUT SPACE JPIXELS) 

OSPACH HOE. SPACING BETWEEN OUTPUT-SPACE PIXELS 

OSPACV VER. SPACING BETWEEN OUTPUT-SPACE PIXELS 

SCALTM NOMINAL SPACECRAFT ALTITUDE (METERS) 

ISWIDM WIDTH OF INPUT SPACE (METERS) 


SD41 

SDU2 

SD43 

SD44 

SD45 

SD46 


(METERS) 

(METERS) 


BLANK CHARACTER 

SAMPLING-DELAY CORRECTION FOR BAND 4, 
SAMPLING-DELAY CORRECTION FOR BAND 4, 
SAMPLING-DELAY CORRECTION FCR BAND 4, 
SAMPLING-DELAY CORRECTION FOR BAND 4, 
SAMPLING-DELAY CORRECTION FOR BAND 4, 
SAMPLING-DELAY CORRECTION FOR BAND 4, 


2 

12 

22 

32 

42 

52 


1 

11 

21 

31 

41 

51 

61 


Cl 
FI 0 
F10 
F 1 0 
F10 
F 10 
Fl 0 



BLANK CHARACTER 




1 

Cl 

MIRRC 

CONSTANT COEFFICIENT OF MIRROR-VELOCITY 

MODEL 

2 

— 

21 

D2Q 

MIRE1 

DEGEEE-1 COEFFICIENT OF MIRROR-VELOCITY 

MODEL 

22 

— 

41 

D20 

MIRR2 

DEGFEE-2 COEFFICIENT OF MIRROR-VELOCITY 

MODEL 

42 


61 

D2C 

_ 

BLANK CHARACTER 




1 

Cl 

MIRR3 

DEGREE- 3 COEFFICIENT OF MIRROR-VELOCITY 

MODEL 

2 

— 

21 

D20 

MAXMIR 

MAXIMUM SCAN-MIRROR ANGLE (RADIANS) 


22 


41 

D 20 



BLANK CHARACTER 




1 

Cl 

SCSKEW 

SCAN-SKEW PARAMETER 


2 

— 

21 

D20 

TBTWN 

TIME BETWEEN CONSECUTIVE MIRROR SWEEPS 


22 

— 

41 

D20 

TACTIVE 

TIME FOR AN ACTIVE MIRROR SCAN (SECONDS) 


42 


61 

D20 

— 

BLANK CHARACTER 




1 

Cl 

A 

SEMIMAJOR AXIS OF EARTH-SPHEROID MODEL | 

[METER S) 

2 

- 

21 

D20 

B 

SEMIMINOR AXIS OF EARTH-SPHEROID MODEL ( 

METERS) 

22 

— 

41 

D20 

ECCON 

EARTH-CURVATURE CONSTANT 


42 


6n 

D20 



BLANK CHARACTER 





1 



1 

Cl 

SD11 

SAMPLING-DELAY 

CORRECTION 

FOR 

BAND 

1, 

DETECTOR 

2 

- 

11 

F1Q 

SD1 2 

SAMPLING-DELAY 

CORRECTION 

FOR 

BAND 

1, 

DETECTOR 

2 

12 

— 

21 

FI 0 

SD13 

SAMPLING-DELAY 

CORRECTION 

FOR 

BAND 

1, 

DETECTOR 

3 

22 

— 

31 

F10 

SD14 

SAMPLING-DELAY 

CORRECTION 

FOR 

BAND 

1, 

DETECTOR 

4 

32 

— 

41 

V IB 

SD15 

SAMPLING-DELAY 

CORREC TION 

FOR 

BAND 

1, 

DETECTOR 

5 

42 

— 

51 

FI 0 

SD 1 6 

SAMPLING-DELAY 

CORRECTION 

FOR 

BAND 

1/ 

DETECTOR 

6 

52 

— 

61 

F10 

« 

m 

« 

* 

* 

t 








a 

• 

• 


• 






1 

Cl 

DETECTOR 

1 

2 

— 

11 

F 10 

DETECTOR 

2 

12 

— 

21 

F1C 

DETECTOR 

3 

22 


31 

FI 0 

DETECTOR 

4 

32 

— 

41 

F 1C 

DETECTOR 

5 

42 

- 

51 

Fl C 

DETECTOR 

6 

52 


61 

Fl 0 





1 

Cl 

DETECTOR 

1 

2 


11 

F 1 0 

DETECTOR 

2 

12 

- 

21 

F 1 0 



22 


31 

FIG 



32 

— 

41 

F 10 



42 

— 

51 

Fl 0 



52 

— 

61 

Fin 



62 


71 

F 1 0 


BLANK CHARACTER 

SD51 SAMPLING-DELAY CORRECTION FOR BAND 5, 

SD52 SAMPLING-DELAY CORRECTION FOR BAND 5, 

BBOFF5 BAND-TO-BAND OFFSET FOR BAND 5 

BBOFF6 BAND-TO-BANE OFFSET FOR BAND 6 

BBOFF7 BAND-TO-BAND OFFSET FOR BAND 7 

BBOFF8A BAND-TO-BAND OFFSET FOR BAND 8A 

BBOFF8B BAND-TO-BAND OFFSET FOR BAND 8B 
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TABLE B-18 


ANCILLARY DATA SET FIELDS ( CCT ANCILLARY RECORD NUMBER 2 DATA ) 


PARAMETER DESCRIPTION COLUMNS FORMAT 


TITLE 

CCT ANCILLARY RECORD NUMBER 2 DATA: 

1 

— 

60 

C6Q 

NREC 

NUMBER OF DATA RECORDS IN 

THIS "SECTION" ( = 11 ) 

61 


70 

no 

— 

BLANK CHARACTER 




1 

Cl 

WRSFFM 

WHS FRAME 


2 

— 

11 

F1C 

WRSORB 

WRS ORBIT 


12 

— 

21 

F 10 

WRSLAT 

GEODETIC LATITUDE OF WRS FORMAT CENTER (RADIANS) 

22 

— 

31 

T?10 

WRSLON 

LONGITUDE OF WRS FORMAT CENTER (RADIANS) 

32 

— 

41 

F 10 

FCXCTIH 

SPACECRAFT TIME OF FRAME CENTER 

42 


61 

C20 

— 

BLANK CHARACTER 




1 

Cl 

NADLAT 

GEODETIC LATITUDE OF SPACECRAFT NADIR (RADIANS) 

2 

— 

21 

D20 

NADLON 

LONGITUDE OF SPACECRAFT NADIR (RADIANS) 

22 


41 

D20 

_ 

BLANK CHARACTER 




1 

Cl 

NADI RX 

EARTH-FIXED COORDINATE OF 

SPACECRAFT NADIR (METERS) 

2 

- 

21 

D20 

NADIRY 

EARTH-FIXED COORDINATE OF 

SPACECRAFT NADIR (METERS) 

22 

— 

41 

D20 

NADIRZ 

EARTH-FIXED COORDINATE OF 

SPACECRAFT NADIR (METERS) 

42 

— 

61 

D20 

— 

BLANK CHARACTER 




1 

Cl 

BETA 

HEADING ANGLE 


2 

- 

21 

D 20 

INADLIN 

INPUT-SPACE LINE NUMBER OF 

' SPACECRAFT NADIR 

22 

- 

41 

D20 

INADSFL 

INPUT-SPACE SAMPLE NUMBER 

OF SPACECRAFT NADIR 

42 

— 

61 

D20 

— 

BLANK CHARACTER 




1 

Cl 

DELVV 

NORMALIZED VELOCITY CHANGE OE SPACECRAFT 

2 

— 

21 

D20 

EPVEL 



22 

— 

41 

D20 

ERPARM 

EARTH-ROTATION PARAMETER 


42 

— 

61 

D20 

— 

BLANK CHARACTER 




1 

Cl 

rollC 

CONSTANT-TERM COEFFICIENT 

OF ROLL MODEL 

2 

— 

11 

F 10 

ROLL 1 

DEGEEE-1-TERM COEFFICIENT 

OF ROLL MODEL 

12 

— 

21 

T? 1 

ROLL 2 

DEGREE- 2-TERM COEFFICIENT 

OF ROLL MCDEL 

22 

— 

31 

F10 

PCLL3 

DEGREE-3 -TERM COEFFICIENT 

OF ROLL MODEL 

32 

— 

41 

Fin 

— 

BLANK CHARACTER 




1 

Cl 

PITCHP 

CONSTANT-TERM COEFFICIENT 

OF PITCH MODEL 

2 

— 

11 

no 

PITCH"? 

DEGEEE-1-TERM COEFFICIENT 

OF PITCH MODEL 

12 

— 

21 

Fin 

PITCH2 

DEG REE-2 -TERM COEFFICIENT 

OF PITCH MODEL 

22 

— 

31 

F10 

PITCH3 

DEGREE- 3- TERM COEFFICIENT 

OF PITCH MODEL 

32 

— 

41 

F10 

_ 

BLANK CHARACTER 




1 

Cl 

YAWC 

CONSTANT-TERM COEFFICIENT 

OF YAW MODEL 

2 

— 

11 

F 10 

Y AW 1 

DEGREE- 1-TEPM COEFFICIENT 

OF YAW MODEL 

12 

- 

21 

Fin 

YAW2 

DEGREE- 2- TERM COEFFICIENT 

OF YAW MODEL 

22 

- 

31 

F10 

YAW3 

DEGREE- 3 -TERM COEFFICIENT 

OF YAW MODEL 

32 

— 

41 

FIG 

- 

BLANK CHARACTER 




1 

Cl 

DELHP 

CONSTANT-TERM COEFFICIENT 

OF ALTITUDE -CHANGE MODEL 

2 

- 

11 

F‘1 0 

DELHI 

DEGREE- 1- TERM COEFFICIENT 

OF ALTITUDE-CHANGE MODEL 

12 

— 

21 

F10 

_ 

BLANK CHARACTER 




1 

Cl 

NUMGCP 

NUMBER OF GCPS USED TO GET 

' ERROR MODELS 

2 

— 

11 

Fin 
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TABLE B-19. ANCILLARY CATA SET FIELDS ( HRS POINTS FROM ANC. RECORDS 3 - 19 ) 


PARAMETER 

DESCRIPTION 

COLUMNS 

FORMAT 

TITLE 

HRS POINTS FROM ANCILLARY RECORDS 3-10: 

1 

- 

60 

C60 

NREC 

NUMBER OF DATA RECORDS IN THIS "SECTION” ( = K+1 ) 

61 

-* 

70 

110 

• 

BLANK CHARACTER 




1 

Cl 

HI 

HRS POINT NUMBER 

1 

2 

- 

5 

R4 

H2 

« 

HP S POINT NUMBER 

• 

2 

6 

• 

9 

R4 

m 

• 

HI 6 

• 

HRS POINT NUMBER 

16 

62 

• 

m 

65 

i 

— 

BLANK CHARACTER 




1 

Cl 

H17 

HRS POINT NUMBER 

17 

2 

— 

5 

E4 

H 18 

• 

HRS POINT NUMBER 

• 

18 

6 

• 

o 

F4 

« 

• 

H32 

• 

HRS POINT NUMBER 

32 

62 

• 

• 

65 

i. 

• 

• 



• 


• 

• 

• 

• 

• 



• 

« 


• 

• 


TABLE B-20 

. ANCILLARY DATA 

SET FIELDS ( VRS 

POINTS 

FROM ANC. RECORDS 

3 - 

1C ) 

PARAMETER 

DESCRIPTION 


COLUMNS 

FORMAT 

TITLE 

VRS POINTS FROM ANCILLARY RECORDS 

3-10 

: 1 


60 

C60 

NREC 

NUMBER OF DATA RECORDS IN THIS ‘'SECTION” 

( - L+1 ) 61 

— 

70 

1 10 

— 

BLANK CHARACTER 





1 

Cl 

VI 

VRS PCINT NUMBER 

1 


2 

— 

5 

R4 

H2 

• 

VRS POINT NUMBER 

• 

2 


6 

« 

9 

R 4 

• 

• 

vie 

♦ 

VRS POINT NUMBER 

16 


62 

• 

• 

65 

• 

R4 


BLANK CHARACTER 





1 

Cl 

VI 7 

VRS PCINT NUMBER 

17 


2 

- 

5 

R4 

VI 8 

• 

VRS POINT NUMBER 

• 

18 


6 

• 

9 

E4 

• 

V32 

• 

VRS PCINT NUMBER 

32 


62 

• 

• 

65 

• 

R4 

• 

• 




« 


* 

• 

• 

* 

• 




• 

• 


• 

• 
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B.2.8 


Shadeprint Data Set 


The shadeprint data set is a computer listing on which image data is displayed. 
Gray levels are simulated by overstrike printing. Each printed character posi- 
tion represents a single pixel. An example of a shadeprint data set is shown 
in Figure B-2. 
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Figure B-2. Example of a Shadeprint 


B.2.9 User-Request Data Set 


Each of the SLDMS programs has its own User-Request Data 
to allow a user to supply any optional or necessary data 


Set CURDS). It 
to the program. 


is used 


The general formats of all URDSs are similar. Each will consist of a FORTRAN 
NAMELIST followed by formatted data. The records on a URDS are 80-byte card 
images . 
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B.3 REFORMATTING PROGRAM (MRFT) 


B.3.1 Statement of Problem 

MSS data, as obtained from NASA or EDC, can be in any of several formats. The 
LARSYS System uses a particular Multispectral Image Storage Tape (MIST) format 
that facilitates multispectral processing. This format is defined on pages 5-52 
through 5-58 of the LARSYS System Manual. Reformatting, as defined here, is the 
reorganization of the image data into the MIST format. 

During the reformatting process, it is convenient to perform several other pre- 
liminary processing tasks. These include removing line-length-correction pixels, 
histogramming, striping reduction, and ancillary-data extraction. 


B.3. 2 Data Flow 

There is a severe limit on the amount of temporary disk storage available on the 
LARS computer. This restriction has influenced the way data will flow through 
the reformatting program. In particular, temporary image data sets will be 
stored on tape data sets. This will require somewhat more I/O than would be 
necessary on a system that had available about 35 million bytes of temporary 
disk space. In order to minimize tape-drive utilization, the X-format and 
BSQ-format tapes will be reformatted so as to use a maximum of three tape drives 
at one time. 


B.3. 3 Non-Image Data Considerations 

The non-image data that is contained on the Landsat MSS computer tapes will be 
used by other programs in the SLDMS. For X-format CCTs, the non- image data 
fields are completely defined (see "Generation and Physical Characteristics of 
the Landsat 1 and 2 MSS Computer Compatible Tapes," GSFC , November 1975, 
X-563-75-223) . The format of the BSQ and BIL tapes that will be produced by 
EROS Data Center are apparently not yet officially defined. These will be 
specified in the new edition of the Landsat Data User's Guide (expected to be 
published m January, 1979). Actual coding of some sections of the reformatting 
program cannot be started until those formats are defined. Currently, defini- 
tions of only X format and BSQ format (from MDP) are available to IBM. 


B.3. 4 Program Description 

The design of the Reformatting Program is described in Figure B-3 and the nine 
program specifications that flow. 
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FIGURE B-3. MODULE HIERARCHY FOR LANDS AT MSS REFORMATTING PROGRAM. 
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****** 


****** 


SAF/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME . MBFT 

DESIGN EF STEPHEN W. MDRPHRE Y 

FUNCTION THIS IS THE TOP-LEVEL DRIVER FOR THE LANDSAT MSS REFORMATTING 

EHOGBAM. 

INPUTS 1. AN IMAGE TAPE IN X, BSQ, OR BIL FORMAT. 

2. A DATA SET CONTAINING USER'S REQUESTS (ON CARD IMAGES) . 

OUTPUTS 1. AN IMAGE TAPE IN MIST FORMAT. 

2. AN SLDMS ANCILLARY DATA SET FOR THE SCENE TO BE PROCESSED. 

3. A PPINTER LISTING DATA SET. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... NOT APPLICABLE. 

PDL (PROGRAM OUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE 

READ USEE'S REQUESTS 

PRINT PARAMETERS DEFINING THIS RUN 

CASENTBY 

CASE 1 

REFORMAT X-FCRMAT DATA 
CASE 2 

REFORMAT BS Q-FORMAT DATA 
CASE 3 

REFORMAT BIL-FORMAT DATA 

FNECASE 

IF DATA IS UNCOFRECTED OR PARTIALLY-PROCESSED DATA THEN 
DETERMINE STRIPING-RED UCTION TABLES 
IF STRIPING REDUCTION HAS REQUESTED THEN 

PERFORM STRIPING REDUCTION 
ENDTF 
ENDIF 

PRINT EECG BAM-EXIT MESSAGE 

****** END OF SPECIFICATION **** ** 
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SAB /L AND SAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME MBFIBS 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION REFOEMATS IMAGE DATA IN BSQ FORMAT. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL MRFTBS ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

OPEN FIRST INPUT DATA SET (1,1) 

OPEN FIRST OUTPUT DATA SET JO, 1) 

READ AND PROCESS NON-IMAGE RECORDS 

WRITE ID RECORD 

DO EOR EACH IMAGE RECORD 

READ ONE INPUT RECORD (1,1) 

PROCESS INPUT DATA AND PLACE IT IN OUTPUT BUFFER 
INITIALIZE NON-IMAGE FIELDS IN OUTPUT BUFFER 
WRITE ONE OUTPUT RECORD (0,1) 

ENDDO 

WRITE END-OF-TAPE RECORD 
CLOSE DATA SETS (1,1) AND (0,1) 

DO FOR EACH REMAINING INPUT DATA SET (1,11) 

CPEN NEXT INPUT DATA SET (I,N) 

OPEN LAST OUTPUT DATA SET (0,N-1) FOR INPUT 
CPEN NEXT OUTPUT DATA SET (0,N) 

READ AND PROCESS NON-IMAGE RECORDS 
COPY ID RECORD FROM (O.N-1) TO (0,N) 

DO FOR EACH IMAGE RECORD 

READ ONE INPUT RECORD (I,N) 

READ ONE RECORD (O.N-1) 

PROCESS INPUT DATA AND PLACE IT IN OUTPUT BUFFER 
WRITE ONE OUTPUT RECORD (0,N) 

ENDDO 

WRITE END-OF-TAPE RECORD 

CLOSE DATA SETS (I,N) AND (0,N-1) AND (0,N) 
ENDDO 

****** END OF SPECIFICATION 
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****** 


SAR/LANDSAT DATS MERGING SYSTEM PROGRAM SPECIFICATION 


NAME MRFTXA 

DESIGNER STEPHEN W. MTJRPHREY 

FUNCTION READS AND PROCESSES X-FORMAT NON-IMAGE RECORDS. 

INPUTS 1. SEE CARLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEO. ... CALL MRFTXA { ) 


VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PD I ( EEOGRAM OUTLINE) 

READ -BYTE ID RECORD FROM INPUT DATA SET 
READ 624-BYTE ANNOTATION RECORD FROM INPUT DATA SET 
DETERMINE FORMAT LEVEL 
PRINT FIELDS FROM ID RECORD; 

SCENE/FRAME ID 
CCT SEQUENCE number 
DATA RECORD LENGTH 
IAT IDENTIFICATION 
MSS ADJUSTED LINE LENGTH = LLA 
MSS DATA MODE/CORRECTION CODE 
PRINT FIELDS FROM ANNOTATION RECORD: 

DATE OF ACQUISITION 

FORMAT CENTER LATITUDE AND LONGITUDE 
NADIR LATITUDE AND LONGITUDE 
ORBIT DIRECTION 
NOMINAL PATH AND ROW 

SENSOR AND SPECTRAL BAND SPECIFIC INFORMATION 

DIRECT OR RECORDED TRANSMISSION 

SUN ELEVATION AND AZIMUTH ANGLES 

TYPE OF COPRECTICNS CN DATA 

SCALE OE IMAGE 

MAP PROJECTION 

RESAMPLING ALGORITHM 

PREDICTED OR DEFINITIVE EPHEMERIS DATA 
PROCESSING PROCEDURE 

EARTH IMAGE OR RBV CALIBRATION IMAGE 
SENSOR GAIN OPTIONS 
TYPE OF MSS TRANSMISSION 
TYPE OF MSS DATA PROCESSING 
FRAME IDENTIFICATION NUMBER 
NDPF IDENTIFICATION CODE 
DIRECT OR RECORDED MSS DATA 
MSS DATA ACQUISITION SITE 
TICK-MARK LOCATION DATA 
IF THIS IS FROM FIRST STRIP THEN 

WRITE ID DATA TC ANCILLARY DATA SET 
WRITE ANNOTATION DATA TO ANCILLARY DATA SET 
END IF 

****** END OF SPECIFICATION 


ORIGINAL PAGE IS 
OF POOR QUALITY 


****** 


B-26 



****** 


****** 


SAP /LAN DSA 1 DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME MRFTX 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION REFORMATS X-FORHAT DATA. 

INPUIS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL MRETX ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

OPEN FIRST INPUT DATA SET (1,1) 

OPEN FIFST OUTPUT DATA SET jO,1) 

PE AD AND PROCESS NCN -IMAGE RECORDS 

WRITE ID RECORD 

DO FOR EACH IMAGE RECORD 

READ ONE INPUT RECORD (1,1) 

PROCESS INPUT DATA AND PLACE IT IN OUTPUT BUFFER 
INITIALIZE NON-IMAGE FIELDS IN OUTPUT BUFFER 
WRITE ONE OUTPUT RECORD (0,1) 

ENDDO 

WRITE END-OF-TAPE RECORD 
CLOSE DATA SETS (1,1) AND (0,1) 

DO FOR EACH REMAINING INPUT DATA SET (I, N) 

OPEN NEXT INPUT DATA SET (I,N) 

OPEN LAST OUTPUT DATA SET (0,N-1) FOR INPUT 
OPEN NEXT OUTPUT DATA SET JO.N) 

READ AND PROCESS NON-IMAGE RECORDS 
COPY ID RECORD FROM (O.N-1) TO (0,N) 

DO FOF EACH IMAGE RECORD 
READ ONE INPUT RECORD 
FEAD ONE RECORD (0,N-1) 

PROCESS INPUT DATA AND PLACE IT IN OUTPUT BUFFER 
W^ITE ONE OUTPUT RECORD (0,N) 

ENDDC 

WRITE END-OF-TAPE RECORD 

CLOSE DATA SETS (I,N) AND (0,N-1) AND (0,N) 

ENDDO 

READ AND PROCESS SIAT FILE 

****** END OF SPECIFICATION ****** 
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****** S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME MRFTBA 

DESIGNER STEPHEN W. MUBPHRE V 

FUNCTION READS AND PROCESSES NON-IMAGE RECORDS FROM BS Q-FORM AT OR 

BII-FCRMAT TAPE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL MRFTBSA ( } 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

DO UNTIL END-OF-DATA IS REACHED 
READ ONE NCN-IMAGE RECORD 
END DO 

IF THIS IS FIRST FILE OF INPUT TAPE THEN 
PFINT DIRECTORY RECORD 

WRITE DIRECTORY DATA TO ANCILLARY DATA SET 

EL SE 

PRINT TRAILER RECORD 

WRITE TRAILER DATA TO ANCILLARY DATA SET 
ENDIF 

PRINT HEADER RECORD 

PRINT ANCILLARY RECORDS 

PRINT ANNOTATION RECORDS 

WRITE HEADER DATA TO ANCILLARY DATA SET 

WRITE ANCILLARY DATA TO ANCILLARY DATA SET 

WRITE ANNOTATION DATA TO ANCILLARY DATA SET 

****** END OF SPECIFICATION 


** ** ** 


****** 
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SAF/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME MEFTBII 

DESIGNER STEPHEN W. MDEPHREY 

FUNCTION PROCESSES ONE LINE OF BSQ-FORMAT OR BIL-FORMAT IMAGE DATA. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL MSFTBI ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PD I (PROGRAM OUTLINE) 

MOVE IMAGE DATA INTO CURRENT-BAND AREA OF OUTPUT BUFFER 
IF THIS IS PARTIALLY-PROCESSED BATA THEN 
DO FOR EACH IMAGE PIXEL 

CBTAIN V = PIXEL VALUE 

INCREMENT V COUNTER FOR CURRENT DETECTOR BY 1 
ENDDO 
EN DIF 

****** END OF SPECIFICATION ****** 
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****** SAR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME MRFTBI 

DESIGNER STEPHEN W. MORPHREY 

FUNCTION REFORMATS BIL-FORMAT DATA. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL MRFTBI ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

OPEN INPUT DATA SET 

OPEN OUTPUT DATA SET 

READ AND PROCESS NON -IMAGE RECORDS 

WRITE ID RECORD 

DO FOR EACH IMAGE LINE 

INITIALIZE NON-IMAGE FIELDS IN OUTPUT BUFFER 
DO FOR EACH SPECTRAL EAND 
READ ONE INPUT RECORD 

PROCESS INPUT DATA AND PLACE IT IN OUTPUT BUFFER 
ENDDO 

WRITE ONE OUTPUT RECORD 
ENDDO 

WRITE END-OF-TAPE RECORD 

CLOSE INPUT AND OUTPUT DATA SETS 

****** END OF SPECIFICATION 


****** 


****** 
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B . 4 AUTOMATIC CONTROL-POINT LOCATION PROGRAM 


B.4.1 Statement of Problem 

The objective of this program is to determine the locations of a set of control 
points in a Landsat MSS image. These locations are pixel coordinates m the high- 
frequency- corrected input space (see Figure B-6 in Section B.6). They are to be 
used as inputs to the Geometric Transformation Program. 

In order to locate a particular control point in a given high- frequency-corrected 
input space, it is necessary to know the location of that control point in some 
high-frequency-corrected input space. A search area is defined to be a 128-by- 
128 pixel subimage of the given high-frequency-corrected input space. The con- 
trol point is assumed to be contained in the search area, but its precise 
location within the search area is not known. A window is defined to be a 
32-by-32 pixel subimage of a high-frequency-corrected input space. The control 
point is near the center of the window, and its location within the window is 
precisely known. 

The Automatic Control-Point Location Program uses cross correlation to compare 
the window with all possible 32-by-32 pixel subimages of the search area. This 
results in a 97-by-97 array of regularly spaced points on a correlation surface. 
The location of the peak of the correlation surface represents the location of 
the control point. Therefore it is the function of this program to locate the 
peak of the correlation surface for each control point. 

This program is required to determine control-point locations m an automatic 
process that requires no manual effort by the user. 


B.4.2 Data Flow 

The data flow through the Automatic Control-Point Location Program involves both 
image data and ancillary data. Windows and search areas are obtained and hori- 
zontally resampled to remove high-frequency errors, if necessary. Some ancil- 
lary parameters are used to determine locations of the search areas . The 
results from this program are added to the ancillary data set. 


B.4.3 Program Description 

The design of the Automatic Control-Point Location Program is described m 
Figure B-4, the eight program specifications that follow, and five of the pro- 
gram specifications from the Manual Control-Point Location Program. 
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4. MODULE HIERARCHY FOR AUTOMATIC CONTROL-POINT LOCATION PROGRAM. 
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****** SAR/LANDSAI data merging system program specification ****** 

NAME ACPI 

DESIGNEE STEPHEN W. MUEPHBEY 

FUNCTION THIS IS THE TOP-LEY EL DRIVES EOR THE AUTOMATIC CONTROL POINT 

LOCATION PROGRAM. 

INPUTS 1. CONTROL POINT LIERARY COPY TAPE OR MIST-EORMAT TAPE. 

2. AN SLDMS ANCILLARY DATA SET EOR THE SCENE TO BE PROCESSED. 

3. A DATA SET CONTAINING USER'S REQUESTS ( ON CARD IMAGES }. 

OUTPUTS 1. THE UPDATED SLDMS ANCILLARY DATA SET. 

2. A PRINTER LISTING DATA SET. 


COMMON ELOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... NOT APPLICABLE. 

PD I (PROGRAM CUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE 
READ USER'S REQUESTS 

IF CONTROL-POINT WINDOWS COME FROM CPL COPY TAPE THEN 

EXTRACT 32-3Y-32 PIXEL WINDOWS FROM CPL COPY TAPE INTO MEMORY 

ELSE 

RESAMPLE 32-BY-32 PIXEL WINDOWS FROM A MIST FORMAT TAPE INTO MEMORY 
EN DIF 

FIND INPUT-SPACE COORDINATES OF SEARCH AREAS 

RESAMPLE 128-BY-128 PIXEL SEARCH AREAS FROM THE MIST-FORMAT TAPE INTO MEMORY. 
OBTAIN 5-BY-5 CORRELATION SURFACE CENTERED AT NEAREST PIXEL TO EACH GCP. 
COEFFICIENTS OF DEGREE-4 LEAST SQUARES POLYNOMIAL FIT TO SURFACE OF EACH GCP. 
LOCATE PEAK OF EACH GCP SURFACE BY NEWTON'S METHOD. 

ESTIMATE QUALITY EACH GCP (VALUE OF SURFACE AND MINIMUM CURVATURE AT PEAK) . 
WRITE RESULTS TO ANCILLARY DATA SET 
PRINT PROGRAM- EXIT MESSAGE 

****** END OF SPECIFICATION ****** 
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****** SAR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME ACPLCT 

DESIGNER STEPHEN W. MDRPHRE Y 

FUNCTION EXTRACTS WINDOWS FROM A CONTROL POINT LIBRARY COPY TAPE 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL ACPLCT ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

SET $WRS AND $EOF FLAGS TO FALSE 
OPEN DIRECTORY FILE CF CPL COPY TAPE 
DO UNTIL $WFS IS TRUE OR $EOF IS TRUE 
READ ONE DIRECTORY RECORD INTO BUEEER 
IF A TAPE MARK WAS READ THEN 
SET $£OF FLAG TO TRUE 

PRINT MESSAGE THAT SUITABLE CONTROL POINTS WERE NOT FOUND 
ENDIF 

IF WRS ROW, WRS PATH. AND REGISTRATION TYPE ARE AS REQUESTED THEN 

IF A SPECIFIC LANDSAT MISSION WAS REQUESTED BY THE USER THEN 

IF LANDSAT MISSION IS AS REQUESTED BY THE USER THEN 
SET $WRS FLAG TO TRUE 
ENDIF 

ELSE 

SET $ WR S FLAG TO TRUE 
ENDIF 
ENDIF 
ENDDO 

CLOSE DIRECTORY FILE OF CPL COPY TAPE 
IF $WRS IS TRUE THEN 

PRINT INFORMATION FROM CURRENT DIRECTORY RECORD 
ENDIF 

OPEN CONTROL-POINT FIXE ON REEL 1 OF CPI COPY TAPE 
SET N = 1 

DO FOR EACH CONTPOL-POINT POINTER PI 
K = PI - N 
DO WHILE K > 0 

READ (SKIP) CNE CONTROL POINT RECORD 
IF END-OF-DATA IS REACHED THEN 

CLOSE CONTROL-POINT FILE ON CURRENT REEL OF CPL COPY TAPE 
OPEN CONTROL-POINT FILE CN NEXT REEL CF CPL COPY TAPE 
READ (SKIP) CNE CONTROL POINT RECORD 
ENDIF 
K = K - 1 
N = N + 1 
ENDDO 

READ ONE CONTROL POINT RECORD INTO MEMORY 
IF END-OF-DATA IS REACHED THEN 

CLOSE CONTROL-POINT FILE ON CURRENT REEL OF CPL COPY TAPE 
OPEN CONTROL-FOINT FILE ON NEXT REEL OF CPL COPY TAPE 
READ ONE CONTROL POINT RECORD INTO MEMORY 
ENDIF 
N = N + 1 

PRINT ANCILLARY DATA FROM CONTROL POINT RECORD 
ENDDO 

CLOSE CURRENT CONTROL-POINT FILE 

'****** END OF SPECIFICATION ****** 
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****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION **’ 

NAME ACPLSA 

DESIGNER STEPHEN N . MDRPHREY 

FUNCTION ESTIMAIFS THE LOCATIONS IN THE INPUT SPACE OF THE INPUT MAP 

CCNTRCI POINTS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

COMMON BLOCKS .. ERTHCONS, MSSCONS 
RESTRICTIONS ... 


CALLING SEQ. ... 


CALL ACPLSA (NATEST, ATTTIM, AMSROL, AMSPIT, 
AMSYAW , NALEST, ALTTIM, AMSALT, AMSNLD, AHSNLN, 
DELVV, RNOMLL, NMCP, RLDR, RLNS . ELEV, BETANR, 
XEC, Y EC, ZEC, VN, VAN, VCN, VELROT, VA , VC, 
CTPV, CTPH , ESTL, ESTS) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

NATEST 

1*2 

T 

— 

NUMBER OF AMS ATTITUDE VALUES 

ATTTIM 

P*8 

T 

NATEST 

TIME OF ATTITUDE VALUES 

ATTROL 

F;*8 

I 

NATEST 

AMS ROLL VALUES 

ATTPIT 

R*8 

I 

NATEST 

AMS ROLL VALUES 

ATTYAW 

R* 8 

I 

NATEST 

AMS ROLL VALUES 

NALEST 

T*2 

I 

- 

NUMBER OF ALTITUDE VALUES 

ALTTIM 

R*8 

I 

NALEST 

TIME OF ALTITUDE VALUES 

AMSALT 

Pi* 8 

I 

NATEST 

EPHEMERI S ALTITUDE VALUES 

AMSNLD 

R*8 

I 

- 

NADIR GEOCENTRIC LATITUDE 

AMSNLN 

R*8 

I 

— 

NADIR LONGITUDE 

DELVV 

R*8 

T 

— 

NORMALIZED S/C VELOCITY ERROR 

E.NCMLL 

R*8 

I 

— 

NOMINAL LINE LENGTH IN INPUT 

NMCP 

1*2 

T 

— 

NUMBER OF INPUT MAP CONTROL POINTS 

RLDR 

R*8 

I 

NMCP 

GEODETIC LATITUDE OF CP'S 

RLNR 

R* 8 

I 

NMCP 

LONGITUDE OF CP* S 

ELEV 

R* 8 

I 

NMCP 

ELEVATION OF CONTROL POINTS 

BETANR 

R*8 

I 

— 

S/C HEADING RNGLE 

XEC 

R* 8 

I 

NMCP 

EARTH CENTERED X COORDINATE OF CP 

YEC 

R*8 

T 

NMCP 

EARTH CENTERED Y COORDINATE OF CP 

ZEC 

R*8 

I 

NMCP 

EARTH CENTERED Z COORDINATE OF CP 

VN 

R* 8 

I 

— 

EARTH ROTATION VELOCITY OF NADIR 

VAN 

R*8 

I 

— 

ALONG TRACK COMPONENT OF VN 

VCN 

R*8 

I - 

- 

ACROSS TRACK COMPONENT OF VN 

VELROT 

R* 8 

T 

NMCP 

EARTH ROTATION VEL OF CP'S 

VA 

R*8 

I 

NMCP 

ALONG TRACK COMPONENT OF VELROT 

VC 

R*8 

x 

NMCP 

ACROSS TRACK COMPONENT OF VELROT 

CTPV 

R* 8 

T 

NMCP 

CTP VERTICAL COORDINATE OF CP'S 

CTPH 

R*8 

I 

NMCP 

CTP HORIZONTAL COORDINATE OF CP'S 

ESTL 

R* 8 

I 

NMCP 

ESTIMATED LINE COORDINATE OF CP'S 

ESTS 

R*8 

I 

NMCP 

ESTIMATED SAMPLE COORDINATE 


UNITS 


SECONDS 

RADIANS 

RADIANS 

RADIANS 

SECONDS 

METERS 

RADIANS 

RADIANS 

RADIANS 

SAMPLES 

RADIANS 

RADIANS 

METERS 

RADIANS 

METERS 

METERS 

METERS 

MET/SEC 

MET/SEC 

MET/SEC 

MET/SEC 

MET/SEC 

MET/SEC 

METERS 

METERS 

II LINE 

II SAMP 


PDL (PROGRAM OUTLINE) 


COMPUTE ATTITUDE / ALTITUDE MODELS FROM AMS DATA <ADTAAMOD> 

COMPUTE THE ANGLE BETA (S/C GROUND TRACK INCLINATION) <CPBETA> 
COMPUTE THE EARTH CENTERED COORDINATES OF THE NADIR <LLEC> 

COMPUTE EARTH- ROT ATICN VELOCITY VECTORS AT THE NADIR < VETHROT> 

MAP THE INPUT CONTROL POINTS FROM LAI/LON TC EARTH CENTERED (LLEC> 
MAP THE INPUT CONTROL POINTS FROM EARTH CENTERED TO CTP <ECCTP> 
COMPUTE EARTH-ROTATIONAL VELOCTITY VECTORS EOR THE CONTROL POINTS 
MAP THE INPUT CONTROL POINTS FROM CTP TO UTP <CTPUTP> 

MAP THE INPUT CONTROL POINTS FROM UTP TO II <UTPII> 


****** 


END OF SPECIFICATION 


****** 
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****** S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME ACPLRS 

DESIGNER STEPHEN W. MORPHREY 

FUNCTION RESAMPLES 32-BY-32 OR 128-BY-128 PIXEL AREAS FROM A MIST-FORMAT 

MSS TAPE. THE ONLY CHANGE RESULTING FROM THE RESAMPLING IS THE 
REMOVAL OF HIGH-FREQUENCY, HORIZONTAL GEOMETRIC ERRORS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL ACPLRS ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 

PDL (PROGRAM OUTLINE) 

EXTRACT ALL SUBIMAGES TO COMPUTER MEMORY (MCELEX) . 

DO FOR EACH EX TRACT El SUBIMAGE. 

COMPUTE GRID-POINT -CORRESPONDENCE AND OTHER CONSTANTS (MCPLGP) . 
RESAMPLE SUBIMAGE CORRECTING ONLY HIGH-FREQUENCY ERRORS (MCPLVR) , 
ENEDO 


****** 


END OF SPECIFICATION 


****** 



****** SAR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME ACP1CS 

DESIGNER STEPHEN W. MDRPHREY 

FUNCTION COMPUTES 97- BY-9V CORRELATION SURFACE, AND FINDS THE PEAK THAT 

DETERMINES THE 5-BY-5 CORRELATION SURFACE THAT WILL BE USED TO 
LOCATE THE CONTROL POINT. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL ACPLCS(C,A) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 

C E*8 I 25X1 5-BY-5 CORRELATION MATRIX 

PDL (PROGRAM OUTLINE) 

DO FOR EACH CONTROL ECINT. 

PLACE 32-BY-32 WINDOW IN THE CENTER OF A 128-BY-128 ARRAY. 

INITIALIZE ALL CTHER PIXELS IN THIS ARRAY TO ZERO. 

NORMALIZE WINDOW ARRAY BY SUBTRACTING WINDOW MEAN FROM EACH ELEMENT. 
CALCULATE THE VARIANCE OF THE WINDOW. 

CALCULATE 97-BY-97 ARRAY OF 32-BY-32 SEARCH-AREA VARIANCES. 

DO FOE EACH ROW OF 128-EY-128 SEARCH-AREA ARRAY. 

CALCULATE COOLEY-TUKEY FAST FOURIER TRANSFORM OF CURRENT ROW. 

ENDDO 

DO FOR EACH COLUMN OF 128-BY-128 TRANSFORMED SEARCH-AREA ARRAY. 

CALCULATE COOLEY-TUKEY FAST FOURIER TRANSFORM OF CURRENT COLUMN, 

ENDDO 

DO FOR EACH ROW OF 128-BY-128 WINDOW ARRAY. 

CALCULATE COOLEY-TUKEY FAST FOURIER TRANSFORM OF CURRENT ROW. 

ENDDO 

DO FOR EACH COLUMN OF 128-BY-128 TRANSFORMED WINDOW ARRAY. 

CALCULATE COOLEY-TUKEY EAST FOURIER TRANSFORM OF CURRENT COLUMN. 

ENDDO 

DOT PRODUCT: (SEARCH-AREA TRANSFORM) (COMPLEX-CONJUGATE WINDOW TRANSFORM) 

DO EOF EACH COLUMN OF 128-BY-128 PRODUCT ARRAY. 

CALC. COOLEY-TUKEY INVERSE FAST FOURIER TRANSFORM OF CURRENT COLUMN. 
ENDDO 

DO FOR EACH ROW OF 128-BY-128 TRANSFORMED PRODUCT ARRAY. 

CALCULATE COOLEY-TUKEY INVERSE FAST FOURIER TRANSFORM OF CURRENT ROW. 
ENDDO 

DO FOR CENTER 97-BY-97 ELEMENTS OF FINAL TRANSFORMED ARRAY. 

CORRELATION- SURFACE ELEMENT = RE AL**2 / SAVARIANCE * 16384. 

ENDDO 

SET C = 5- BY- 5 SUBMATRIX CENTERED AT MAXIMUM OF CORRELATION SURFACE. 

ENDDO 

****** END OF SPECIFICATION ****** 
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****** SAR/LANDSAT eat a merging system program specification ****** 

NAME ACPLSF 

DESIGNER STEPHEN W. MUFPHREY 

FUNCTION COMPUTES COEFFICIENTS OF A DEGREE 4 (15-TERM) LEAST-SQUARES 

POLYNOMIAL THAT FITS THE 5-BY-5 CORRELATION SURFACE. THIS IS 
DONE FOR EACH CONTROL °OINT. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 


RESTRICTIONS ... 


CALLING SEQ. 

• • • 

CALL 

ACPLSF (C 

,A) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

M 

R*8 

— 

25X15 

INTERNAL MATRIX 

- 

P 

R*8 

- 

15X25 

INTERNAL LEAST-SQUARES MATRIX 

— 

C 

R*8 

I 

25X1 

5-BY-5 CORRELATION MATRIX 

— 

A 

R*8 

0 

15X1 

POLYNOMIAL COEFFICIENT MATRIX 

- 


PD I (PROGRAM OUTLINE) 

K = 1 

DO FOR I = -2 TO 2 BY 1 


DO FOB J = -2 TO 2 BY 1 


M { 

;k,i ] 

= 1. 

M 1 

K,2 

= J 

M< 

:k,3 

= I 

Ml 

K,4 

= J**2 

M ( 

K,5 

= J * I 

M 1 

K,6 

= t**2 

M l 

K,7 

= j**3 

M ( 

K r8 

- j**2 * I 

M 1 

K,9 

- j * t**2 


K, 10 

= 1**3 

Ml 

K,11 

= j**4 

M i 

[K,12 

= j**3 * I 

Ml 

K, 13 

= j**2 * 1**2 

M 1 

K,14 

= J * 1**3 

M 1 

K,15 

= 1**4 

K 

= K - 

«- 1 

ENDDO 



ENDDO 



T 

-1 

T 


P = { M M ) M ( MATRIX EQUATION; P IS 15-BY-25 MATRIX ) 

DO FOR EACH CONTROL POINT. 

A = P C ( MATRIX EQUATION ) 

ENDDO 

****** END OE SPECIFICATION ****** 
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****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 


NAME ACPLSP 

DESIGNER STEPHEN W. MDRPHRE Y 

FUNCTION COMPUTES LOCATION OF THE CORRELATION-SURFACE PEAK USING 

NEWTON'S METHOD. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

CALLING SEQ. ... CALL ACPLSP (C,A) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

uni: 

NGCP 

1*4 

_ 


NUMBER OF CONTROL POINTS 


— 

F 

R*8 

- 

2X1 

1ST PARTI ALS 


-> 

S 

3*8 

- 

2X1 

ERACTICNAL CORRELATION ’ PEAK COORDS. 

— 

N 

R* 8 

— 

2X2 

2ND PARTIALS 


— 

IH 

R*8 

I 

NGCP 

INTEGER CORRELATION PEAK 

COORDINATE 


IV 

R*8 

I 

NGCP 

INTEGER CORRELATION PEAK 

COORDINATE 

- 

WH 

F*8 

I 

NGCP 

WINDOW COORDINATE OF GCP 

FEATURE 

— 

WV 

R* 8 

I 

NGCP 

WINDOW COORDINATE OF GCP 

FEATUPE 

- 

H 

R*8 

0 

NGCP 

INPUT-SPACE GCE COORDINATE 

— 

V 

3*8 

0 

NGCP 

INPUT-SPACE GCP COORDINATE 

— 


PDL (PROGRAM OUTLINE) 

DO FOB EACH CONTROL POINT L = 1 IC NGCP BY 1 
INITIALIZE X = Y = 0 . 

DO FIVE TIMES. 

2 2 3 2 2 3 

F (1) = A +2A X +A Y +3A X +2A XY +A Y + 4A X +3A X Y+2A XY +A Y 
1 3 4 6 7 8 10 11 12 13 


2 2 3 2 2 3 

F (2) = A +2 A X+A Y + 3A X + 2A XY+A Y +4A X +3A X Y +2A XY +A Y 

2 4 5 7 8 9 11 12 13 14 

2 2 

N (2,2) = DFlDX = 2A + 6A X + 2A Y + 12A X + 6A XY + 12A Y 

367 10 11 12 

2 2 

N (1,2) = -DF1DY = - A - 2A X - 2A Y - 3A X - 4A XY - 3A Y 

4 7 8 11 12 13 

N(2 r 1) = -DE2DX = -DEI DY 

2 2 

N ( 1 , 1) = DF2DY = 2A + 2A X + 6A Y + 2A X + 6A XY + 12A Y 

5 8 9 12 13 14 

D = 1. / ( DFlDX * DF2DY - DF1DY * DF2DX ) 

'S = S - D N F ( MATRIX EQUATION ) 

x = s r 

Y = S (2) 

END DO 

H (I) = S (1) + IH (L) + WH (L) 

V(L) - S(2) + IV (L) + WV(L) 

* END OF SPECIEICATION ****** 
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****** SAB/LANDSAT data merging system program specification 

NAME ACPLPQ 

DESIGNER STEPHEN W. MDRPHREY 

FUNCTION COMPUTES ESTIMATES OF THE QUALITY OF THE CORRELATION- 

PEAKS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS 
CALLING SEQ. . 


CALL ACPLPQ ( ) 


VARIABLE 

TYPE 

I/O 

DIM 

NGCP 

1*4 

— 

— 

MAXCOR 

R*8 

— 

NGCP 

MINCURV 

B*8 

-* 

NGCP 

N 

P*8 

I 

2X2 

S 

R*8 

I 

2X2 


DESCRIPTION 


UNITS 


PDL (PROGRAM OUTLINE) 

DO FOR EACH CONTROL POINT 
E 
c 


= DF1DX = N(2,2) 

= DF1DY = -N (1 ,2) 

T = DF2DY * N (1 , 1) 

U = SQRT ( (R-T)**2 + (2*S)**2 ) 
MINCURV = ( R + T + U ) / 2.0 


NUMBER OF CONTROL POINTS 
VALUE OF SURFACE AT PEAK 
MINIMUM CURVATURE AT SURFACE PEAK 
2ND PARTI ALS 

FRACTIONAL CORRELATION PEAK COORDS. - 


L ** 1 TO NGCP BY 1 


: m\ 

MAXCOR = A 


X = 
Y 


A Y 
9 

PRINT RESULT 
ENDDO 


A X + A Y 
1 2 


+ A X 

10 


+ A X + 
3 

3 

A X Y + 
11 


A XY + 
4 


A Y 
5 


+ A X 
6 


A X Y + A 
7 l 


2 2 34 

A XY +A XY +A i 
12 13 14 


****** 


END OF SPECIFICATION 


****** 


SURFACE 


2 

XY + 


****** 
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B.5 MANUAL CONTROL-POINT LOCATION PROGRAM 


B.5.1 Statement of Problem 

The objective of this program is to produce digital image data sets in which the 
scale has been, changed or horizontal, high-frequency geometric errors have been 
removed. These two geometric changes are optional, and no other kinds of geo- 
metric changes will be performed. 

The geometric changes are produced by resampling the digital image data using 
cubic convolution. The concept of resampling is discussed in detail in Section 
B.7. The algorithm described there is also used in this program. The scale 
change is required to be a digital enlargement by an integral factor (input- 
pixel spacing is an integral multiple of output-pixel spacing) . This simplifies 
the resampling process and guarantees that output-pixel coordinates can easily 
be converted to input-pixel coordinates . 


B.5. 2 Data Flow 

One pass through an input MIST-format data set is made by this program. During 
this pass, all subimages are extracted and stored in the computer's (virtual) 
memory. This obviously places some limit on the amount of data that may be pro- 
cessed. The program is designed to process a maximum of 100 subimages. The 
buffer that contains the input subimage data is obtained dynamically at execu- 
tion time. Therefore, the effective maximum number of subiraages depends on the 
amount of virtual memory available and the size of the subimages being processed. 

It will be fewer than 100 in some cases . 

If histograms of the subimages have been requested by the user, input (unresampled) 
data will be used. These histograms may be used to generate character sets for 
the line-printer displays (shadeprints) of the subimages. 


B.5. 3 User-Specified Inputs to MCPL 

The user- specified input parameters for MCPL are contained m a data set consist- 
ing of 80-byte records. This data set contains the following: 

a. A FORTRAN NAMELIST 

b. A shadeprmt character set (optional) 

c. One or more Sub image Descriptor Records. 

These are described in Tables B-21, B-22, and B-23. 
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k TAEL E E- 2 1 . MCPL NAMELIST 


PARAMETER 


DE SCRIPTION 


UNITS DEFAULT 


NWID 

NLEN 

NSTBK 

BANDS 

$HIST 

$ R DC HE 

SSHDPT 

SRSFL 

$HFCOR 

$1 SUBW 

$OSUBW 

HSPACO 

VS PA CO 

HSPACI 

VSPACI 

HSCFAC 

VSCFAC 


WIDTH OF A SUBIMAGE 

LENGTH OF A SUBIMAGE 

NUMBER OE OVERSTRIKES TO BE USED FCB SHADEPRINTS 

LIST OF SPECTRAL BANDS TO BE PROCESSED 

INDICATES THAT HISTOGRAMS OF SUBIMAGES ARE WANTED 

INDICATES THAT USER WILL SUPPLY SHADEPRINT CHARACTER 
SET 

INDICATES THAT SHADEPRINTS 0? SUBIMAGES ARE WANTED 

INDICATES THAT RESAMPLING WILL EE PERFORMED ON EACH 
SUBIMAGE 

INDICATES THAT HIGH-FREQUENCY HORIZONTAL ERRORS 
WILL BE CORRECTED 

INDICATES THAT A TAPE COPY OF UNCOPRECTED SUBIMAGES 
WILL BE MADE 

INDICATES THAT A TAPE COPY OF RESAMPLED SUBIMAGES 
WILL BE MADE 

HORIZONTAL SPACING BETWEEN PIXEL CENTERS OF 
RESAMPLED SUBIMAGE 

VERTICAL SPACING BETWEEN PIXEL CENTERS OF 
RESAMPLED SUBIMAGE 

HORIZONTAL SPACING BETWEEN PIXEL CENTERS OF 
INPUT IMAGE 

VERTICAL SPACING BETWEEN PIXEL CENTERS OF 
INPUT IMAGE 

HORIZONTAL SCALE FACTOR = HSP ACI/HSPACO 

VERTICAL SCALE FACTOR = VSPACI/VSPACO 


PIXELS 

PIXELS 

4 

- 1,2,3 


F 

F 

F 


METERS 

METERS 

METERS 

METERS 
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TABLE B-22. SHADEPRINT CHARACTER SET 


ORIGINAL PAGE IS 
OF POOR QUALITY 


PARA METE B 

DESCRIPTION 



COLUMNS 

FORMAT 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 1ST OVERSTRIKE 

VALUES 

0-63 

1 

- 

64 

64C1 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 1ST OVERSTRIKE 

VALUES 

64 - 127 

1 

- 

64 

64C1 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 1ST OVERSTRIKE 

VALUES 

128 - 191 

1 

- 

64 

64C1 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 1ST OVERSTRIKE 

VALUES 

192 - 255 

1 

- 

64 

64C1 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 2ND OVERSTRIKE 

VALUES 

0-63 

1 

- 

64 

64C1 

$CHAP 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 2ND OVERSTRIKE 

VALUES 

64 - 127 

1 

- 

64 

64C1 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 2ND OVERSTRIKE 

VALUES 

128 - 191 

1 

- 

64 

64C 1 

$CHAR 

PRINTABLE CHARACTERS FOR PIXEL 
USED FOR 2ND OVERSTRIKE 

VALUES 

192 - 255 

1 

- 

64 

64C1 

• 

• 

« 

• 

m 

• 




• 

• 

• 


• 

♦ 


TABLE B-23. SUBIMAGF EESCRIPTOR RECORD 


PARAMETER 

DESCRIPTION 

COLUMNS 

FORMAT 

CENSAM 

INPUT-SPACE HORIZONTAL COORDINATE OF CENTER PIXEL 
OF SUBIMAGE 

1 - 5 

14 

CENLIN 

INPUT-SPACE VERTICAL COORDINATE OF CENTER PIXEL 
OF SUBIMAGE 

6-10 

14 

WIDTH 

WIDTH OF SUBIMAGE IE PIXELS 

11 - 15 

14 

LENGTH 

LENGTH OF SUBIMAGE IN PIXELS 

16-20 

14 

BANDS 

LIST OF SPECTRAL BAND NUMBERS 

21-25 

5X1 

TITLE 

PRINTABLE TITLE FOR THE SUBIMAGE 

26 - 17 . 

47C1 
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The NAMELIST allows the user to define the functions to be performed and to 
specify certain constants. The shadeprint character set may be specified by 
the user, if desired. There must be a Subimage Descriptor Record for each sub- 
image to be processed. 


B.5.4 Program Description 

The design of the Manual Control-Point Location Program is described in Fig- 
ure B-5 and the program specifications that follow. 
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MCPLRD 

MCPLEX 

MCPLHI 

MCPLAH 

MCPLPH 

MCPLCS 

HCPLGP 

EiCPLVR 

MCPLWR 

MCPLPS 




IS 


OF POOR 


MCPLHF 

HCPLHR 


FIGURE B-5. MODULE HIERARCHY FOR MANUAL CONT EOL-POINT LOCATION PROGRAM. 
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****** 


****** 


SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME MCPL 

DESIGNER STEPREN W. MDRPHREY 

FUNCTION THIS IS THE TOP LEVEL DRIVER FOR THE MANUAL CONTROL -POINT 

LOCATION PROGRAM. 

INPUTS 1. A DATA SET CONTAINING USER'S REQUESTS (ON CARD IMAGES). 

2. A MIST-FORMAT DATA SET 

OUTPUTS 1. A ERIN TER LISTING DATA SET. 

2. MIST-FORMAT TAPE(S) CONTAINING EXTRACTED SUBIMAGES. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... NOT APPLICABLE. 

PDL (PROGRAM OUTLINE) 

PRINT SEPARATOR PAGE TO INDICATE START OF RUN. 

READ USER'S REQUESTS AND INITIALIZE CONSTANTS AND ARRAYS. 

EXTRACT AIL SUBIMAGES TO COMPUTER MEMORY. 

IF TAPE COPY OF EXTRACTED DATA NAS REQUESTED THEN 
WRITE ALL EXTRACTED SUBIMAGES TO TAPE 
IN DIF 

DO FOR EACH EXTRACTED SUBIMAGE. 

IF HISTOGRAM WAS REQUESTED THEN 
ACCUMULATE HISTOGRAM DATA. 

PRINT HISTOGRAM STATISTICS. 

ENDIE 

IF SHADEPRINT WAS REQUESTED AND CHARACTER SET NOT GIVEN THEN 
CREATE A CHARACTER SET FROM THE HISTOGRAM. 

END IF 

IF RESAMPLING WAS REQUESTED THEN 

COMPUTE GRID-POINT CORRESPONDENCE AND OTHER CONSTANTS. 

RESAMPLE SUBIMAGE {CHANGING SCALE OR CORRECTING HIGH-FREQUENCY ERRORS) 
IF TAPE COPY OF RESAMPLED DATA WAS REQUESTED THEN 
WRITE RESAMPLED DATA SET TO TAPE. 

END IF 
ENDIF 

IF SHADEPRINT WAS REQUESTED THEN 

DO FOR EACH REQUESTED SPECTRAL BAND 
PRINT THE SHADEPRINT. 

ENDDO 

ENDIF 

ENDDO 

PRINT SEPARATOR PAGE TO INDICATE END OF RUN. 

****** END OF SPECIFICATION ****** 
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****** SA.R/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME MCPLEX 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION EXTRACTS ALL SUBIMAGES FROM INPUT MIST-FORMAT DATA SET 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


CALLING SEQ. ... CALL MCPLEX (CHGFLG ,CHGLIN, SUEPTR, LINPTR, MAXLIN ,T BYTES, NUMSUB) 

VARIABLE TYPE I/O DIM DESCRIPTIION UNITS 

CHGFLG 1*2 0 2X200 SUBIMAGE INDEX / CHANGE TYPE 

CHGLIN 1*4 0 200 LINE NUMBER OF CORRESP. CHANGE PIXELS 

SUBPTS 1*4 0 100 ADDRESS OF SOBIMAGE IN BUFFER 

LINPTR 1*4 0 100 ADDRESS OF NEXT LINE IN BUFFER 

MAXLIN 1*4 0 - NUMBER OF INPUT RECORDS TO USE 

TBYTES 1*4 0 - SIZE OF EUFFER BYTES 

NUMSUB 1*4 0 - NUMBER OF SUBIMAGES TO EXTRACT 

PDL (PROGRAM OUTLINE) 

OBTAIN VIRTUAL MEMORY FOR SUBIMAGE BUFFER ( GETMAIN ). 

DO FOR EACH SUBIHAGE ( N = 1. NUMSUB ) 

SUBPTR (N) = SUBPTR (N) + ADDRESS OF BUFFER 
LINPTR (NJ = SUBPTR (N) 

ENDDO 

OPEN INPUT MIST-FORMAT DATA SET 
BEAD CCNTROL RECORD 
SET LINE#, NLIVE TO 0 
DO WHILE LINE# < MAXLIN 
LINE# = LINE# + 1 
READ ONE IMAGE RECORD 

UPDATE ARRAYS THAT CONTROL PROCESSING CHANGES 

DO WHILE PROCESS CHANGES AT THIS LINE 

IF A SUBIMAGE BEGINS THEN 
NLIVE = NLIVE + 1 

COPY ELEMENT FROM CHGFLG ARRAY TO BOTTOM OF SEQLST ARRAY 

ELSE 

FIND ELEMENT OF SEQLST ARRAY THAT IS NO LONGER NEEDED. 

COPY BOTTOM SEQLST ELEMENT INTO VACANT SPOT. 

NLIVE = NLIVE - 1 
ENDIF 

INCREMENT RELATIVE POINTER TC NEXT ELEMENT OF CHANGE ARRAYS. 

K = NLIVE 
DO WHILE K > 0 

MOVE LINE NUMBER INTO SUEIMAGE BUFFER { 2 BYTES ) . 

MOVE 2 BYTES OF ZERO FILL INTO SUBIMAGE BUFFER. 

DC FOR EACH SPECTRAL BAND 

HOVE IMAGE DATA INTO SUBIMAGE BUFFER. 

MOVE 6 BYTES OF ZERO FILL INTO SUBIMAGE BUFFER. 

ENDDO 

MOVE INPUT IMAGE DATA INTO BUFFER UNDER CONTROL OF ARRAYS 
ENDDO 

****** END OF SPECIFICATION ****** 
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****** 


SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME MCPLKD 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION READS USER'S REQUESTS AND INITIALIZES CONSTANTS AND ARRAYS 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET, 


CALLING SEQ. ... CALL MCPLRD (CHGFLG ,CH GLIN. BAND S , LFTS AM.LINPTR, LENGT H,W IDTH, 

RECLEN,SUBPTR,TOPLIN r $TITLE, DEFALT, SFLAGS) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 

CHGFLG 1*2 0 2X200 SUBIHAGE INDEX / CHANGE TYPE 

CHGLIN 1*4 0 200 LINE NUMBER OF CCRRESP. CHANGE PIXELS 

BANDS 1*4 0 100 EACH BIT IS FLAG FOR THAT BAND 

LFTS AM 1*4 0 100 COORDINATE OF LEFT SAMPLE PIXELS 

LINFTB 1*4 O 100 ADDRESS CF NEXT IINE IN BUFFER 

LENGTH 1*4 0 100 VERTICAL SIZE OF SUBIHAGE PIXELS 

WIDTH 1*4 0 100 HORIZONTAL SIZE OF SUBIMAGE PIXELS 

SECXEN 1*4 O 100 LENGTH OF ONE SUBIMAGE RECORD BYTES 

SUBPTR 1*4 0 100 ADDRESS OF SUBIMAGE IN BUFFER 

TOPLIN 1*4 0 100 COORDINATE OF TCP LINE PIXELS 

STITLE L*1 0 47X100 TITLE OF SUBIMAGE 

DEFALT 1*4 0 10 VARIOUS DEFAULT VALUES . 

$FLAGS L*4 0 7 VARIOUS PROCESSING FLAGS 

PDL (PROGRAM OUTLINE) 

. SET DEFAULT VALUES FCR USER'S OPTIONS. 

' READ USER'S OPTIONS 

IF SHADEPRINTS WERE REQUESTED THEN 

COUNT THE (DEFAULT) NUMBER OF SPECTRAL BANDS TO BE SHADEPRINTED 
IF SHADEPRINT CHARACTER SET IS SUPPLIED BY THE USER THEN 
* READ THE CHARACTER SET. 

PRINT THE CHARACTER SET. 

END IF 
ENDIF 

INITIALIZE N t MAXLIN r SUBPTR (11, LINPTR{1) TO 0 
DO UNTIL END-OF-DATA IS REACHED ON CARD-IMAGE INPUT DATA SET 
N = N + 1 

READ ONE SUBIMAGE-DESCRIPTOR CARD. 

IF END-OF-DATA WAS NOT REACHED THEN 

IF WIDTH (N), LENGTHEN). OR BANDS (N) WERE NOT GIVEN THEN 

USE THE APPROPRIATE DEFAULT VALUE (S) FOR THOSE NOT GIVEN 
ENDIF 

SET SPECTRAL-BAND FLAGS (IN ARRAY EIEMENT) ACCORDING TO BAND LIST. 
TOPLIN (N) = MAXIMUM ( 1 , CENLIN - LENGTH (N)/2 ) 

LFTSAMjN = MAXIMUM ( 1 , CENLIN - WIDTH (N) /2 ) 

SET BOTH INDICES IN CHGFLG ARRAY TO 4 * N - 4 
CHGLIN (1,N) = TOPLIN (N) 

CHGLIN (1,N) = TOPLIN (N) + LENGTH (N) 

RECLEN(N) = 4 + ( WIDTH (N) + 6 ) * N BANDS 
SUBPTR iN + 1) = SUBPTR (N) + RECLEN(N) * LENGTH (N) 

LINPTP (N+-1) « SUBPTR (N+1) 

MAXLIN = MAXIMUM ( MAXLIN , TOPLIN (N) + LENGTH (N) - 1 ) 

ENDIF 
N = N - 1 
ENDDO 

TOTAL BYTES IN BUFFER = SUBPTR (N + 1) 

SORT BOTH CHANGE ARRAYS SO THAT THEY ARE IN ASCENDING ORDER OF CHGLIN ARRAY. 
****** END OF SPECIFICATION ****** 


B-50 



ORIGINAL PAGE IS 
OF POOR QUALITY 


****** SAF./LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME MCPLWI 

DESIGNER STEPBEN W. MURPHREY 

FUNCTION WRITES ALL EXTRACTED SUBIMAGES TO A SINGLE COMPUTER TAPE 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 

CALLING SEQ. ... CALL MCPLWI (BANDS, LENGTH ,RECIEN, SUBPTR, IDREC) 

VARIABLE TYPE I/O DIM DESCRIPTIION UNITS 

BANDS 1*4 I 100 EACH BIT IS FLAG FOR THAT BAND 

LENGTH 1*4 I 100 VERTICAL SIZE OF SUBIMAGE PIXELS 

RECLEN 1*4 I IOC LENGTH OF ONE SUBIMAGE RECORD BYTES 

SUBPTR 1*4 I 100 ADDRESS OF SUBIMAGE IN BUFFER 

IDREC 1*4 I 2C0 ID RECORD FROM INPUT DATA SET 


PDL (PROGRAM OUTLINE) 

INITIALIZE FIELDS IN ID-RECORD BUTEER: 

RUN NUMBER 
CONTINUATION CODE 
FLIGHT-LINE ID 

MONTH, DAY, YEAR, & TIME DATA WAS TAKEN 
ALTITUDE OF PLATFORM 
HEADING OF VEHICLE 
DATE OF THIS TAPE 

DO FOR EACH EXTRACTED SUEIMAGE { I = 1,NUMSUB ). 

OPEN OUTPUT TAPE DATA SET FOR FILE I. 

FILL IN REMAINING FIELDS IN ID-RECOFD BUFFER; 

LARS TAPE NUMBER 
FILE NUMBER 

NUMBER OF DATA CHANNELS 
NUMBER OF SAMPLES PER CHANNEL 
NUMBER OF LINES IN THE IMAGE 
WAVELENGTH RANGE OF EACH SPECTRAL BAND 
WRITE THE ID RECORD ONTO THE TAPE. 

OBTAIN THE ADDRESS OF THE SUBI MAGE DATA FROM THE SUBPTR ARPAY. 
OBTAIN THE WIDTH OF ONE RECORD FROM THE RECLEN ARRAY. 

DO FOR EACH LINE IN THE SUBIMAGE. 

WRITE THE CURRENT IMAGE RECORD ONTO THE TAPE. 

INCREMENT ADDRESS OE SUBIMAGE DATA BY THE WIDTH OF ONE RECORD. 
ENDDO. 

CLOSE THE OUTPUT TAPE DATA SET. 

ENDDC. 

OPEN THE OUTPUT TAPE DATA SET FOR FILE NUMSUB + 1. 

SET FIELDS IN THE END-OF-T APE-RECORD BUEFERS 
LARS TAPE NUMBER 
FILE NUMBER 
CONTINUATION CODE 

WRITE THE END-CF-TAPE RECORD ONTO THE TAPE. 

CLOSE THE OUTPUT TAPE DATA SET. 

****** END OF SPECIFICATION 


****** 


** 


B-51 



****** SIB /LAND SAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME MCPLAH 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION , , ACCUMULATES HISTOGPAM DATA FOR A SINGLE EXTRACTED SUBIMAGE . 

INPUTS 1. SEE CABLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 

CALLING SEQ. ... CALL MCPLAH (B ANDS , LENGTH ,RECI-EN, SUBPTR, HIST) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 

EANDS 1*4 I 100 EACH BIT IS FLAG FOR THAT BAND 

LENGTH 1*4 I 10O VERTICAL SIZE OF SUBIMAGE PIXELS 

RECLEN 1*4 I 100 LENGTH OF ONE SUBIMAGE RECORD BYTES 

SUBPTR 1*4 I 100 ADDRESS OF SUBIMAGE IN BUFFER 

HIST 1*4 0 256X5 HISTOGRAM ARRAY 


PDL (PROGRAM OUTLINE) 

ZERO THE HISTOGRAM ARRAY. 

INITIALIZE THE DATA POINTER. 

DO FOR EACH LINE IN THE SUBIMAGE. 

INITIALIZE POINTER TO THE BAND-1 HISTOGRAM ARRAY. 

INITIALIZE POINTER TO THE BAND-1 DATA ARRAY. 

DO FOR EACH SPECTRAL BAND. 

DO FOR EACH PIXEL IN THE LINE. 

OBTAIN THE PIXEL VALUE V. 

INCREMENT HISTOGRAM-ARRAY-ELEMENT V BY 1 . 

ENDDO. 

INCREMENT HISTOGRAM POINTER TO THE NEXT BAND 
INCREMENT DATA POINTER TO THE NEXT BAND 
ENDDO. 

INCREMENT DATA POINTER TO THE NEXT LINE. 

ENDDO. 

****** END OF SPECIFICATION ****** 
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****** SAR/LAND sat data merging system program specification ****** 

NAME MCPLPH 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION PRINTS HISTOGRAMS AND THE RESULTING STATISTICAL PARAMETERS FOR 

ALL BANDS OF A SINGLE EXTRACTED SUBIMAGE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING BATA SET. 

CALLING SEQ. ... CALL MCPLPH (BANDS , HIST, CUMPCT) 

VARIABLE TYPE I/O DIM DESCRIPTITON UNITS 

BANDS 1*4 I ICO EACH BIT IS FLAG FOR THAT BAND 

HIST 1*4 I 256X5 HISTOGRAM 

CUMPCT F*4 0 256X5 CUMULATIVE HISTOGRAM PERCENT 


PDL (PROGRAM OUTLINE) 


DO FOR EACH SPECTRAL BAND. 

FIND SMALLEST PIXEL VALUE THAT OCCURS IN THIS BAND. 
FIND LARGEST PIXEL VALUE THAT OCCURS IN THIS BAND. 
COMPUTE THE TOTAL NUMBER OF PIXELS IN THIS BAND. 
INITIALIZE THE CUMPCT ARRAY TO ZEROS. 

INITIALIZE MENTOT, SMSTOT, S C ARRAYS TO ZEROS. 
ENDDO 

DO FOR EACH SPECTRAL EAND ( B ) . 

DO FROM K = 1 TO 256 BY 1 

HIS PCT (K, B) = HIST (K, B) * 100. / TOTAL (B) 

C (B) = C (B) + HI ST (K ,B) 


(K,B) 


CiJm(K,B) 

CUMPCT JK,B) = 
IF HISPCT(K,B) 
MODE (B) = 

ENDIF 

IF CUMPCT (K,B) 
MEDIAN (B) = 
ENDIF 
MENTOT (B) 

RMSTOT (B) 


ijli + H 

:tJM (K,B) 
> HISPC 


(K,B) 


* 100 . 

T (MODE ( 

THEN 


. / TOTAL' 


RMSTOT 
ENDDO 
MEAN (B) 

RMS (B) 
STDEV(B) 
MEDIAN (B) 
MODE(B) 

SUM = 0 
DO FROM K 
SUM 
COUNT 
L 

PRINT 


MENTOT 

RMSTOT 


T (B) + (K-1) 
T (B) + (K-1) 


* HIST 
**2 * H 




B) **2 


= MENTOT (B) / TOTAL (B) 
= SQRT ( RMSTOT (B) / 

= SQRT ( RMS(B)**2 - 
= MEDIAN (B) - 1 

= MODE (B) - 1 

= MINIMUM VALUE TO M 
= SUM + HIST (K, B) 

= MAXIMUM ( 1 , HIST 
= MAXIMUM { 1 , 10. 


/ TOTAL (B) ) 
- MEAN (B) **2 


MAXIMUM VALUE 


j 1 10f'* (K £ol;i!) (COUNT) 

T(K,B) , SUM . HISPCT ( K, 


K-1 , HIST (K ,B) , SUM . HISPC 

LINES (1,L) , LINES (2 ,L) , ... 

NUMBANDS , MEAN (B) , STDE V (B) 


RMS (B) , MEDIAN (B) 


****** 


END OF SPECIFICATION 


****** 
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SAB /LAND SAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


•Jesfe 


NAME MCPLCS 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION CREATES CHARACTER SETS (BASED ON HISTOGRAMS) 

A SINGLE EXTRACTED SUBIMAGE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 

CALLING SEQ. ... CALL MCPLCS (BANDS, $CH AH, CUMPCT) 

VARIABLE TYPE I/O DIM DESCRIPTION 

BANDS 1*4 I 100 EACH BIT IS FLAG FOR THAT BAND 

$CHAR L*1 0 256X5 SHADEPRINT CHARACTER SETS 

CUMPCT F*4 I 256X5 CUMULATIVE HISTOGRAM 


PDL (PROGRAM OUTLINE) 

DO FOR EACH SPECTRAL EAND. 

DO FROM 1=1 TO 256 BY 1 

INDEX = MINIMUM ( 16 , 1 + 16. * CUMPCT(I) / 100. 

DO FOR EACH STRIKE ( S ). 

$CHAR(I,S) = PRINT (INDEX, S) 

ENDDO. 

I? INDEX DOES NOT EQUAL I THEN 

PRINT K-1 , PPINT(INDEX,1) , PRINT (INDEX , 2) , 

ENDIF. 

ENDDO. 

ENDDO. 

jt****# END OP SPECIFICATION 


FOR ALL BANDS OF 


UNITS 


PERCENT 
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****** - SAP/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME MCPLGP 

DESIGNER STEPHEN W. MURPHRE Y 

FUNCTION THIS MODULE COMPUTES THE GRID-POINT LOCATIONS AND OTHER 

CONSTANTS FOR THE RESAMPLING PORTION OF THE MANUAL CONTROL- 
POINT LOCATION PROGRAM. CORRECTIONS FOR HIGH-FREQUENCY. 
HORIZONTAL MSS GEOMETRIC ERRORS ARE INCORPORATED INTO THE 
GRID-POINT LOCATIONS. 

INPUTS 1, SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON SLOCKS . . 


RESTRICTIONS ... 


CALLING SEQ. ... MCPLGP ( ) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

LENGTH 

1*4 

I 

— 

LENGTH OF INPUT-SPACE SUBIMAGE 

PIXELS 

WIDTH 

1*4 

I 

— 

WIDTH OF INPUT-SPACE SUBIMAGE 

PIXELS 

LFTSAM 

1*4 

I 

— 

COORDINATE OF LEFT SAMPLE 

PIXELS 

HSFAC 

1*4 

I 

— 

HORIZONTAL SCALE FACTOR 

PIXELS 

VSFAC 

1*4 

I 

— 

VERTICAL SCALE FACTOR 

PIXELS 

OLEN 

T*4 

0 

— 

LENGTH OF OUTPUT SPACE 

PIXELS 

OWID 

1*4 

0 

— 

WIDTH OF OUTPUT SPACE 

PIXELS 

T 

T* 4 

0 

— 

NUMBER OF GRID-PCINT SEGMENTS 

PIXELS 

N 

1*4 

0 

I 

# HYB.-SP. POINTS IN A SEGMENT 

PIXELS 

DH 

R*4 

0 

I 

INPUT-SPACE DISTANCE BETWEEN 
ADJACENT PIXELS IN A SEGMENT 

PIXELS 

H 

R*4 

0 

1X26 

GRID-POINT COORDINATES 

PIXELS 

D 

1*4 

0 

7SFAC 

VERTICAL RESAMPLING DISPLACEMENTS 

PIXELS 


PDL (PROGRAM OUTLINE) 


DO FROM M = 1 TO VSFAC BY 1 
D (M) = ( M - 1 ) / VSFAC 

ENDDO 

CLEN = VSFAC * ( LENGTH - 3 ) 

OWID = HSFAC * WIDTH - 3 ) 

ISWID = G( 50 / HSFAC ) 

CSWID = ISWID * HSFAC 
1=1 
C = 1 ft 

DO WHILE C LESS THAN OR EQUAL TO WIDTH - 2 
DO FOR EACH DETECTOR ( J = 1 TO 26 BY 1 ) 

IF HIGH-FREQUENCY CORRECTION WAS REQUESTED THEN 

COMPUTE F{ J , C + LFTSAM ) = HIGH-FREQUENCY CORRECTED COORD. 
H(I,J) = F( J , C+LFTSAM ) - LFTSAM 

EL SE 


H (I, J) 
ENDIF 


C 


ENDDO 

NIPIX = MINIMUM ( ISWID r WIDTH - 3 - C ) 
N (I) = NIPIX * HSFAC 
DH (I) = NIPIX / N (I) 

1 = 1+1 


C = DB (M) + NIPIX 
ENDDO 
1 = 1-1 


****** 


END OF SPECIFICATION 


*^**** 
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****** S AR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME MCPLHF 

DESIGNER STEPHEN H. MURPHREY 

FUNCTION THIS MODULE CORRECTS HORIZONTAL COORDINATES FOP HIGH-FR 

GEOMETRIC ERRORS. IT CONVERTS A HIGH-FREQUENCY-CORRECT 
COORDINATE TO AN UNCORRECTED-INPUT-SP ACE COORDINATE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. ... MCPLHF (HFC,LLF ,DET ,UIS) 


VARIABLE 

TYPE 

I/O DIM 

DESCRIPTION 

UNITS 

HFC 

R*4 

T 

HIGH -FREQUENCY -CORRECTED COORD. 

PIXELS 

ILF 

R*4 

I 

LINE-LENGTH FACTOR 


DET 

1*4 

I 

DETECTOR NUMBER 


UIS 

R*4 

0 

UNCORRECTED-INPUT-SPACE COORD. 

PIXELS 


PDL (PROGRAM OUTLINE) 

IF THIS IS FIRST EXECUTION OF THIS MODULE THEN 
COMPUTE COEFFICIENTS OF CORRECTION POLYNOMIALS. 
ENDIF 

EVALUATE CORRECTION POLYNOMIAL FOR THIS DETECTOR. 
****** END OF SPECIFICATION 


****** 


EQUENCY 

ED 


** **** 
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i 

NAME 

DESIGNER .... 
FUNCTION .... 

INPUTS 

OUTPUTS 

COMMON BLOCKS 
RESTRICTIONS 
CALLING SEQ. 


SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 
... MCPLVR 

... STEPHEN W. MURPHREY 

... THIS MODULE PERFORMS VERTICAL RESAMPLING AS PART OF THE 
MANUAL CONTROL-POINT LOCATION PROGRAM. 

... 1. SEE CALLING SEQUENCE. 

... 1. SEE CALLING SEQUENCE. 


. . . MCPLVR ( ) 


aK* 


PDL (PROGRAM OUTLINE) 

ADD FI P S T FOUR LINES TO THE HYBRID-SPACE BUFFER. 

INITIALIZE THE FOUR POINTERS TO THE HYBRID-SPACE BUFFER. 

DO FOR EACH HYBRID-SPACE LINE. 

DO FROM I = 1 TO K BY 1 f K * VERTICAL-SCALE -INCREASE FACTOR ) 
DO FOR EACH SPECTRAL BAND. 

DO FOR EACH OUTPUT SPACE PIXEL ( IN ONE LINE ) 

PERFORM CUBIC CONVOLUTION USING D (I) 

INCREMENT HYBRID-SPACE POINTERS BY 1 
ENDDO 

INCREMENT HYBRID-SPACE POINTERS TO NEXT SPECTRAL BAND. 
ENDDO 

WRITE ONE LINE OF OUTPUT DATA. 

ENDDO 

ADD NEXT HYBRID-SPACE LINE TO HYBRID-SPACE BUFFER. 

RESET THE FOUP POINTERS TO THE HYBRID-SPACE BUFFER. 

ENDDO 


Jit*#*## 


END OF SPECIFICATION 
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****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME MCPLHR 

DESIGNER STEPHEN 8. MDRPHPE Y 

FUNCTION THIS MODULE PERFORMS HORIZONTAL RESAMPLING AS PART OF THE 

MANUAL CONTROL-POINT LOCATION PROGRAM. GRID POINTS HAVE BEEN 
CORRECTED FOR HIGH-FREQUENCY, HORIZONTAL GEOMETRIC ERRORS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. 

• • • 

MCPLHR ( ) 


♦ 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

T 

-L. 

1*4 

I 

- 

NUMBER OF GEID-PCINT 

SEGMENTS 

N 

T* 4 

I 

I 

# HYB. -SP. POINTS IN 

A SEGMENT 

DH 

R*4 

I 

I 

INPUT-SPACE DISTANCE 

BETWEEN 





ADJACENT PIXELS IN A 

SEGMENT 

H 

R*4 

T 

1X26 

GRID-POINT COORDINATES 


PDL (PPOGRAM OUTLINE) 

DO FOR EACH SPECTRAL BAND. 

INITIALIZE POINTER TO JTH SET OF GRID-POINT COORDINATES. 

DO FOP EACH GRID-POINT SEGMENT ( M = 1 TC I BY 1 ) . 

G( * ) = GREATEST INTEGER NOT LARGER THAN * 

D = H (M, J) - G ( H (M , J) ) 

H = HjM.jj 

P = EOISTER TC INPUT-SPACE-BUFFER PIXEL G( H(M,J) ) - 1 
DO FOR EACH HYBRID-SPACE PIXEL IN CURRENT GRID-POINT SEGMENT. 

PERFORM CUBIC CONVOLUTION USING D. 

'E = D + EH (M) 

DO WHILE D >= 1.0 
D = D - 1.0 
P = P + 1 
FNDDO 
ENDDO 
' ENDDC 
ENDDO 

****** END OF SPECIFICATION ****** 
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****** SAB/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME MCPLPS 

DESIGNER STEPHEN fl. MDRPHREY 

FUNCTION PRINTS A SHADEPRINT OF A SINGIE SUBIMAGE THAT RESIDES IN THE 

(VIRTUAL) MEMORY OF THE COMPUTER. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. A PRINTER LISTING DATA SET. 

CALLING SEQ. ... CALL MCPLPS (BAND ,NSTRKS, OFFSET r SUBINFAD, CHARS AD) 

VARIABLE TYPE I/O DIM DESCRIPTIION UNITS 

BAND 1*4 X - SPECTRAL BAND NUMBER 

NSTRKS 1*4 I - NUMBER OF STRIKES TO BE PRINTED 

OFFSET 1*4 I - OFFSET INTO SUBIMAGE ARRAYS BYTES 

SUBINFAD 1*4 I - POINTER TO SUBIMAGE ARRAYS 


i-i -l. ix -is -t- —i i v a li u nu jj nnuti a 

CRAPSAD 1*4 I - POINTER TO CHARACTER-SET 

PDL (PROGRAM OUTLINE) 

INITIALIZE AND SAVE POINTER TO THE SHADEPRINT CHARACTER SET. 

INITIALIZE AND SAVE POINTER TO THE IMAGE DATA. 

INITIALIZE AND SAVE POINTER TO THE TITLE. 

INITIALIZE LEFTS AMPLE AND RIGHTS AMPLE- 1 POINTERS 
DO FROM LEFTSAMPLE TO HIGH TSA MPLE-1 BY 120 

PRINT HEADINGS AND SAMPLE NUMBERS ACROSS TOP OF LISTING. 

INITIALIZE LINE COUNTER 
INITIALIZE "FROM" POINTER 

INITIALIZE THE POINTER TO THE SHADEPRINT CHARACTER SET. 

BLANK OUT THE PRINT BUFFER 
DO FOP EACH LINE IN THE SUBIMAGE. 

PLACE THE LINE NUMBER IN THE LEET AND RIGHT EDGES 0? THE PRINT BUFFER. 
SET THE CARRIAGE-CONTROL CHARACTER TO A BLANK. 

DO FOP EACH STRIKE. 

MOVE AND TRANSLATE THE IMAGE DATA INTO THE PRINT BUFFER. 

PRINT THE DATA IN THE PRINT BUFFER. 

SET THE CARRIAGE-CONTROL CHARACTER TO A + . 

INCREMENT THE POINTER TO THE SHADEPRINT CHARACTER SET. 

ENDDO 

PRINT SAMPLE NUMBERS ACROSS THE BOTTOM OF THE LISTING. 

ENDDO. 

****** END OF SPECIFICATION ****** 
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B. 6 GEOMETRIC TRANSFORMATION PROGRAM 


B.6.1 Statement of Problem 

The objective of this program is to determine the mathematical transformation 
between two MSS image spaces. In particular, this program must produce data 
that controls the operation of the Resampling Program. 

The transformation is always between an input-image space and a resampled-image 
space. The input space can represent either an uncorrected or a previously 
resampled, image data set. There are four cases to be considered: 

a. Scene processing (with GCPs) of uncorrected or partially processed 
MSS data. 

b. Systematic processing (no GCPs) of uncorrected or partially processed 
MSS data. 

c. Completion processing of partially processed MSS data. 

d. Scale/orientation change of fully processed data. 

In the case of scene processing, a geometric error assessment from the geometric 
error models is also computed. In each case, the geometric transformation is 
specified in tabular form. It consists of input-space grid points (used for 
horizontal resampling) and hybrid-space grid points used for vertical resampling. 

If the input data is uncorrected or partially processed MSS data, there are three 
image spaces involved in determining the geometric transformation: 

a. Input space — MSS data that has not been geometrically corrected or 
resampled in any way. 

b. High- frequency-corrected input space — MSS data that has been hori- 
zontally resampled to remove all known high-frequency horizontal geo- 
metric errors. 

c. Corrected output space — MSS data (at a user-specified scale and 
orientation) that has been resampled to remove all known geometric 
errors . 

The transformations between these image spaces are illustrated in Figure B-6. 

The geometric transformation to be determined by this program is a tabular 
representation of the function H of Figure ; B-6. 
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Input Space High-Frequency-Corrected Corrected Output Space 

Input Space 


Figure B-6, Image Spaces 


If the input data is fully processed MSS data, there are only two image spaces 
involved in determining the geometric transformation: 

a. Fully processed input space — MSS data that has been resampled to 
remove all known geometric errors. 

b. Corrected output space — MSS data at a user-specified scale and 
orientation. 

The transformation between these image spaces is a tabular representation of a 
straightforward rotation and scale change. 


B.6.2 Data Flow 

The only data used by the Geometric Transformation program is ancillary m 
nature. This data comes from either the Ancillary data set or the User-Request 
Data Set. The resulting transformation is written into the Ancillary Data Set 
for use by the Reformatting program. 


B.6.3 Program Description 

The design of the Geometric Transformation Program is described in Figure B-7 
and the program specifications that follow. 
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GTRN 


GTHANC 

GTEURQ 


GTEFMD 


GTOSDF 


GTGDPT 


GTEBAS 

GTWANC 


I IOTP 

MAPLL 

LLEC 

ECCTP 

DELTAS 

AAEIT 

AAFIT2 ADTAA 

CTPEC 

ECLL 

ILMAP 

MAPOI 

INITOS 

OIMAP 

HAP1L 

LLEC 

ECCTP 

CTPUTP 

OTPII 

IIOTP 

OTPCTP 


FIGURE B-7. MODULE HIERARCHY FOR GEOMETRIC TRANSFORMATION PROGRAM. 
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****** SAB/LMDS AT DATA MERGING SYSTEM PROGRAM SPECIFIC A TION ****** 

'NAME GTEN 

DESIGNER STEPHEN W. HURPHREY 

FUNCTION THIS IS THE TOP LEVEL DRIVER FOR THE LANDSAT MSS GEOMETRIC 

TRANSFORMATION PROGRAM. 

INPUTS 1. AN SLDHS ANCILLARY DATA SET FOR THE SCENE TO BE PROCESSED. 

2. A DATA SET CONTAINING USER*S REQUESTS (ON CARD IMAGES) . 

OUTPUTS 1. THE UPDATED SLDMS ANCILLARY DATA SET. 

2. A PRINTER LISTING RATA SET. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... NOT APPLICABLE. 

PDL (EEOGFAM OUTLINE) 

PRINT SEPARATOR PAGE TO INDICATE START OF RUN. 

READ USER'S REQUESTS. 

FEAD ANCILLARY DATA SET TO INITIALIZE DATA AREAS. 

CAS ENTRY 

CASE 1: SCENE PROCESSING OR SYSTEM PROCESSING. 

PERFORM ERROR MODELING (ATTITUDE/ALTITUDE FIT) . 

DEFINE THE OUTPUT SPACE. 

CREATE INTERPOLATION GRID POINTS ( FROM ERROR MODELS ) . 

PERFORM ERROR ASSESSMENT ( FROM EPRCP MODELS ) , 

WRITE GRID POINT CORRESPONDENCE TO ANCILLARY DATA SET. 

CASE 2: COMPLETION PROCESSING. 

DEFINE THE OUTPUT SPACE. 

CREATE INTERPOLATION GRID POINTS ( FROM INPUT HRS/VRS COORDINATES ). 
PERFORM ERROR ASSESSMENT ( FROM INPUT ANCILLARY DATA ) . 

WFITE GRID POINT CORRESPONDENCE TO ANCILLARY DATA SET. 

CASE 3: SCALI/ORIENTATION CHANGE. 

DEFINE THE OUTPUT SPACE. 

CREATE INTERPOLATION GRID POINTS ( BASED ON SCALE AND ORIENTATION ) . 
PERFORM ERROR ASSESSMENT ( FROM INPUT ANCILLARY DATA ). 

WRITE GRID POINT CORRESPONDENCE TO ANCILLARY DATA SET. 

ENLC AS E 

PRINT SEPARATOR PAGE TO INDICATE END OF RUN. 

****** END OF SPECIFICATION ****** 
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FUNCTION 
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SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 
... GTRANC 

... STEPHEN W . MURPHREY 

... READS THE SLDMS ANCILLARY DATA SET 


INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. MSSCCNS 
RESTRICTIONS ... 

CALLING SEQ. CALL GTRANC () 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

PRINT PROGRAM-ENTPY MESSAGE. 

DO UNTIL END-OF-DATA ON ANCILLARY DATA SET. 

READ NEXT RECORD (ID RECORD) . 

IF IT IS "CCT-HEADER DATA", THEN 

READ NEXT TWO RECORDS (DATA RECORDS). 

PRINT CCT- HEADER DATA. 

ENDIF 

IF IT IS "SIAT DATA " , THEN 

READ NEXT TWENTY-CNE RECORDS (DATA RECORDS) 
PRINT SIAT DATA. 

ENDIF 
END DO 

PRINT PROGRAM-EXIT MESSAGE. 


****** 


END OF SPECIFICATION 


****** 


****** 
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****** 


****** 


SAR /L AND SA 1 DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME GTEUBQ 

DESIGNEB STEPHEN W. HUKPHREY 

FUNCTION HEADS USEE’S REQUESTS EOS "LMSSGT" PROGRAM. 


INPUTS 

OUTPUTS 


1 . 

1 . 

2 . 


SEE CALLING SEQUENCE. 

SEE CALLING SEQUENCE. 

A PRINTER LISTING DATA SET. 


COMMON BLOCKS . . 

RESTRICTIONS ... 

CALLING SEQ. ... CALL GTRURQ () 

VARIABLE TYPE I/O DIM DESCRIPTION 


UNITS 


PDX (PROGRAM OUTLINE) 

PPINT PROGRAM-ENTRY MESSAGE. 

INITIALIZE NAMELIST PARAMETERS. 

READ NAMELIST OF USEE'S REQUESTS AND PARAMETERS. 

T>T?CTJT*rS 

PRINT USER'S REQUESTS AND DEFINING PARAMETERS. 
PPINT PROGRAM-EXIT MESSAGE. 


**** ** 


END OF SPECIFICATION 


****** 
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S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 


NAME GTERMD 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION IT COMPUTES THE 14 COEFFICIENTS WHICH MODEL THE ATTITUDE AND 

ALTITUDE OVER A SCENE. THE COMPUTATION IS DONE BY A WEIGHTED 


LEAST-SQUARES FIT TO CONTROL-POINT AND AMS/EPREMERIS DATA. AS 
A BY-PRODUCT OF THE FIT, COVARIANCE MATRICES WHICH ARE LATER 
USED TO ESTIMATE THE ACCURACY OF THE COMPUTED COEFFICIENTS, ARE 
COMPUTED. 


INPUTS 1. SEE CALLING SEOUENCE. 

OUTPUTS ........ 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS . 
RESTRICTIONS . . 


CALLING SEQ. 


CALL GTERMD (NMCP, NATEST, NALEST, KWALTY, AMSFOL, 
AMSPIT, AMSYAW, AMSALT, ATTTIM, ALTTIM, TIMECP, 
CTPH , DELTAV, DELTAH, ALTNOM, AAMOD, PINV, QINV) 


VARIABLE 

TYPE 

I/O 

DIM 

NMCP 

1*2 

I 


NATEST 

1*2 

I 

— 

NALEST 

1*2 

I 

- 

KWALTY 

T*2 

T 

NMCP 

AMSROLL 

R*8 

I 

NATEST 

AMSPIT 

B*8 

I 

NATEST 

AMSYAW 

R*8 

I 

NATEST 

AMSALT 

R* 8 

I 

NALE ST 

ATTTIM 

R*8 

I 

NATEST 

ALTTIM 

R*8 

I 

NALEST 

TIMECP 

R*8 

T 

NMCP 

CTPH 

R*8 

I 

NMCP 

DELTAV 

F*8 

I 

NMCP 

DELTAH 

F *8 

T 

NMCP 

ALTNOM 

P*8 - 

I 

— 

AAMOD 

R*8 

0 

14 

PINV 

R*8 

0 

6X6 

QINV 

R*8 

0 

8X8 


DE SCRIPTIION UNITS 

# OF CONTROL POINTS 

# OF ATTITUDE ESTIMATES 

# OF ALTITUDE ESTIMATES 
STATUS OF CONTROL POINTS 

AMS ROLL VALUES RADIANS 

AMS PITCH VALUES RADIANS 

AMS YAW VALUES RADIANS 

AMS ALTITUDE VALUES METERS 

TIME ASSOCIATED WITH ROLL, SECONDS 

PITCH AND YAW VALUES 

TIME ASSOCIATED WITH ALTITUDE SECONDS 
TIME ASSOCIATED WITH SECONDS 

CONTROL POINTS 

HORIZONTAL COORDINATES METERS 

OF CONTROL POINTS IN CTP 

VERTICAL RESIDUAL OF METERS 

CONTROL POINTS IN CTP 

HORIZONTAL RESIDUALS OF METEPS 

CONTROL POINTS 

NOMINAL ALTITUDE OF SPACECRAFT METERS 
COEFFICIENTS OF ROLL, PITCH, 

YAW, AND ALTITUDE 
COVARIANCE MATRIX {ROLL AND 
ALTITUDE) 

COVARIANCE MATRIX (PITCH AND YAW) 


PDL (PROGRAM OUTLINE) 


MAP OBSERVED INPUT-SPACE CP COORDINATES TO UNCORRECTED TANGENT -S PACE. 

MAP NOMINAL MAP CP COORDINATES TO CORRECTED TANGENT SPACE. 

COMPUTE DIFFERENCES OF UNCORRECTED AND CORRECTED TANGENT-SPACE COORDINATES, 

IF THERE APE TWO OR MORE CONTROL POINTS THEN 
INITIALIZE I SING FLAG TO BE NONZERO. 

ATTEMPT TO PERFORM ATTITUDE/ ALTITUDE FIT WITH "AAFIT" ALGORITHM. 

IF "AAFIT" IS UNSUCCESSFUL THEN 
SET ISING FLAG TO ZERO. 

FHDIF 

ELSE 

INITIALIZE ISING FLAG TO BE ZERO. 

ENDIF 

IF ISING FLAG IS ZERO THEN 

PERFORM ATTITUDE/ALTITUDE FIT WITH " AAFIT2" ALGORITHM. 

ENDIF 

***** END OF SPECIFICATION ****** 
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****** SAR/L AND SAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

JNAME IIUTP 

DESIGNEE STEPHEN W. MDBPHREY 

FUNCTION CONVERTS AN AEBAY OF INPUT-SPACE COOPDINATES (LINE AND SAMPLE) 

TO UNCOBREC TED-TANG ENT -SP ACE COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. MSSCONS 
RESTRICTIONS ... 


CALLING SEO. 


CALL IIUTP (NMCP, VRA, VRC, VRAN, 
CLINE, CSAMP, KWALTY, UTPV, UTPH, 


VPCN , DELVV 
DTI ME, CTIM 


k 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

NMCP 

1*2 

I 

_ 

# OF POINTS 

— 

VRA 

R*8 

T 

NMCP 

EARTH ROTATION VELOCITY 
(ALONG -TRACK) AT INPUT POINT 

METERS/ 

SEC 

VRC 

R*8 

I 

NMCP 

EARTH ROTATION VELOCITY 
(ACROSS-TRACK) AT INPUT POINT 
EARTH ROTATION VELOCITY 
(ALONG-TRACK) AT NADIR 

n 

VRAN 

F*8 

I 


«« 

VRCN 

R* 8 

I 

” 

EARTH ROTATION VELOCITY 
(ACROSS-TRACK) AT NADIR 

it 

DELVV 

R* 8 

T 

— 

SPACECRAFT VELOCITY ERROR 

it 

CLINE 

■R*8 

I 

NMCP 

LINE OF CONTROL POINT 

PIXELS 

CSAMP 

R*8 

I 

NMCP 

SAMPLE OF CONTROL POINT 

PIXELS 

KWALTY 

1*2 

I 

NMCP 

STATUS OF POINTS 

PIXELS 

UTPV 

.R*8 

0 

NMCP 

VERTICAL COORDINATE OF POINT 
IN UTP 

METERS 

UTPH 

?*8 

0 

NMCP 

HORIZONTAL COORDINATE OF 
POINT IN UTP 

METERS 

CTIME 

P *8 

0 

NMCP 

TIME (CONTINUOUS) ASSOCIATED 
WITH POINT 

SECONDS 

DTIME 

R*8 

0 

NMCP 

TIME (DISCRETE) ASSOCIATED 

WITH POINT 

SECONDS 


PDL (PFOGRAM OUTLINE) 

PRINT SUBROUTINE-ENTRY MESSAGE. 

PRINT HEADINGS. 

DO FOR EACH CONTROL POINT: 

CONVERT TO IMAGS-CENTER-OFIGI N COORDINATES. 

COFRECT HORIZONTAL COORDINATE FOR MIRROR-VELOCITY ERRORS. 

CORRECT FOR DIFFERENTIAL INPUT SCALES. 

CORRECT VERTICAL COORDINATE FOR SCAN-SKEW ERRORS. 

CORRECT VERTICAL COORDINATE FOR SPACECRAFT-VELOCITY ERRORS. 

CORRECT FOR DISTORTIONS DUE TO EARTH ROTATION. 

CORPECT HORIZ. COORDINATE FOR EA PTH -CURVATURE / PANORAMIC-PROJECTION ERRORS. 
PRINT RESULTS OF ABOVE ERROR CORRECTIONS. 

ENDDO 

PRINT SUBFOUTINE-EXIT MESSAGE. 

****** END OF SPECIFICATION ****** 
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****** SAF/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME DELTAS 

DESIGNER STEPHEN W. MORPHREY 

FUNCTION COMPUTE THE DIFFERENCES OF THE UNCORRECTED-TA NGENT-SPACE 

COORDINATES AND THE CORRECTED-T ANGENT-SPACE COORDINATES OF 
THE CONTROL POINTS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 


RESTRICTIONS ... 

CALLING SEQ. ... CALL DELTAS (NMCP, KWALTY, CTPV, CTPH, UTPV, UTPH, 

DELTAV, DELTAH) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTIION 

UNITS 

NMCP 

1*2 

T 

— 

# OF CONTROL POINTS EQUAL 

— 

KWALTY 

1*8 

I 

NMCP 

STATUS OF CONTROL POINTS 


CTPV 

F*8 

I 

NMCP 

CTP VERTICAL COORDINATES 

METERS 

CTPH 

R*8 

T 

NMCP 

CTP HORIZONTAL COORDINATES 

METERS 

UTPV 

R* 8 

T 

NMCP 

UTP VERTICAL COORDINATES 

METERS 

UTPH 

R*8 

I 

NMCP 

UTP HORIZONTAL COORDINATES 

METERS 

DELTAV 

R* 8 

0 

NMCP 

DIFFERENCES BETWEEN VERTICAL 
UTP AND CTP COORDINATES 

METERS 

DELTAH 

R*8 

0 

NMCP 

DIFFERENCES BETWEEN HORIZONTAL 
UTP AND CTP COORDINATES 

METERS 


fc>DL (PROGRAM OUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE. 

PRINT HEADINGS. 

DC FOR EACH CONTROL POINT 

IF THIS IS A GOOD CONTROL POINT (AS DETERMINED BY THE KWALTY ARRAY) 

SET DELTAV = UTPV-CTPV FOR THIS POINT 
SET DELTAH = UTPH-CTPH FOP THIS POINT 
INDIF 

PRINT RESULTS. 

ENDDO 

PPINT PROGRAM-EXIT MESSAGE. 

****** - END 0? SPECIFICATION ****** 
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****** SAP /LAND SAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME AAFIT 

DESIGNER STEPHEN N. MDRPHFEY 

FUNCTION COMPUTE THE 14 COEIEICIENTS WHICH MODEL THE ATTITUDE AND 

ALTITUDE OF EE A SCENE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. ... CALL AAFIT (NMCP, NATEST, NALEST, KWALTY, AMSROL, 

AMSPIT, AMSYAW, AMSALT, ATTTI M , ALTTIM, TIMECP, 
CTPH, DELTAV, DELTAH, ALTNOM, AAMOD, PINV, QINV) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPT I ION 

UNITS 

NMCP 

1*2 

I 


# OF CONTROL POINTS 

_ 

NATEST 

1*2 

I 

— 

# OF ATTITUDE ESTIMATES 

- 

NALEST 

1*2 

I 

— 

# OF ALTITUDE ESTIMATES 

— 

KWALTY 

1*2 

T 

NMCP 

STATUS OF CONTROL POINTS 

- 

AMSROLL 

R*8 

i 

NATEST 

AMS ROLL VALUES 

RADIANS 

AMSPIT 

P*8 

i 

NATEST 

AMS PITCH VALUES 

RADIANS 

AM SY AW 

R* 8 

T 

NATEST 

AMS YAW VALUES 

RADIANS 

AMSALT 

R*8 

I 

NALEST 

AMS ALTITUDE VALUES 

METERS 

ATTTIM 

R*8 

I 

NATEST 

TIME ASSOCIATED WITH ROLL, 
PITCH AND YAW VALUES 

SECONDS 

ALTTIM 

R*8 

T 

u. 

NALEST 

TIME ASSOCIATED WITH ALTITUDE 

SECONDS 

TIMECP 

R* 8 

I 

NMCP 

TIME ASSOCIATED WITH 
CONTROL POINTS 

SECONDS 

CTPH 

R*8 

T 

NMCP 

HORIZONTAL COORDINATES 
OF .CONTROL POINTS IN CTP 

METERS 

DELTAV 

R*8 

I 

NMCP 

VERTICAL RESIDUAL OF 
CONTRCL POINTS IN CTP 

METERS 

DE1TAH 

R*8 

I 

NMCP 

HORIZONTAL RESIDUALS OF 
CONTROL POINTS 

METERS 

ALTNOM 

R* 8 

I 

— 

NOMINAL ALTITUDE OF SPACECRAFT 

METERS 

AAMOD 

R* 8 

0 

14 

COEFFICIENTS OF ROLL, PITCH, 
YAW, AND ALTITUDE 

** 

PINV 

R*8 

0 

6X6 

COVARIANCE MATRIX (ROLL AND 
ALTITUDE) 

— 

QINV 

R* 8 

0 

8X8 

COVARIANCE MATRIX (PITCH AND YAW) 


PDL (PROGRAM OUTLINE) 


PRINT PROGRAM-ENTRY MESSAGE. 

PRINT HEADINGS. 

INITIALIZE "P n , "Q", "5", AND "S" MATRICES TO ZEROS. 

CALCULATE ATTITUDE DIEFERENCES. 

CALCULATE ALTITUDE DIEFERENCES. 

CALCULATE TIME CONSTANTS. 

COMPUTE ELEMENTS OF THE "R» MATRIX. 

COMPUTE ELEMENTS OF THE "S" MATRIX. 

COMPUTE ELEMENTS OF THE "P" MATRIX. 

COMPUTE ELEMENTS OF THE »Q" MATRIX. 

PRINT THE "P", "Q", "S'*, AND «S" MATRICES. 

SOLVE THE TWO SYSTEMS OF EQUATIONS. 

PRINT FROGBAM-EXIT MESSAGE.' 

****** END OF SPECIFICATION ****** 
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SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME AAFIT2 

DESIGNER STEPHEN W. MDRPHRE Y 

FUNCTION COMPUTE THE 14 COEFFICIENTS WHICH MODEL THE ATTITUDE AND 

ALTITUDE OVER A SCENE WHEN "AAFIT" IS NOT APPLICABLE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. 


CALL AAFIT2 {NMCP, NATEST. NALEST, KWALTY, AMSROL, 
AMSPIT, AMSYAW, AMSALT, ATTTIM, ALTTI M, TIMECP, 
CTPH, DELTA?, DELTAH, ALTNOM, AAMOD, PINV, QINV) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPT I ION 

UNITS 

NMCP 

1*2 

T 

— 

# OF CONTROL POINTS 


NATEST 

1*2 

I 

— 

# OF ATTITUDE ESTIMATES 

- 

NALEST 

T*2 

I 

— 

# OF ALTITUDE ESTIMATES 

- 

KWALTY 

1*2 

I 

NMCP 

STATUS CF CONTROL POINTS 

— 

AMSPOLL 

R*8 

T 

NATEST 

AMS ROLL VALUES 

RADIANS 

AMSPIT 

F* 8 

I 

NATEST 

AMS PITCH VALUES 

RADIANS 

AMSYAW 

R* 8 

I 

NATEST 

AMS YAW VALUES 

RADIANS 

AMSALT 

P.*8 

I 

NALEST 

AMS ALTITUDE VALUES 

METERS 

ATTTIM 

R*8 

I 

NATEST 

TIME ASSOCIATED WITH ROLL, 
PITCH AND YAW VALUES 

SECONDS 

ALTTIM 

P*8 

T 

NALEST 

TIME ASSOCIATED WITH ALTITUDE 

SECONDS 

TIMECP 

R* 8 

I 

NMCP 

TIME ASSOCIATED WITH 
CONTROL POINTS 

SECONDS 

CTPH 

p*8 

I 

NMCP 

HORIZONTAL COORDINATES 
OF CONTROL POINTS IN CTP 

METERS 

DELTAV 

R* 8 

I 

NMCP 

VERTICAL RESIDUAL OF 
CONTROL POINTS IN- CTP 

METERS 

DELTAH 

R*8 

I 

NMCP 

HORIZONTAL RESIDUALS OF 
CONTROL POINTS 

METERS 

ALTNOM 

R*8 

T 

— 

NOMINAL ALTITUDE OF SPACECRAFT 

METERS 

AAMOD 

P* 8 

6 

14 

COEFFICIENTS OF ROLL, PITCH, 
YAW, AND ALTITUDE 

"" 

PINV 

R*8 

0 

6X6 

COVARIANCE MATRIX {ROLL AND 
ALTITUDE) 

** 

QINV 

R*8 

0 

8X8 

COVARIANCE MATRIX (PITCH AND YAW) 


PDL (EROGRAM OUTLINE) 


PRINT PROGP AM- ENTRY MESSAGE. 

PFINT HEADINGS. 

COMPUTE AMS ATTITUDE AND ALTITUDE MODEL COEFFICIENTS. 
ADJUST GCP DATA FOR AMS RATE AND ALTITUDE. 

COMPUTE ROLL BIAS TERM. 

COMPUTE PITCH AND YAW BIAS TERMS. 

COMPUTE NEW PREVIOUS EIASES AND VARIANCES. 

WRITE QUALITY DATA. 

PRINT PROGRAM-EXIT MESSAGE. 




END OF SPECIFICATION 


. -r 




****** 


****** 
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k ****** S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

IsAME ADTAA 

DESIGNER STEPHEN H. MURPHREY 

FUNCTION TO COMPUTE THE 14 ATTITUDE / ALTITUDE COEFFICIENTS BASED 

AMS ESTIMATES OF ROLL, PITCH, YAH, AND ALTITUDE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1 . SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 

RESTRICTIONS ... THE NUMBER OF AMS ATTITUDE VALUES AND THE NUMBEP OF 

EPHEMEPIS ALTITUDE VALUES MUST BOTH BE GREATER THAN 1. 


CALLING SEQ. ... CALL ADTAA (NATEST, NALEST. AHSROL. AMSPIT, 

AMSYAW, AMSALT, ATTTIM , ALTTIM, ALTNOM, AAAMOD) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

NATEST 

1*2 

T 

— 

# OF ATTITUDE ESTIMATES 

— 

NALEST 

1*2 

I 

— 

# OF ALTITUDE ESTIMATES 

— 

AMSPOL 

R*8 

I 

NATEST 

AMS ROIL VALUES 

RADIANS 

AMSPIT 

R*8 

I 

NATEST 

AMS PITCH VALUES 

RADIANS 

AMSYAW 

P*8 

I 

NATE ST 

AMS YAW VALUES 

RADIANS 

AMSALT 

R*8 

T 

NALEST 

AMS ALTITUDE VALUES 

METERS 

ATTTIM 

R*8 

I 

NATEST 

TIME ASSOCIATED WITH ROLL, 
PITCH, YAW VALUES 

SECONDS 
FROM IC 

ALTTIM 

R* 8 

I 

NALEST 

TIME ASSOCIATED WITH 
ALTITUDE VALUES 

SECONDS 
FROM IC 

AITNCM 

R*8 

I 

— 

NOMINAL ALTITUDE 

METERS 

AAAMOD 

R* 8 
( 

0 

14 

COEFFICIENTS OF POLL. PITCH 
YAW AND ALTITUDE DEVIATION 



PDL (PROGPAM OUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE. 

INITIALIZE AP, AQ, RM, PM, YM, AND HM MATRICES TO ZEROS. 
CALCULATE AP, AQ , RM, PM, YM, AND HM ELEMENTS. 
PRECONDITION AP MATRIX. 

INVERT AND TEST AP MATRIX. 

COMPUTE ROIL, PITCH, AND YAW MODEL COEFFICIENTS. 
PRECONDITION AQ MATRIX. 

INVERT AND TEST AQ MATRIX. 

COMPUTE ALTITUDE MODEL COEFFICIENTS. 

PRINT RESULTS. 

PRINT PROGRAM-EXIT MESSAGE. 

****** END OF SPECIFICATION 


steifc 


ON 
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****** S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME GTOSDF 

DESIGNER STEPHEN W. MURPHBE Y 

FUNCTION THIS SUBROUTINE DEFINES THE OUTPUT SPACE. IT COMPUTES THE 

OUTPUT-SPACE PIXEL COORDINATES OF THE NOMINAL FORMAT CENTER. 

INPUTS 1. SEE CALLING SEQUENCE.- 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL GTOSDF { ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE. 

COMPUTE T?NGENT- SPACE COORDINATES OF FORMAT CENTER. 

TRANSFORM FORMAT-CENTER COORDINATES TO GEODETIC LATITUDE AND LONGITUDE 
TPANSEORM FORMAT-CENTER COORDINATES TO MAP COORDINATES. 

INPUT NOMINAL-FORMAT -CENTER MAP COORDINATES. 

COMPUTE OUTPUT-SPACE COORDINATES OF NOMINAL FORMAT CENTER. 

SCALE AND ROUND THE NCMINAL-FORM AT -CENTER COORDINATES. 

PRINT RESULTS. 

PRINT PROGRAM-EXIT MESSAGE. 


***:£** 


END OF SPECIFICATION 


****** 
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****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME CTPEC 

DESIGNER STEPHEN W. MURPHREY 


FUNCTION TC MAE EOINTS FROM THE CORRECTED TANGENT SPACE TO THE EARTH- 

CENTERED CARTESIAN COORDINATE SYSTEM. THE MAPPING IS, IN 
ORDER, A POTATION (FROM THE TANGENT-SPACE ORIENTATION TO THE 
EARTH-CENTERED ORIENTATION), A TRANSLATION (FROM THE TANGENT- 
SPACE ORIGIN TO THE EARTH -CENTERED ORIGIN) , AND A PROJECTION 
(FROM THE TANGENT SPACE ONTO THE EARTH SURFACE) . 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 

COMMON BLOCKS .. 

RESTRICTIONS ... 


CALLING SEQ. ... CALL 


CTPEC (XEX, XEY , XECNAD, YEX, YEY , YECNAD , ZEX, 
ZECNAD, A SQR , BSQR, A4TH , B4TH, V, NUM, 


Z.-EI, ZEUNAU, A by K , BbU'f, 1 

CTPV , CTPH, XEC, YEC, ZEC) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UIITS 

XEX 

R* 8 

T 

— 

ROTATION MATRIX ELEMENT 

— 

XEY 

p.*e . 

I 

- 

ROTATION MATRIX ELEMENT 

- 

XECNAD 

R* 8 

T 

— 

NADIR EARTH CENTERED X COORDINATE 

METERS 

YEX 

R*e 

I 

- 

ROTATION MATRIX ELEMENT 

- 

YEX 

R*8 

T 

— ' ' 

ROTATION MATRIX ELEMENT 

- 

YECNAD 

R*8 

I 

- 

NADIR EARTH CENTERED Y COORDINATE 

METERS 

ZEX 

R*8 

I 

- 

ROTATION MATRIX ELEMENT 

- 

ZEY 

R*8 

I 

— 

ROTATION MATRIX ELEMENT 

— 

ZECNAD 

B* 8 

I 

— 

NADIR EARTH CENTERED Z COORDIANTE 

METERS 

ASQR 

R* 8 

I 

— 

PRECOMPUTED EARTH A**2 

— 

BSQR 

R*8 

I 

- 

PRECOMPUTED EARTH B**2 

- 

A4TH 

p*8 

I 

- 

PRECOMPUTED EARTH A**4 

- 

B4TH 

E*8 

T 

— 

PRECOMPUTED EARTH B**4 

- 

V 

P*8 

I 

- 

PROJECTION PARAMETERS 

- 

NUM 

2*8 

T 

- 

NUMBER OF POINTS TC MAP 

— 

CTPV 

R* 8 

I 

NUM 

CTP VERTICAL COORDINATES 

METERS 

CTPH 

R* 8 

I 

NUM 

CTP HORIZONTAL CCORDIANTES 

METERS 

XEC 

B*8 

T 

NUM 

EARTH CENTERED X OF THE POINTS 

METERS 

XEC 

P*8 

T 

NUM 

EARTH CENTERED Y OF THE POINTS 

METERS 

XEC 

R*8 

I 

NUM 

EARTH CENTERED Z OF THE POINTS 

MET ERS 


PDL (EROGRAM OUTLINE) 


DC FOR EACH ICINT 

ROTATE AND TRANSLATE THE TANGENT SPACE COORDINATES OF THE POINT 
INTO EARTH CENTERED COORDINATES 
X ECT = CV*XEX+CH*XEY+ XECNAD 
YECT = CV*YEX+CH*YEY+ YECNAD 
Z ECT = CV*ZFX+CH*ZEY+ ZECNAD 

COMPUTE PARAMETERS USED TO PROJECT THE POINT FROM THE TANGENT 
SPACE TO THE EARTH SURFACE 

PI = BS0R*(V(1)*XECT+V(2)*YECI)+A4TH/BSQR*V(3)*ZECT 
P2 = ASQR*BSQE+V (3) *V (3) * (A4TH/B4TH-1.0D0) 

P3 = BSQR*7XECT*XECT+YECT*yECT) + ASQR* (ZECT*ZECT-BSQE) 

S = (P1-DSQET(P1*P1-P2*P3) )/P2 
PROJECT THE POINT FPCM THE TANGENT SPACE TO THE EARTH SURFACE 
XEC = XECT-V(1)*S 
YEC = YECT-V (2) *S 
ZEC = ZECT-V (3) *S*ASQR/BSQR 

ENDDO 


_ * ;$c 


END OF SPECIFICATION 


****** 
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****** 


SAB/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFIC A TT OB 


NAME ECLL 

DESIGNED STEPHEN W. MURPHREY 

FUNCTION MAP POINTS FROM EARTH-CENTERED COORDINATES TO GEODETIC 

LATITUDE AND LONGITUDE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. 

• • • 

CALL 

ECLL 

(AELLPS , BELLPS, NUM, XEC, YEC, ZEC, 

RD LA TR , 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTITON 

UNITS 

AELLPS 

R*8 

I 

— 

MAJOR AXIS OF EARTH ELLIPSOID 

METERS 

BELLPS 

R*8 

I 

- 

MINOR AXIS OF EARTH ELLIPSOID 

METERS 

NUM 

1*2 

I 

- 

NUMBER OF POINTS TO TRANSFORM 

- 

XEC 

R*8 

I 

NUM 

EARTH CENTERED X COORDINATE 

METERS 

YEC 

R*8 

T 

NUM 

EARTH CENTERED Y COORDINATE 

METERS 

ZEC 

E*8 

I 

NUM 

EARTH CENTERED Z COORDINATE 

METERS 

RDLATR 

5*8 

0 

NUM 

GEODETIC LATITUDE OF INPUT POINTS 

RADIANS 

RDLONR 

R*8 

0 

NUM 

LONGITUDE OF INPUT POINTS 

RADIANS 


PD1 (PROGRAM OUTLINE) 

DO FOR EACH POINT 

LATITUDE=ARCTAN ( (ASQP/BSQR) *ZEC/SQRT f XEC*XEC+ YEC*YEC) ) 
LONGITUDE=ARCSIN (DABS (Y EC) /SQRT (XEC*XEC+YEC*YEC) ) 

IE (XEC > OR = 0) 

IF (YEC < 0) 

LONGITUDE= -LONGITUDE 
ENDIF 
ELSE 

IF (YEC < U) 

LONGITUDE=LONGITUDE-PI 

ELSE 

LONGITUDE= PI -LONGITUDE 
ENDIF 
ENDIF 

PRINT RESULTS FOR THIS POINT 
ENDDO 

****** END OF SPECIFICATION 


SLONGR) 


Jcic £5$: 
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****** SAB /LAND SAT DATA HEDGING SYSTEM PROGRAM SPECIFICATION ****** 

^1AME XIMAP 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION MAP GEODETIC LATITUDE AND LONGITUDE COORDINATES TO (UTM) MAP 

COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS ........ 1. SEE CALLING SEQUENCE. 

2. A EEINTER LISTING DATA SET. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. 

• • « 

CALL 

LLMAP 

( AELLPS , BELLPS , NUM . IZONED , 
DNORTH , DEAST , IZONE } 

PDLATR , RLONGR 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRI PTIION 

UNITS 

AELLPS 

p*8 

I 

— 

MAJOR AXIS OF EARTH ELLIPSOID 

METERS 

BEILPS 

B*8 

I 

— 

MINOR AXIS OF EARTH ELLIPSOID 

METERS 

NUM 

1*4 

I 

— 

NUMBER OF POINTS TO TRAN SFOPM 


IZONED 

1*4 

I 


UTM ZONE DESIRED 

— 

RDLA.TR 

R*8 

T 

NTJM 

GEODETIC LATITUDE OF INPUT POINTS 

RADTANS 

RDLONR 

R* 8 

i 

NUM 

LONGITUDE OF INPUT POINTS 

RADIANS 

DNOP TH 

R*8 

0 

NUM 

UTM NORTHING COORDINATE 

METERS 

DEAST 

R*8 

0 

NUM 

UTM EASTING COORDINATE 

METERS 

IZONE 

1*4 

0 

NUM 

UTM ZONE 


COMMENTS ... 

• • • ♦ 

THIS 

MODULE 

IS A MODIFICATION OF A FORTRAN PROGRAM 


OETAINED FROM THE UNITED STATES GEOLOGICAL SURVEY, 
p****** END OF SPECIFICATION ****** 
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SAR /LAND SAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 


'NAME MAPOI 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION CONVERT (UTM) MAP COORDINATES TO OUTPUT-SPACE PIXEL 

COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 


RESTRICTIONS ... 


CALLING SEQ. 

* * * 

CALL 

MAPOI 

( NUM , T ZONED , CNORTH , CEAST 

r BETA 





DNORTH , DEAST , IZONE , HO 

, VO 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPT I ION 

UNITS 

NUM 

T*4 

I 

- 

NUMBER OF POINTS TO TRAN SFOPM 

— 

EETA 

R*8 

I 

- 

HEADING ANGIE 

RADIANS 

IZONED 

1*4 

I 


UTM ZONE DESIRED 


CNORTH 

p*8 

I 

NTJM 

UTM NORTHING OF FORMAT CENTER 

METERS 

CEAST 

B* 8 

I 

NUM 

UTM EASTING C? FORMAT CENTER 

METERS 

DNCRTH 

R* 8 

I 

NUM 

UTM NORTHING COORDINATE 

METERS 

DEAST 

P*8 

I 

NUM 

UTM EASTING COORDINATE 

METERS 

IZONE 

1*4 

I 

NUM 

UTM ZONE 


HO 

R* 8 

I 

NUM 

OUTPUT SPACE COORDINATE 

PIXELS 

VO 

R*8 

T 

NUM 

OUTPUT SPACE COORDINATE 

PIXELS 


k PDL (PPOGRAM CUTLINE) 

ALPHA =3 * PI / 2 - BETA 
SA = SINE ( ALPHA ) 

CA = COSINE { ALPHA ) 

DC FOR EACH PCINT 


Horn = 
VO (I) = 
ENDDO 

(SA * ( 
(CA * 1 

' DNOPTH ( 
1 DNORTH ( 

I) - CNORTH ) - CA * ( 
I) - CNORTH ) + SA * ( 

DEAST (I 
DEAST (I 

) - CEAST ) 

) - CEAST j 

1 * 

HSCALE 

VSCALE 

****** 



END OF SPECIFICATION 




****** 
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k *****'*' ' SAR/LANDSA1 DATA MERGING SYSTEM FROGRAM SPECIFICATION ****** 

'Name gtgdpt 

DESIGNER STEPBEN W. MURPHREY 

FUNCTION THIS SUBROUTINE CALCULATES THE GRID-POINT CORRESPONDENCE THAT 

IS THE GEOMETRIC TRANSFORMATION BETWEEN THE OUTPUT SPACE AND 
THE INPUT SPACE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. PRINTER LISTING. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. . . . CALL GTGDPT ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGPAM OUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE. 

INITIALIZE OUTPUT-SPACE ARRAY OF GRID POINTS. 

CONVERT OUTPUT-SPACE COORDINATES TC MAP COORDINATES. 

CONVERT MAP COORDINATES TO GEODETIC COORDINATES. 

CONVERT GEODETIC COORDINATES TO EARTH -CENTERED COORDINATES. 

CONVERT EARTH-CENTERED COORDINATES TO CORRECTED-TANGENT-SPACE COORDINATES. 

DO WHILE (# ITERATIONS < MAXIMUM) AND ( I NOMINAL-COMPUTED \ > MAXIMUM). 

CONVERT CORBECTED-TANGENT-SPACE COORDINATES TO UNCORRECTE D-T ANGENT SPACE. 
CCNVERT UNCO ERECTED -TANG ENT -SPACE COORDINATES TO INPUT SPACE. 

CONVERT INPUT-SPACE COORDINATES TC UNCORRECTED-TANGENT SPACE. 

CCNVERT UNCORRECTED-TANGENT-SPACE COORDINATES TO CORRECTED-TA NGE NT SPACE. 
COMPUTE DIFFERENCES OF NOMINAL AND COMPUTED TANGENT-SPACE COORDINATES. 

ENDDO. 

PRINT RESULTS. 

PRINT PROGRAM-EXIT MESSAGE. 

****** END OF SPECIFICATION ****** 
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S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


'NAME INITOS 

DESIGNER STEPHEN H. MURPHREY 

FUNCTION INITIALIZE A RECTANGULAR LATTICE OF OUTPUT-SPACE GRID-POINT 

CCORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE . 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. 

» m * 

CALL 

INITOS 

{ NUM 

, HO , VO ) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

NUM 

I*« 

I 

— 

NUMBER 

OF POINTS TO TRANSFORM 

- 

HO 

R*8 

I 

NUM 

OUTPUT 

SPACE COORDINATE 

PIXELS 

VO 

R*8 

I 

NUM 

OUTPUT 

SPACE COORDINATE 

PIXELS 


PDL (PROGRAM OUTLINE) 

TO BE DETERMINED 

****** IND OF SPECIFICATION ****** 
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4 ***** 


SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME 


GIMAP 


DESIGNER STEPHEN W. MURPHREY 

FUNCTION CONVERT OUTPUT-SPACE PIXEL COORDINATES TO (UTH) MAP 

COORDINATES. 


INPUTS 

OUTPUTS 

COMMON BLOCKS .. 
RESTRICTIONS ... 


1. SEE CALLING SEQUENCE. 

1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


CALLING SEQ. 

* • • 

CALL 

OIMAP 

( NUM , CNORTH , CEAST , BETA 

DNORTH , DEAST , HO r VO 

f 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

NUM 

T* i) 

T 

- 

NUMBER CF POINTS TO TRANSFORM 

— 

BETA 

R*8 

i 

— 

HEADING ANGLE 

RADIANS 

CNORTH 

R*8 

i 

NUM 

UTM NORTHING OF FORMAT CENTER 

METERS 

CEAST 

?*8 

T 

NUM 

UTM EASTING CF FORMAT CENTER 

METERS 

DNORTH 

R*8 

I 

NUM 

UTM NORTHING COORDINATE 

METERS 

DEAST 

R*8 

I 

NUM 

UTM EASTING COORDINATE 

METERS 

HO 

K*8 

I 

NUM 

OUTPUT SPACE COORDINATE 

PIXELS 

VO 

F*8 

T 

-i- 

NUM 

OUTPUT SPACE COORDINATE 

PI XELS 


PD I (PROGRAM OUTLINE) 


SB = SINE ( BETA ) 

CB = COSINEC BETA ) 

DO FOR EACH POINT 

HOC = HO (I) - HSZO / 2 

voc = vom - vszo / 2 

DNORTH! I) = CNORTH - HOC * HSPACI * CB 

DE AST (I) = CEAST + HOC * HSPACI * SB 

ENDDO 


VOC * VSPACI * SB 
VOC * VSPACI * CB 


****** 


END OF SPECIFICATION 


****** 
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SAF/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


'NAME MAPI! 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION CONVERTS (UTM) MAP COORDINATES TO GEODETIC LATITUDE AND 

LOGNITUDE COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 


RESTRICTIONS 

> * • • 





CALLING SEQ. 

• « • 

CALL 

MAPLL 

( AELLPS , BE LLPS , NUM , IZONED 
DNORTH , DEAST , IZONE ) 

, RDLATR , 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

AELLPS 

R*8 

I 


MAJOR AXIS OF EARTH ELLIPSOID 

METERS 

BELLPS 

R*8 

I 

— 

MINOR AXIS OF EARTH ELLIPSOID 

METERS 

NUM 

i*a 

I 

— 

NUMBER OF POINTS TO TRANSFORM 

- 

IZCNED 

1*4 

I 


UTM ZONE DESIRED 


EDLATR 

E*8 

0 

NUM 

GEODETIC LATITUDE OF INPUT POINTS RADIANS 

RDLONR 

F* 8 

0 

NUM 

LONGITUDE OF INPUT POINTS 

RADIANS 

DNOFTH 

p*8 

I 

NUM 

UTM NORTHING COORDINATE 

METERS 

DEAST 

R*8 

I 

NUM 

UTM EASTING COORDINATE 

METERS 

IZONE 

1*4 

I 

NUM 

UTM ZONE 

— 

COMMENTS ... 
r****** 

• • • • 

THIS MODULE ISA MODIFICATION OF A FORTRAN 
OETAINED FROM THE UNITED STATES GEOLOGICAL 

END OF SPECIFICATION 

PROGRAM 

SURVEY. 


RLONGR , 


****** 
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****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

^ AME LLEC 

DESIGNEE STEPHEN H. MURPHREY 

FUNCTION MAP POINTS EXPRESSED IN LATITUDE AND LONGITUDE TO EARTH- 

CENTERED CARTESIAN COORDINATES (X, Y, Z) . 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 


RESTRICTIONS ... 


CALLING SEO. 


CALL 

LLEC 

(N, RLDR, RLNR, ELEV, FISHA, FISHB, 

XEC, YEC 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRI PTIION 

UNITS 

N 

1*2 

T 

— 

# OF POINTS 

— 

RLNR 

R*8 

I 

N 

LONGITUDE OF POINTS 

RADIANS 

RLDR 

P*8 

I 

N 

LATITUDE OF POINTS 

RADIANS 

ELEV 

R*8 

I 

N 

TERRAIN ELEVATION OF POINTS 

METERS 

FISH A 

R*8 

I 

“ 

AXIS OF EARTH MODEL 
(FISHER ELLIPSOID) 

METERS 

FISHB 

R* 8 

I 

— 

n 

METERS 

XEC 

R* 8 

0 

N 

EARTH-CENTERED X-COORDINATE 
OF POINT 

METERS 

YEC 

R*8 

0 

N 

EARTH-CENTERED Y-COORDINATE 

METERS 

ZEC 

p*8 

0 

N 

EARTH-CENTERED Z-COORDINATE 

METERS 


PDL (PROGRAM OUTLINE) 


^Sc***** 


END OF SPECIFICATION 


****** 
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k ****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

'NAME ECCTP 

DESIGNER STEPHEN W. MORPHREY 

FUNCTION MAPS CONTROL POINTS EXPRESSED IN EARTH-CENTERED COORDINATES TO 

CORRECTED-TANGENT-SPACE COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS . . 
RESTRICTIONS ... 


CALLING SEQ 

• • • • 

CALL ECCTP 
ZECNAD, BETA! 

(NMCP. XEC, YEC, ZEC, XECNAD, YECNAD, 
I, RNADLD, RNADLN, CTPV, CTPH) 

VARIABLE 

TYPE 

I/O 

DIM 

DESCRTPTIION 

UNITS 

NMCP 

1*2 

I 

— 

# OF CONTROL POINTS 


XEC 

B*8 

I 

NMCP 

X-COORDINATE (EARTH-CENTERED) 

METERS 

YEC 

P * 8 

I 

NMCP 

Y-COORDINATE " " 

METERS 

ZEC 

R*8 

I 

NMCP 

Z -COORDINATE « » 

METERS 

XECNAD 

R* 8 

I 

— 

X-COORDINATE-NADIR " » 

METERS 

YEXNAD 

R*8 

I 

— 

Y-COORDINATE-NADIR " " 

METERS 

ZECNAD 

R*8 

T 

- 

Z -COORDINATE-NADIR " " 

METERS 

BETAN 

R* 8 

I 

— 

SPACECRAFT HEADING ANGLE 

RADIANS 

RNADLD 

R*8 

I 

— 

GEODETIC LATITUDE OF NADIR 

RADIANS 

RNADLN 

R*8 

I 

• 

LONGITUDE OE NADIR 

RADIANS 

CTPV 

R*8 

0 

NMCP 

C TP VERTICAL COORDINATE 

METERS 

CTPH 

R*8 

0 . 

NMCP 

CTP HORIZONTAL COORDINATE 

METERS 


r PDL (EBOGBAM OUTLINE) 


****** END OF SPECIFICATION ****** 
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****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME CTPOTP 

DESIGNER STEPHEN ¥J » MDRPHREY 

FUNCTION MAP POINTS FROM THE CORRECTED-TANGENT-SPACE TO THE UNCORRECTED- 

TANGENT-SPACE. THIS MAPPING CONSISTS OF APPLYING THE ATTITUDE 
AND ALTITUDE CORRECTIONS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL CTPUTP (NUM, TIMEPT, UTPV, UTPH, A A MOD , 

ALTNOM, APRXCH, CTPV, CTPH) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

NUM 

1*2 

I 


NUMBER OE POINTS TO TRANSFORM 

- 

TIMEPT 

F* 8 

I 

NUM 

TIME OF EACH POINT 

SECONDS 

UTPV 

R*8 

I 

NUM 

UTP VERTICAL COORDINATES 

METERS 

UTPH 

R* 8 

I 

NUM 

UTP HORIZONTAL COORDINATES 

METERS 

AA MOD 

F* 8 

I 

14 

A/A MODEL COEFFICIENTS 

— 

ALTNCM 

R*8 

I 

- 

NOMINAL S/C ALTITUDE 

METERS 

CTPV 

P* 8 

I 

NUM 

CTP VERTICAL COORDINATE 

METERS 

CTPH 

R*8 

I 

NUM 

CTP HORIZONTAL COORDINATE 

METERS 


PDI (PROGRAM OUTLINE) 

DO FOR EACH POINT 

EVALUATE ROLL, PITCH, YAH, AND DELTA ALT POLYNOMIALS AT TIMEPT (I) 


UTPV (I 

) = CTPV (I) 

i + YAW * CTPH (I) + PITCH * ALTNOM 

1 + DH * CTPH (I) / ALTNOM 

- ROLL * ( ALTNOM + CTPH(I)**2 / ALTNOM ) 


UTPH (I 
ENDDO 

) = CTPH (I) 


****** 


END OF SPECIFICATION 

****** 
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SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


'NAME UTPII 

DESIGNER STEPHEN W. MDSPHREY 

FUNCTION CONVERTS AN ARRAY OF UNCORRECTED-TANGENT- SPACE COORDINATES 

TO INPUT-SPACE COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS ........ 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. MSSCONS 
RESTRICTIONS . . . 


CALLING SEQ 

• • * • 

CALL 

UTPII 

(NMCP, VRA, VRC, VRAN, VRCN, DELVV, 



CLINI 

!, CSAMP, KWALTY, UTPV, UTPH, DTIME, CTIME) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPT I ION 

UNITS 

NM CP 

1*2 

T 

— 

# OF POINTS 

- 

VRA 

S*8 

I 

NMCP 

EAPTH ROTATION VELOCITY 

METERS/ 





(ALONG-TRACK) AT INPUT POINT 

SEC 

VRC 

F.*8 

I 

NMCP 

EARTH ROTATION VELOCITY 
(ACROSS-TRACK) AT INPUT POINT 
EARTH ROTATION VELOCITY 

ii 

VR AN 

R*8 

I 

- 

ii 





(ALONG-TRACK) AT NADIR 


VRCN 

R*8 

T 

— 

EARTH ROTATION VELOCITY 
(ACROSS-TRACK) AT NADIR 

ii 

DELVV 

R*8 

I 

— 

SPACECRAFT VELOCITY ERROR 

tt 

CLINE 

R*8 

I 

NMCP 

LINE OF CONTROL POINT 

PIXELS 

CSAMP 

R*8 

I 

NMCP 

SAMPLE OF CCNTRCI POINT 

PIXELS 

KWALTY 

1*2 

T 

NMCP 

STATUS OF POINTS 

PIXELS 

DTPV 

R* 8 

6 

NMCP 

VERTICAL COORDINATE OF POINT 
IN UTP 

METERS 

UTPH 

R* 8 

0 

NMCP 

HORIZONTAL COORDINATE OF 
POINT IN UTP 

METERS 

CTIME 

R*8 

O 

NMCP 

TIME (CONTINUOUS) ASSOCIATED 
WITH POINT 

SECONDS 

DTIME 

R*8 

O 

NMCP 

TIME ( DISCRETE) ASSOCIATED 

WITH POINT 

SECONDS 


PDL (PROGRAM OUTLINE) 

PRINT SUBROUTINE-ENTRY MESSAGE. 

PRINT HEADINGS. 

DO FOR EACH CONTROL EOINT : 

CORRECT HORIZ. COORDINATE FOR EARTH-CURVATURE / PANORAMIC-PROJECTION ERRORS. 
CORRECT FOE DISTORTIONS DUE TO EARTH ROTATION. 

CORRECT VERTICAL COORDINATE FOR SPACECRAFT-VELOCITY ERRORS. 

CORRECT VERTICAL COORDINATE FOR SCAN-SKEW ERRORS. 

CORRECT FOR DIFFERENTIAL INPUT SCALES. 

CORRECT HORIZONTAL COORDINATE FOR MIRROR-VELOCITY ERRORS. 

CONVERT TO IMAGE-CENTER-ORIGIN COORDINATES. 

PRINT RESULTS OF ABOVE ERROR CORRECTIONS. 

ENDDO 

PRINT SUBROUTINE-EXIT MESSAGE. 

****** END of SPECIFICATION ****** 
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SAF/LANDSAT data merging system Program specification 


NAME IIUTP 

DESIGNER STEPHEN W. MDRPHREY 

FUNCTION CONVERTS AN ARRAY OF INPUT-SPACE COORDINATES (LINE AND 

TO UNCORFECTED-TANGENT-SPACE COORDINATES. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. MSSCCNS 
RESTRICTIONS ... 


CALLING SEQ. 

• • • 

CALL 

CLINE 

IIUTP (NMCP, VRA, VRC, VRAN, VRCN, DELVV, 
!, CSAMP, KWALTY, UTPV, UTPH, DTIME, CTIME) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTIION 

UNITS 

NMCP 

1*2 

I 

- 

# OF POINTS 

_ 

VRA 

R*8 

I 

NMCP 

EARTH ROTATION VELOCITY 
(ALONG -TRACK) AT INPUT POINT 

METERS/ 

SEC 

VEC 

F* 8 

I 

NMCP 

EARTH ROTATION VELOCITY 
(ACROSS-TRACK) AT INPUT POINT 

ii 

VRAN 

R*8 

I 


EARTH ROTATION VELOCITY 
(ALONG-TRACK) AT NADIR 

n 

VP CN 

R* 8 

I 


EARTH ROTATION VELOCITY 
(ACROSS-TRACK) AT NADIR 

ii 

DELVV 

R* 8 

T 

- 

SPACECRAFT VELOCITY ERROR 

ii 

CLINE 

R*8 

I 

NMCP 

LINE OF CONTROL POINT 

PIXELS 

CSAMP 

R* 8 

I 

NMCP 

SAMPLE OF CONTROL POINT 

PIXELS 

KWALTY 

1*2 

I 

NMCP 

STATUS OF POINTS 

PIXELS 

UTPV 

R*8 

0 

NMCP 

VERTICAL COORDINATE OF POINT 
IN UTP 

METERS 

UTPH 

R*8 

0 

NMCP 

HORIZONTAL COORDINATE OF 
POINT IN UTP 

METERS 

CTIME 

R*8 

0 

NMCP 

TIME (CONTINUOUS) ASSOCIATED 
WITH POINT 

SECONDS 

DTIME 

R*8 

0 

NMCP 

TIME (DISCRETE) ASSOCIATED 

WITH POINT 

SECONDS 


PDL (PROGRAM OUTLINE) 

PRINT SUBROUTINE-ENTPY MESSAGE. 

PRINT HEADINGS. 

DO FOR EACH CONTROL POINT: 

CONVERT TO I MAGE -CENTER -ORIGIN COORDINATES. 

CORRECT HORIZONTAL COORDINATE FOR MIRROR-VELOCITY ERRORS. 

CORRECT FOR DIFFERENTIAL INPUT SCALES. 

CORRECT VERTICAL COORDINATE FOR SCAN-SKEW ERRORS. 

CORRECT VERTICAL COORDINATE FOR SPACECRAFT-VELOCITY ERRORS. 

CORRECT FOR DISTORTIONS DUE TO EARTH ROTATION. 

CORRECT HORIZ. COORDINATE FOR EARTH-CURVATURE / PANORAMIC- PROJECTION 
PRINT RESULTS OF ABOVE ERROR CORRECTIONS. 

ENDDO 

PRINT SUBFOUTINE-EXI 2 MESSAGE. 

****** END OF SPECIFICATION 


SAMPLE) 


ERRORS. 


****** 
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ORIGINAL PAGE $ 
OF POOR QUALITY 


****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME UTPCTP 

DESTGNEP STEPHEN W. MURPHREY 

FUNCTION MAP POINTS FROM THE UNCORRECTED-TANGENT- SPACE TO THE CORRECTED- 

TANGENT-SPACE . THIS MAPPING CONSISTS OF APPLYING THE ATTITUDE 
AND ALTITUDE CORRECTIONS. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL UTPCTP (NUM , TIMEPT. UTPV, UTFH, AAMOD, 

ALTNOM, APRXCH, CTPV, CTPH) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTION 

UNITS 

NUM 

1*2 

I 

— 

NUMBER OF POINTS TO TRANSFORM 

- 

TIMEPT 

R*8 

I 

NUM 

TIME OF EACH POINT 

SECONDS 

UTPV 

F* 8 

I 

NUM 

UTP VERTICAL COORDINATES 

METERS 

UTPH 

R*8 

I 

NUM 

UTP HORIZONTAL COORDINATES 

METERS 

AAMOD 

R*8 

I 

14 

A/A MODEL COEFFICIENTS 

- 

ALTNOM 

R*8 

I 

— 

NOMINAL S/C ALTITUDE 

METERS 

APRXCH 

R*8 

I 

NUM 

APPROXIMATE CTP HORIZONTAL 

METERS 

CTPV 

F*8 

I 

NUM 

CTP VERTICAL COORDINATE 

METERS 

CTPH 

R*8 

I 

NUM 

CTP HORIZONTAL COORDINATE 

METERS 


PDL (PROGRAM OUTLINE) 

DO FOR EACH POINT 

COMPUTE THE TIME FOP THE POINT 

COMPUTE ROLL, PITCH. YAW. AND DELTA ALT BY EVALUATING THE A/A 
POLYNOMIALS AT THE TIME OF THE POINT 
COMPUTE THE DELTA CORRECTIONS FOR THE TANGENT PLANE COORDINATES 
BY PUTTING THE ATTITUDE VALUES INTO THE DISTORTION MODELS 
ADD THE DELTA CORRECTIONS TO THE UNCORRECTED VALUES 
FNDDO 

****** SND OF SPECIFICATION ****** 
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SAR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 


NAME GTERAS 

DESIGNER STEPHEN W. MURPHREY 

FUNCTION ESTIMATES THE ACCURACY TO WHICH GEODETIC COORDINATES (LAT/LON) 

CAN BE ASSIGNED TO POINTS IN A SCENE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. MSSCCNS 
RESTRICTIONS ... 

CALLING SEO. ... CALL GTERAS (NMCP, RACOV, PYCOV, CTPV, CTPH, SIGH, AESTV,AESTR) 


VARIABLE 

TYPE 

I/O 

DIM 

DESCRIPTIION 

UNITS 

NMCP 

1*2 

I 

— 

# OF CCNTBOL POINTS 

_ 

RACOV 

R*8 

I 

6X6 

COVARIANCE MATRIX (ROLL AND 
ALTITUDE DEVIATION) 

— 

PYCOV 

R* 8 

T 

8X8 

COVARIANCE MATRIX (PITCH AND YAW) 

- 

CTPV 

R*8 

I 

NMCP 

VERTICAL COORDINATE OF POINT 
IN CORRECTED TANGENT PLANE 

METERS 

CTPH 

R*8 

I 

NMCP 

HORIZONTAL COORDINATE OF POINT 
IN CORRECTED TANGENT PLANE 

METERS 

AESTV 

R*8 

0 

NMCP 

ACCURACY ESTIMATE OF VERTICAL 
COORDINATE 

METERS 

AESTH 

R*8 

0 

NMCP 

~ ACCURACY ESTIMATE OF HORIZONTAL 
COORDINATE 

METERS 


PDL (PROGRAM OUTLINE) 


DO FOR EACH INPUT POINT 

COMPUTE THE TIME OE THE POINT 

COMPUTE MATRIX Ml WHICH RELATES THE A/A MODEL COVARIANCES TO 
HORIZONTAL GEODETIC COVARIANCES 
COMPUTE MATRIX M2 WHICH RELATES THE A/A MODEL COVARIANCES TO 
VERTICAL GEODETIC COVARIANCES 


FROM THE TIME OF THE POINT, COMPUTE THE MATRIX T1 WHICH ALLOWS Ml 
TO BE EVALUATED AT THE TIME OF THE GIVEN POINT 
FROM THE TIME OF THE POINT, COMPUTE THE MATRIX T2 WHICH ALLOWS M2 
TO BE EVALUATED AT THE TIME OF THE GIVEN POINT 
COMPUTE FI = M1*T1 WHICH PELATES THE A/A MODEL COVARIANCES TO 
HORIZONTAL GEODETIC COVARIANCES AT THE TIME OF THE POINT 
COMPUTE F2 = M2*T2 WHICH RELATES THE A/A MODEL COVARIANCES TO 
VERTICAL GEODETIC COVARIANCES AT THE TIME OF THE POINT 
COMPUTE THE VARIANCE OF THE GEODETIC ACCURACY OF THE HORIZONTAL 
COORDINATE OF THE GIVEN POINT (H VAR = F2 *PYCOV*F2T) 

COMPUTE THE variance o? the redtie'ptc arnn-Rarv n e the t. rnm 


COORDINATE OF THE GIVEN POINT (H VAR = F2 *PYCOV*F2T) 

COMPUTE THE VARIANCE OF THE GEODETIC ACCURACY OF THE VERTICAL COOR- 
DINATE OF THE GIVEN POINT (V VAR = FI *RACOV*FlT) 

COMPUTE THE HORIZONTAL "ERROR ELLIPSE RADII"= SIGM*SORT(H VAR) 
COMPUTE THE VERTICAL "ERROR ELLIPSE RADII"= SIGM* SORT (H VAR) 

ENDDO ' 


****** 


END OF SPECIFICATION 


****** 
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ORIGINAL PAGE lb 
OF POOR QUALITY 


****** SAR/LANDSAT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME GTWANC 

DESIGNER ....... STEPHEN W. MORPHREY 

FUNCTION WRITES GEOMETRIC-TRANSFORMATION PARAMETERS TO THE SLDMS 

ANCILLARY DATA SET. 

INPUTS 1. SEE CALLING SEQUENCE, 

OUTPUTS 1. SEE CALLING SEQUENCE. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... CALL GTWANC () 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

PRINT PROGRAM-ENTRY MESSAGE. 

WRITE ERROR MODELS TO THE SLDMS ANCILLARY DATA SET. 

WRITE ERROR ASSESSMENT TO THE SLDMS ANCILLARY DATA SET. 

WRITE OUTPUT-SPACE DEFINITION TO THE SLDMS ANCILLARY DATA SET. 
WRITE GRID-POINT CORRESPONDENCE TO THE SLDMS ANCILLARY DATA SET. 
PRINT PROGRAM-EXIT MESSAGE. 

****** END OF SPECIFICATION 


***** 


*****3 
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B . 7 RESAMPLING PROGRAM 


B.7.1 Statement of Problem 

Resampling, as defined here, is the creation of a two-dimensional array of 
eight-bit words (called output-space pixels) from another two-dimensional array 
of eight-bit words (called input-space pixels) by the process indicated in Fig- 
ure B-8. It involves two main steps: the creation of an intermediate two- 
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Figure B-8. Horizontal and Vertical Resampling 


dimensional array (called hybrid-space pixels) directly from the input-space 
array; and the creation of the output-space array from the hybrid-space array. 
The first of these steps is called horizontal resampling, and the second is 
called vertical resampling. 


B.7.2 Data Flow 

The three resampling spaces are shown in Figure B-9. Hybrid space is created by 
a one-dimensional resampling of input space. This is a horizontal resampling. 
Output space is created by a one-dimensional resampling of hybrid space. This 
is vertical resampling. 
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INPUT SPACE HYBRID SPACE OUTPUT SPACE 


Figure B-9 . Resampling Spaces 


In each, case, the resultant image is created one line at a time. The flow is 
from the top of the image to the bottom. A resampled line is created in seg- 
ments, as shown in Figure B-10. 

In order to resample one line's worth of data, the location of each point in 
the input data array must be obtained. This is done, within each segment, by 
a linear interpolation scheme. Consider Figures B-10 and B-ll. A rectangular 
lattice of grid points (usually, but not necessarily, equally spaced) is set up 
in the coordinate system of the image being created. These interpolation grid 
points are then mapped to the corresponding input-data coordinate system. The 
mapped grid points are the SVO and LHO points shown in Figures B-10 and B-ll. 

In each case, another set of grid points with a finer mesh is created by linear 
interpolation between the mapped grid points . These are the SV and LH points 
shown in Figures B-10 and B-ll. There is a set of SV grid points for each line 
in hybrid space, and there is a set of HV grid points for each column in output 
space. In the case of horizontal resampling, some high-frequency geometric 
errors that are a function of the line number are corrected by adjusting the SV 
grid points by the magnitude of the errors. The resulting points are the SV' 
grid points shown in Figure B-10. Finally, the location of each point in the 
input data array is calculated by linear interpolation between the SV' or LH 
grid points . 

B.7.3 Inputs to Resampling Process 

The following data is required to perform the resampling function: 

a. Hybrid-space grid-point locations -and corresponding input-space grid- 
point locations (SVO grid points). 
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Figure B-10. Construction of Hybrid Space 
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b. Constants required to compute horizontal, high-frequency correction 
coefficients. 

c. Output-space grid-point locations and corresponding hybrid-space 
grid-point locations (LHO grid points). 

d. Input-space array of data values (assume 8-bit pixels). 

e. Constants required to evaluate resampling polynomial. 

B.7.4 Algorithm Considerations 

The general formula for a six-point, one-dimensional resampling algorithm may 
be stated as follows : 



where 

Ijj = the intensity of pixel N 

W_ w = the weight for pixel N in the term of degree I (a known constant) 
K = a known constant 

D = the (positive) distance between pixel 0 and the output pixel. 

The relationships are illustrated in figure B-12. 
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Figure B-12. Six-point Resampling 


If. equation (1) is implemented in hardware or software in such a way that the 
weights Wj ^ and the constant K may be easily changed, then a large class of 
resampling ’algorithms” are included in the single implementation. In particular, 
all six-point and four-point cubic -convolution resamplers are included, as is 
the four-point quadratic-convolution algorithm. Some common resampling-algorithm 
constants are given in Table B-24. 

It is obvious from equation (1) that, the single-implementation approach would 
be a very poor strategy for a software resampler on an IBM 370 computer. Using 
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Table B-24. Typical Resampling Constants (Weights) 
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equation (1) for a four-point quadratic convolution would use more than twice 
as much CPU time as would he required by the equation 


V = D(D(I_ a - I p - I r + I 2 ) + C-l _ 1 + ~ I 2 )) + I 0 , C2) 

which would eliminate all of the unnecessary multiplications by zero. A good 
software approach is to code each resampling algorithm separately to minimize CPU 
utilization. This approach will be used for the SLDMS Resampling Program. The 
classic 4-point cubic convolution algorithm will be used. 


B.7.5 Handling of Edges in Resampling Via Hybrid Space 

The boundary of a digital image that has been resampled is usually not a rec- 
tangle that is aligned with the rows and columns of pixels. There is normally 
some non-image fill data at the edges of a digital image. In order to simplify 
the handling of edges during resampling, the input-space data will be artificially 
enlarged by placing fill pixels around the edges. This padding will be suffi- 
ciently large to guarantee- that input-space data exists for every hybrid-space 
pixel and that hybrid-space data exists for every output-space pixel. 

Figures B-10 and B-ll assume no such enlargement of input space. Figures B-13 
and B-14 assume the enlargement described above. In this case, there are no SV' 
points (or LH points) that are not at least three columns (or rows) interior to 
the enlarged input space (or enlarged hybrid space). Therefore, no special case 
will have to be made for output-space pixels that would be mapped outside the 
input-space image data. All hybrid-space and output-space pixels are created m 
the same way by the resampling function. A scale drawing of the three resampling 
spaces is shown in Figure B-15. 


B.7.6 Program Description 

Although the resampling problem is conceptually straightforward, computer software 
that performs resampling is somewhat complicated. This program, the largest CPU 
user in the SLDMS, will be written primarily in assembler language to minimize 
computer utilization costs. 

The design of this program is described in Figure B-16 and the program specifi- 
cations that follow. These describe the cubic convolution algorithm. The 
nearest-neighbor resampling algorithm will use the same design. 

The design of the resampling program has a few significant features. The resamp- 
ling is performed via an intermediate hybrid space. This method was chosen to 
minimize computer utilization. The four point resampling algorithm will he imple- 
mented in fixed-point arithmetic to minimize computer use. The alternative, 
floating-point arithmetic, would use considerably more CPU. 

Another significant feature is the large hybrid-space buffer. In order to mini- 
mize the I/O processing of the resampling program, the output image will be 
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Figure B-15. Resampling Spaces (Drawn to Scale) 






created in full-line increments. This requires hybrid-space data from several 
lines to be available at essentially the same time. This availability will be 
provided by a large buffer in the computer memory. The size of this hybrid- 
space buffer is dependent on the size of the output space and the particular 
geometric transformation between the input and output spaces. The buffer will 
be allocated dynamically during program execution. This will permit the program 
to be run in small regions where possible. 

Due to the generally high running cost of the resampling program, it will have 
a checkpoint/restart capability. A checkpoint is taken for each row of LHO grid 
points. In the event of a computer failure during execution of the resampling 
program, processing can be resumed at the point at which the last checkpoint 
was taken. 
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FIGURE B-16. MODULE HIERARCHY FOR RESAMPLING PROGRAM. 
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****** 


****** 


S AS/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 


NAME RSPL 

DESIGNEE STEPHEN W. MUBPHREY 

FUNCTION THIS IS THE TOP-LEVEL QRIVEB FOR THE RESAMPLING PROGRAM. 


INPUTS 1. AN SLDMS ANCILLARY DATA SET FOR THE SCENE TO BE PROCESSED. 

2. AN IMAGE DATA SET IN MIST FORMAT. 

3. A DATA SET CONTAINING USEE'S REQUESTS (ON CARD IMAGES). 

OUTPUTS 1. AN IMAGE DATA SET IN MIST FORMAT. 

2. A PRINTER LISTING DATA SET. 


COMMON BLOCKS .. 
RESTRICTIONS ... 


CALLING SEQ. ... NOT APPLICABLE. 

PDL (PROGRAM OUTLINE) 

PRINT PROGRAM- ENTRY MESSAGE 
START CPU TIMING 

IF THIS IS A CHECKPOINT/RESTART RUN THEN 
PERFORM RESTART INITIALIZATION 
POSITION DATA SETS 
PRINT RESULTS 

ELSE 


0 ^«SSL* 


READ GRID-POINT CORRESPONDENCE FROM ANCILLARY DATA SET 
INITIALIZE CONSTANTS AND OTHER DATA AREAS 
PRINT RESULTS 

OBTAIN HYBRID SPACE BUFFER AREA IN (VIRTUAL) MEMORY 
CPEN DATA SETS 
PRINT RESULTS 
ENDIF 

TAKE CPU-TIMING READING 
PERFORM VERTICAL RESAMPLING 
TAKE CPU-TIMING READING 
CLOSE DATA SETS 
PRINT RESULTS OF CPU TIMINGS 
PRINT PROGRAM-EXIT MESSAGE 


* * ** ** 


END OF SPECIFICATION 


****** 
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****** SAR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME RSPLV 

DESIGNER STEPHEN W . MURPHREY 

FUNCTION PERFORMS VERTICAL RESAMPLING (INCLUDES CONTROL OF HORIZONTAL 

RESAMPLING) . 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

DO FROM J = 1 TO N - 1 <N= NUMBER OF RCWS OF LHO GRID POINTS ) 

COMPUTE LHK GRID POINTS FOR EACH COLUMN (K) OF PIXELS 
COMPUTE DUK FOR EACH COLUMN OF PIXELS 

DO FROM I " C TO MJ - 1 ( MJ = # ROWS OF PIXELS IN SEGMENT J ) 

DETERMINE LOCATION OF FIRST PIXEL IN OUTPUT-SPACE BUFFER 
DO FROM K = 1 TO L ( L = WIDTH OF ONE OUTPUT-SPACE ROW ) 

C = LHK + I * DUK 

IC = G(C) + 2 ( G(C) = MAXIMUM INTEGER <= C \ 

DO WHILE LINE ID IS NOT IN THE HYBRID-SPACE BUFFER 

HORIZONTALLY RESAMPLE ONE LINE OF INPUT-SPACE DATA 
ENDDO 

OBTAIN THE FOUR PIXEL INTENSITIES STARTING WITH PIXEL IC 
£ _ C _ XC + 2 

PERFORM THE FOUR-POINT RESAMPLING (EQUATION 2 ) 

INSERT RESULT INTO THE OUTPUT-SPACE BUFFER 
DETERMINE LOCATION OF NEXT PIXEL IN OUTPUT-SPACE BUFFER 
ENDDO 

WRITE ONE OUTPUT-SPACE RECORD 
ENDDO 

TAKE A CHECKPOINT 
ENFDO 

****** END OF SPECIFICATION ****** 
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ORIGINAL PAGE IS 
OP POOR QUALITY 


****** S AH/LANDS AT BATA MERGING SYSTEM PROGRAM SPECIFICATION ****** 

NAME RSPLH 

DESIGNER ....... STEPHEN W. MURPHREY 

FUNCTION HORIZONTAL RESAMPLING OF ONE LINE OF INPUT-SPACE DATA. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. ... • 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM CUT1INE) 

READ NEXT RECORD INTC INPUT-SPACE BUFFER 

READ NEXT SET OF SV' GRID-POINT COORDINATES 

DETERMINE LOCATION OF NEXT LINE IN HYBRID-SPACE BUFFER 

DO FROM J = 1 TO N - 1 (N= NUMBER OF SV' GFID POINTS IN ONE ROW ) 

DETERMINE DVJ ' 

DO FROM I = 0 TO MJ - 1 ( MJ = NUMBER OF PIXELS IN SEGMENT J ) 

C = SV'J + I * DVJ 

IC = G (C) - 1 ( G(C) = MAXIMUM INTEGER <= C ) 

OBTAIN THE FOUR PIXEL INTENSITIES STARTING WITH PIXEL IC 
D = C - IC - 1 

PERFORM THE FOUR-POINT RESAMPLING (EQUATION' 2) 

INSERT THE RESULT INTO THE HYBRID-SPACE BUFFER 

DETERMINE THE LOCATION OF THE NEXT PIXEL IN HYBRID-SPACE BUFFER 
ENDDO 
ENDDO 

****** END OF SPECIFICATION ****** 
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****** S AR/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME MRFTXT 

DESIGNER STEPHEN W. MORPHPEY 

FUNCTION PROCESSES X-FORMAT IMAGE DATA. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS . . 

RESTRICTIONS ... 

CALLING SEQ. ... CALL MRFTXI < ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 


IF 


THIS IS NOT A THERMAL-BAND RECORD THEN 
DO FOR EACH 8-BYTE GROUP OF INPUT DATA 


MOVE BYTES 1 
MOVE EYTES 3 
MOVE BYTES 5 
MOVE BYTES 7 
ENDDO 

IF THIS IS INPUT DATA FOR 


AND 

AND 

AND 


TO BAND-1 AREA OF OUTPUT BUFFER 
TO BAND -2 AREA OF OUTPUT BUFFER 
TO BAND-3 AREA OF OUTPUT BUFFER 


AND 8 TO BAND-4 AREA OF OUTPUT BUFFER 


STRIP 4 


THEN 


D « ILC / { ILA - LLC - 6 ) = LENGTH OF SEGMENT W/O DUPLICATED 
' ”C + 6 ) / LLA = LINE-] 

INITIALIZE "FROM* 


LENGTH FACTOR (SAVED IN OUTPUT 


LLF = { LLd + 6 ) / LLA 

a{ID „ T01I aDI)BESSE 3 

DO FOP EANDS 1 TO 4 
R = LLA 

DO WHILE R > D 
MOVE D PIXELS 
INCREMENT "FROM" BY D + 1 
INCREMENT "TO" BY D 
R = E - ( D + 1 ) 

ENDDO 

IF R > 0 

MOVE R PIXELS 
ENDIF 

FILL REMAINDER OF BUFFER AREA WITH FILL CHARACTERS 
EC FOR EACH IMAGE PIXEL 

OBTAIN V = PIXEL VALUE 

INCREMENT V COUNTER FOR CURRENT DETECTOR BY 1 
ENDDO 
ENDDO 
ENDIF 


ELSE 

MOVE ALL IMAGE DATA TO BAND-5 AREA OF OUTPUT BUFFER 
IF THIS IS INPUT DATA SET FOR STRIP 4 
REMOVE LINE-LENGTH CORRECTION PIXELS 
FILL REMAINDER OF BUEFEP AREA WITH FILL CHARACTERS 
DO FOR EACH IMAGE PIXEL 

OBTAIN V = PIXEL VALUE 

INCREMENT V COUNTER FOR CURRENT DETECTOR BY 1 
ENDDO 
ENDIF 
ENDIF 


****** 


END OF SPECIFICATION 


****** 


PIXEL 

BUFFER) 


****** 
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****** SAB/LANDS AT DATA MERGING SYSTEM PROGRAM SPECIFICATION 

NAME MRFTXS 

DESIGNER STEPHEN S. MUEPHREY 

FUNCTION READS AND PROCESSES SIAT FILE FROM X-FORMAT MSS TAPE. 

INPUTS 1. SEE CALLING SEQUENCE. 

OUTPUTS 1. SEE CALLING SEQUENCE. 


COMMON BLOCKS .. 

RESTRICTIONS ... 

CALLING SEQ. . . . CALL MRFIXS ( ) 

VARIABLE TYPE I/O DIM DESCRIPTION UNITS 


PDL (PROGRAM OUTLINE) 

DETERMINE FORMAT LEVEL 
READ ALL 7 OR 8 SIAT RECORDS 
IF FORMAT LEVEL 2 OH 3 THEN 
PRINT CALIBRATION MODIFIERS 
ENDIF 

PRINT SPACECRAFT PERFORMANCE DATA 

PRINT ANNOTATION BLOCK DATA 

PRINT MSS COMPUTATIONAL DATA 

PRINT IMAGE LOCATION DATA 

SPITE SIAT DATA TO ANCILLARY DATA SET 

******* END OF SPECIFICATION 


y 
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